Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 12.0% intergenic (+46/+56) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950CG12.0% 1770.0 / 93.8 576intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (299/200);  new base G (9/68);  total (308/268)
Fisher's exact test for biased strand distribution p-value = 2.78e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 585 total aligned reads displayed.

CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608546‑3608842
                                                                                                                                                     |                                                                                                                                                   
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                    >  2:7488729/1‑150 (MQ=255)
   cGCTGTTTGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                  <  1:966475/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                 <  1:5508283/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                 <  2:6545026/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                 <  1:1151664/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                 >  1:2232268/1‑150 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                             <  1:9146032/150‑1 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                             <  1:2972112/150‑1 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                             <  2:8793365/150‑1 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTAGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                             <  1:4408515/150‑1 (MQ=255)
        gtAGGGATGCCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGGGACAGCGAATGCGTTAAGATTGTTGCTCAGGAGATAAGGCAAACCTATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                             <  2:3152116/149‑1 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                            <  2:7371041/150‑1 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                            >  1:4802880/1‑150 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                           <  1:7126877/150‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGAGTGTTGCGCAGGAGAGAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                           >  2:10466694/1‑150 (MQ=255)
          aGGGATGACCGGGTTGGGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                           <  2:3254069/150‑1 (MQ=255)
           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                          <  2:8186057/150‑1 (MQ=255)
           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                          <  2:3946993/150‑1 (MQ=255)
            ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                         >  1:10870045/1‑150 (MQ=255)
                 aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                    <  1:2524760/150‑1 (MQ=255)
                 aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                    >  2:6290483/1‑150 (MQ=255)
                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                   >  2:2270697/1‑150 (MQ=255)
                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                            >  2:1772779/1‑140 (MQ=255)
                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                  >  1:5733965/1‑150 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                <  1:9959387/150‑1 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                >  1:2847163/1‑150 (MQ=255)
                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                               <  2:727880/150‑1 (MQ=255)
                       ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                              <  2:10579421/150‑1 (MQ=255)
                          tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                           >  1:4244970/1‑150 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                         <  1:10669526/150‑1 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                         >  1:10303478/1‑150 (MQ=255)
                              tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                           <  2:598267/146‑1 (MQ=255)
                               ggcAGTGCTGGCAGATAGGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                      <  2:8908509/150‑1 (MQ=255)
                               ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTACCCTCTCCCCAGAGGGGCGAggggac                                                                                                                                       >  1:3235813/1‑130 (MQ=255)
                                  aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                                                   >  1:731764/1‑150 (MQ=255)
                                  aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCACCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                                                   >  2:9395687/1‑150 (MQ=255)
                                    tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                 >  2:5546861/1‑150 (MQ=255)
                                    tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                 >  1:10529411/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAAGGCGGCATCGa                                                                                                              >  2:6441355/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                              >  1:3952880/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                              >  1:10006556/1‑150 (MQ=255)
                                         gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                     <  2:5850165/141‑1 (MQ=255)
                                                  ggggCGCCGGTGCTGGTGACCGCTAATGCTTTATGTTTGTTGGGCTGGTGATACGCAAACCCGATCGCAACTTTGCGCGCGATGGGTCCCCTGCCCCCTGGGGGGTGTGTTTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGc                                                                                                   <  2:10280168/150‑1 (MQ=255)
                                                  ggggCGCCGGTGCTGGTGACCGCGAAGGCTTTATCTTTGTTGCGCCGGAGATACGGCACCCCGAGCGAAACTTGGCGCGCGATCGGTCCCCTCGCCCCGCGGGGGAGTGGTTTAGGGGGCGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                   <  2:9745065/150‑1 (MQ=255)
                                                    ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                 >  2:5055230/1‑150 (MQ=255)
                                                    ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                 >  2:10732928/1‑150 (MQ=255)
                                                       gcgggTGCTGGTGACCGCGATTGCTTTAAGCTTGTGGCGCAGGTTTTACGCCAACCCGTGCGCAACATTGCGCGCGATCGGCCCCCCCGCCCCTCTGGGGTGCGGGTTTGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                              <  1:7963171/147‑1 (MQ=255)
                                                            ggcgggTGCCAGCGATGGCGTCTCGTTTGTTGCGAAGGCGACAGGGCAAACGGACCGCAACATTGAGGGCGTGGGGCCCCCTCGCCCCTCTGGGGTTTGGGTGAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                         <  2:2047791/146‑1 (MQ=255)
                                                             gCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGACCCTATCAAcc                                                                                        >  2:7779349/1‑150 (MQ=255)
                                                              cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                       >  2:5329923/1‑150 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                      >  1:6687801/1‑150 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                      >  1:3368715/1‑150 (MQ=255)
                                                                 gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    >  2:2384175/1‑150 (MQ=255)
                                                                  tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAGGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                   >  1:2629609/1‑150 (MQ=255)
                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                               >  2:2594574/1‑150 (MQ=255)
                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                               >  2:2812058/1‑150 (MQ=255)
                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                               >  1:5185101/1‑150 (MQ=255)
                                                                        cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                             >  2:4001334/1‑150 (MQ=255)
                                                                              cggtAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCCTCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAAt                                                                       >  2:1624661/4‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCTGGAGATACGGCAAACCGTTCGCAACATTGAGGGCGTGGGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                      <  2:4539147/150‑1 (MQ=38)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                      >  1:7390605/1‑150 (MQ=255)
                                                                                ttATGATTGTGGCGCTGGAGATACGCCAAACCGACCGCAACTTTGCGCGCGATCGGCCCCCTCGCCCGTGTGGGTAGTGGTTTTGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                     <  1:1427098/150‑1 (MQ=255)
                                                                                 tAAGATTGTTGCGCCGGAGATAAGGCAAACCGATCGCAACTTTGCGCGCGATGGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                    <  2:7978494/150‑1 (MQ=255)
                                                                                 tAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                    >  1:8350125/1‑150 (MQ=255)
                                                                                   tttttgttgCGCTGTTGTTAAGGAAAACCGATCGAAATTTTGGGGGGGTTGGGTCCCCGCGCCCGTGGGGGGTGTGGTTTTGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                  <  1:4179465/147‑1 (MQ=255)
                                                                                       tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                              >  1:7704331/1‑150 (MQ=255)
                                                                                         ttgCGCAGGAGATAAGGAAAACCGATCGCAACATTCAGCGCGATCGGTCCCATCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCCATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                            >  1:374961/1‑150 (MQ=255)
                                                                                              cAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                       >  2:8450264/1‑150 (MQ=255)
                                                                                               cggAGATACGGCAACCCTAGCGCAACTTGGCGCGCGCGCGGGCCCCGCGCCCCTGGGGGGTGTGGTTTAGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                      <  2:4204584/149‑1 (MQ=255)
                                                                                                  agaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                   >  1:9472462/1‑150 (MQ=255)
                                                                                                    aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAATGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                 >  2:8587202/1‑150 (MQ=255)
                                                                                                    aTAAGGCAAACCGATCGCAACATAGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                 >  2:6860852/1‑150 (MQ=255)
                                                                                                         gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                            >  1:8929283/1‑150 (MQ=255)
                                                                                                           accccGATCGAAAATTGGCGCGCGAGCGGCCCCCGCGCCCCGCTGGGGTGTGGGTTTGGGTGAGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                          <  1:1344021/147‑1 (MQ=255)
                                                                                                            aaCCGATCGAAAATTTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGATTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                         <  1:9684746/150‑1 (MQ=255)
                                                                                                             aCCTATCGCAACATTGCGCGCCCACGGCCCCCTCGCCCCCCTGGGGAGAGGGTTAGGGTGTGGGCAAAACGCCGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                        <  1:6441355/150‑1 (MQ=255)
                                                                                                                gATCGCAAATTTGAGCGCTATGGGCCCCCTCCCCCCTCTGGGTTGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                     <  2:9467254/150‑1 (MQ=255)
                                                                                                                     cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                >  1:1817544/1‑150 (MQ=255)
                                                                                                                           ggAGCGCGATCGGCCCCCCCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                          <  1:2375570/149‑1 (MQ=255)
                                                                                                                             cgcgcGATCGGTCCCCTCGCCCCTGGGGGGTGTGGTTTGGGGTGAGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                        <  2:4907654/149‑1 (MQ=255)
                                                                                                                             cgcgcGATCGGGCCCCGCGCCCCGCGGGGGTGTGGGTTCGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGATTCTAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCTGTTCGTGTTcc                        <  2:9118584/149‑1 (MQ=255)
                                                                                                                             cgcgcGAGGGGGCCCCTCGCCCCTCGGGGGAGCGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCACATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                        <  2:1247096/149‑1 (MQ=255)
                                                                                                                               cggtaTCGGTCCCCTCGCCCCTCTGGGGTGTGTGTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                      <  2:5683925/146‑1 (MQ=255)
                                                                                                                                gcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCATTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCaga                     >  1:6947910/1‑148 (MQ=255)
                                                                                                                                  gATCGGCCCCCCCGCCCCTCTGGGTTGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                   <  1:208001/150‑1 (MQ=255)
                                                                                                                                   aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  >  2:10930463/1‑150 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCCGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                >  1:5611873/6‑150 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                >  1:6031323/6‑150 (MQ=255)
                                                                                                                                      ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa               >  2:3457404/1‑150 (MQ=255)
                                                                                                                                      ggTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCACCCAGTTCGTGTTCCATAAAGGtc               >  1:8675141/1‑149 (MQ=255)
                                                                                                                                      ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAGAAAGGTa               >  2:9609079/5‑150 (MQ=255)
                                                                                                                                       gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc              >  2:9165861/1‑150 (MQ=255)
                                                                                                                                       gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc              >  1:930208/1‑150 (MQ=255)
                                                                                                                                        tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa             >  2:7136189/1‑150 (MQ=255)
                                                                                                                                          cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa           >  2:6668804/1‑150 (MQ=255)
                                                                                                                                          cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa           >  2:5430215/1‑150 (MQ=255)
                                                                                                                                           ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa          >  2:1833733/1‑150 (MQ=255)
                                                                                                                                               gCCCCTGTGGGTAGTGTTTTTGGGTTTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  2:5232430/150‑1 (MQ=255)
                                                                                                                                                ccccTCTGGGGAGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg     <  2:3670714/150‑1 (MQ=255)
                                                                                                                                                 cccTGTGGGGTGTGGTTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg    <  2:4399852/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608546‑3608842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: