Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,703 A→G 7.9% intergenic (+22/‑206) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7030AG 7.9% 229.7 / 8.7 76intergenic (+22/‑206)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (4/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTA  >  NC_000913/3334558‑3334846
                                                                                                                                                 |                                                                                                                                                
gACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑g                                                                                                                                              >  4:199394/1‑149 (MQ=255)
   gcgTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                           >  1:169765/1‑149 (MQ=255)
        aGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGAt                                                                                                                                      >  4:692147/1‑149 (MQ=255)
         gCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATg                                                                                                                                     >  4:647235/1‑149 (MQ=255)
          cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGC‑GCTGCGCGGGATga                                                                                                                                    <  4:203016/149‑2 (MQ=255)
           gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGt                                                                                                                                   <  3:692147/149‑1 (MQ=255)
                      gaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGtt                                                                                                                        >  2:157713/1‑149 (MQ=255)
                       aagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTc                                                                                                                       <  2:169765/149‑1 (MQ=255)
                        agCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTcc                                                                                                                      <  4:436945/149‑1 (MQ=255)
                         gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                     <  3:28599/149‑1 (MQ=255)
                         gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                     <  3:657190/149‑1 (MQ=255)
                         gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCCCACATCGCTGTTCAACGCGATCGTTAATCCCCCCCCCTCTTCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                     <  1:38311/149‑1 (MQ=255)
                           aGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                     <  3:614607/147‑1 (MQ=255)
                           aGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                     >  4:614607/1‑147 (MQ=255)
                             aCAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                 >  2:731587/1‑149 (MQ=255)
                                     aTCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                         >  3:543478/1‑149 (MQ=255)
                                       cGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                       <  1:419895/149‑1 (MQ=255)
                                        gCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCggg                                                                                                                                        <  1:518035/115‑1 (MQ=255)
                                        gCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCggg                                                                                                                                        >  2:518035/1‑115 (MQ=255)
                                        gCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTc                                                                                                      <  1:165039/149‑1 (MQ=255)
                                          aGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                                                    <  1:157713/149‑1 (MQ=255)
                                              ttACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctct                                                                                                <  4:543478/149‑1 (MQ=255)
                                                  aGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                            <  3:647235/149‑1 (MQ=255)
                                                    gTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                                              >  4:404192/1‑145 (MQ=255)
                                                    gTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                                              <  3:404192/145‑1 (MQ=255)
                                                          ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaa                                                                                    <  3:386893/149‑1 (MQ=255)
                                                          ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaa                                                                                    >  4:657319/1‑149 (MQ=255)
                                                                acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                              <  3:199394/149‑1 (MQ=255)
                                                                  cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAggg                                                                            >  4:201720/1‑149 (MQ=255)
                                                                   ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAac                                                                                           <  3:106091/133‑1 (MQ=255)
                                                                   ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAac                                                                                           >  4:106091/1‑133 (MQ=255)
                                                                   ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                           >  1:243919/1‑149 (MQ=255)
                                                                     ttGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                                         >  4:49309/1‑149 (MQ=255)
                                                                      tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAACATAGGGAtgt                                                                        <  3:49309/149‑1 (MQ=255)
                                                                           agaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                              >  4:127286/1‑90 (MQ=255)
                                                                           agaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                              <  3:127286/90‑1 (MQ=255)
                                                                                 aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtata                                                             >  1:506068/1‑149 (MQ=255)
                                                                                 aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtata                                                             >  1:98867/1‑149 (MQ=255)
                                                                                 aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtata                                                             >  3:461311/1‑149 (MQ=255)
                                                                                   tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                           >  2:545401/1‑149 (MQ=255)
                                                                                   tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                           >  2:377451/1‑149 (MQ=255)
                                                                                     aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                         >  3:569136/1‑149 (MQ=255)
                                                                                      tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                        <  2:302589/149‑1 (MQ=255)
                                                                                      tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                        <  4:569136/149‑1 (MQ=255)
                                                                                      tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                        >  3:391506/1‑149 (MQ=255)
                                                                                      tCGGCCTCGCGAACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                        >  3:455945/1‑149 (MQ=255)
                                                                                       cGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTGATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGGAAAATAGGGATGTGCATCTTTATATATTCt                                                       >  1:526788/1‑149 (MQ=255)
                                                                                        ggCCTAGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                      >  2:496617/1‑149 (MQ=255)
                                                                                         gCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGa                                                     <  1:731587/149‑1 (MQ=255)
                                                                                          ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                    <  3:201720/149‑1 (MQ=255)
                                                                                                gcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAaca                                                                                                             <  2:404142/86‑1 (MQ=255)
                                                                                                gcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAaca                                                                                                             >  1:404142/1‑86 (MQ=255)
                                                                                                 cacaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                             <  4:461311/149‑1 (MQ=255)
                                                                                                  acaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcaca                                            >  3:373363/1‑149 (MQ=255)
                                                                                                    aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACAct                                          >  3:637843/1‑149 (MQ=255)
                                                                                                        cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                      <  2:506068/149‑1 (MQ=255)
                                                                                                        cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                      <  1:496617/149‑1 (MQ=255)
                                                                                                           ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACa                                   >  1:104303/1‑149 (MQ=255)
                                                                                                           ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACa                                   <  4:391506/149‑1 (MQ=255)
                                                                                                                 gcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACt                             >  1:97250/1‑149 (MQ=255)
                                                                                                                   gATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATGTTCACACTCTTTAc                                    <  2:410864/140‑1 (MQ=255)
                                                                                                                   gATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATGTTCACACTCTTTAc                                    >  1:410864/1‑140 (MQ=255)
                                                                                                                   gATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttt                           >  4:558079/1‑149 (MQ=255)
                                                                                                                     tCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttt                           <  2:725020/147‑1 (MQ=255)
                                                                                                                     tCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttt                           >  1:725020/1‑147 (MQ=255)
                                                                                                                        ttAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTaga                      <  4:637843/149‑1 (MQ=255)
                                                                                                                               cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg               <  2:371771/149‑1 (MQ=255)
                                                                                                                                ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg              >  2:422861/1‑149 (MQ=255)
                                                                                                                                  cccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGcc            >  3:88740/1‑149 (MQ=255)
                                                                                                                                     accATCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                        >  1:152475/3‑134 (MQ=255)
                                                                                                                                     accATCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                        <  2:152475/132‑1 (MQ=255)
                                                                                                                                      ccATCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGAGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa        >  4:11766/2‑149 (MQ=255)
                                                                                                                                        aTCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATCTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAga      >  4:397086/1‑148 (MQ=255)
                                                                                                                                        aTCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg      >  2:353724/1‑149 (MQ=255)
                                                                                                                                          cccGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAg    <  3:657319/149‑1 (MQ=255)
                                                                                                                                            cgcgCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAgaa  <  1:377451/149‑3 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTA  >  NC_000913/3334558‑3334846

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑