Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,721 A→C 9.9% intergenic (+72/+30) 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,7210AC 9.9% 374.8 / 7.1 133intergenic (+72/+30)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (86/31);  new base C (2/14);  total (88/45)
Fisher's exact test for biased strand distribution p-value = 3.74e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

GGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGG  >  NC_000913/3608573‑3608869
                                                                                                                                                    |                                                                                                                                                    
ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGTTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                                     <  1:2051431/150‑1 (MQ=255)
ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                                     <  1:2217926/150‑1 (MQ=255)
  ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                                                   <  2:192266/150‑1 (MQ=255)
      cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcgg                                                                                                                                               >  2:479757/1‑150 (MQ=255)
      cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcgg                                                                                                                                               >  2:2077660/1‑150 (MQ=255)
      cAGTGCTGGCAGAGACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATGGTTGCGCAGGAGATAAGGCAAACCGATCGCAAGATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAATggcgg                                                                                                                                               <  2:1762297/150‑1 (MQ=39)
            tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                         >  1:1907786/1‑150 (MQ=255)
             ggCTGATACGGGGGCTTCGGTGCTGGTGACAGCGAATGCTTTAAGATTGTTGCGCTGGAGATAAGGCAAACCGATCGAAACATTGAGCGCGATCGGCCCCCTCGCCCCTGTGGGGTGATGTTTAGGGTGTGGGGAAAAGGCGGCATCGaa                                                                                                                                        <  2:2137431/150‑1 (MQ=255)
             ggCTGATACGGGGGCGACGGTGCTGGTGACAGCGATTGCTTTAAGATTGTTGCGCAGGATATAAGGCAACCCGATCGCAACATTGAGCGCTTTCGGGCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGTGGGGAAAAGGCGGCATCGaa                                                                                                                                        <  2:1204735/150‑1 (MQ=255)
              gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                                       >  2:1427063/1‑150 (MQ=255)
                aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAAGGCGGCATCGAAGcc                                                                                                                                     >  1:2043895/1‑150 (MQ=255)
                         ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                            >  2:1367547/1‑150 (MQ=255)
                           cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                          >  2:667645/1‑150 (MQ=255)
                           cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                          >  2:1097834/1‑150 (MQ=255)
                             ccGGTGGTGGTGACCGCTATTGCTTAAATTTTTTTGCGCAGGATATAAGGCAACCCTATCGAAACATTGAGCGCTATGGGTCCCCCCGCCCCTCTGGGTAGAGGGTTTGGGTTAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAt                                                                                                                        <  2:1708134/149‑1 (MQ=255)
                                 tGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                    >  2:1807815/1‑150 (MQ=255)
                                  gccggggACAGCGCACGCGTTAAGATTGTGGCGCAGGAGAGACGGCAACCCGATCGCAACATTGAGCGCGAACGGCCCCCTCGCCCCTCCGGGGAGAGGGTAAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                                                   <  2:594812/144‑1 (MQ=37)
                                  atctttGACAGCTAATGGTTTATTTTTTTTGCGCTGTAGATACGGCAACCCGACCGAAACTTGGCGCGCCATGGGCCCCCGCGCCCCTCTGGGGTGTGGTTTCGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                                                   <  2:36663/145‑1 (MQ=255)
                                    tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                 >  1:2177422/1‑150 (MQ=255)
                                    tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                 >  1:1653748/1‑150 (MQ=255)
                                    tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                 >  1:795442/1‑150 (MQ=255)
                                     ggTTACAGCGATTGCTTTATGTTTGTTGCGAAGTATATAAGGAAAACCTATCGAAAAATTGAGCGCTATCGGACCCCTCGCCCCTCTGGGTAGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                                                <  1:1773922/150‑1 (MQ=255)
                                       tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCgagat                                                                                                              >  1:36663/1‑146 (MQ=255)
                                         aCAGCTAATGCGTTATGTTTGTTGCGCTGTATATAAGGAAAACCTATCGCAACATTGTGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGGGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                                            <  2:142306/150‑1 (MQ=255)
                                         aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtttc                                                                                                            >  2:1087829/1‑149 (MQ=255)
                                                    gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  2:1057861/1‑150 (MQ=255)
                                                    gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCATCTGGGGAAAGGTTTAGGTTGAGGGCAAAAGCCGCAAACTAAGCCAAAAAGCCCAAAACAACCGCCTTTacaaaacaaaaa                                                                                                 >  2:626124/1‑148 (MQ=255)
                                                      tatttttgtggCGCTGTATATACGGCAAACAGATCCCAAATTTGCGGGCGTTGGGTCCCCGCGCCCGTGTGGGTTGTGGTTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                                               <  1:1687601/145‑1 (MQ=255)
                                                               tgCGCAGGATATATGGCAAACCGATCGAAACATTGCGCGCGATGGGTCCCCGCGCCCCTGTGGGGTGTGGGTTTGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                      <  1:338013/150‑1 (MQ=255)
                                                                    tgtttaTAAGTCAAACGGTTCCAAACTTTGCGCGCGCTGGGTCCCCCCGCCCGTTTGGGGTGTGGTTTTGGGTGGGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                                                 <  1:2077660/145‑1 (MQ=255)
                                                                       cgaTACGGCAAACGGAGCGCAACTTGGCGCGGGATGGCCCCCCTCGCCCCCGTGGGGTGTGGTTTTGGGTGTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                                              <  2:1624829/149‑1 (MQ=255)
                                                                         aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                                            >  1:1871014/1‑150 (MQ=255)
                                                                         aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                                            >  1:1171941/1‑150 (MQ=255)
                                                                            cgccacaCCTATCGAAACTTTGCGCGCGAGGGGTCCCCTCGCCCCTGTGGGGAGTGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                                                         <  2:1149743/144‑1 (MQ=255)
                                                                            cgccAACCCGATCGCACCATGGAGCGCGATCGGGCCCCGCGCCCCTCTGGGTTGAGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                                                         <  2:2140592/147‑1 (MQ=255)
                                                                            aGGCAACCCGATCGCAACTTTGAGCGCGATCGGCCCCCTCGCCCCTGTGGGGATTGGTTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                                                         <  1:804539/150‑1 (MQ=255)
                                                                                aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  1:718785/1‑150 (MQ=255)
                                                                                aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  1:718848/1‑150 (MQ=255)
                                                                                aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  1:864855/1‑150 (MQ=255)
                                                                                aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  1:984480/1‑150 (MQ=255)
                                                                                aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  1:1112002/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  1:927011/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  1:890103/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  1:2155260/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  1:162697/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  1:2290781/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  1:1435168/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTGTAGCa                                                                  >  1:2124232/1‑150 (MQ=255)
                                                                                   ccGATCGCAACATTGAGCGCGATAGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  2:1663960/1‑150 (MQ=255)
                                                                                         gCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                                            >  1:2041611/1‑150 (MQ=255)
                                                                                         gCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                                            >  1:1412910/1‑150 (MQ=255)
                                                                                          cAACTTTGCGCGCGTCCGGCCCCCTCGCCCCTGTGGGTTTTGGGTTCGGGGGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                           <  1:667645/150‑1 (MQ=255)
                                                                                          cAACATGGAGCGCGATCGGGCCCCTCGCCCCTCTGGGTTGCGGGTTTGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                           <  1:1097834/150‑1 (MQ=255)
                                                                                           aaCTTTGAGCGCGAGCGGTCCCCCCGCCCCTCTGGGGAGTGGGTTCGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                          <  1:942196/150‑1 (MQ=255)
                                                                                           aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                          >  1:2379501/1‑150 (MQ=255)
                                                                                           aaCATGGCGCGCTATCGGCCCCCCCCCCCCTCTGGGGAGAGGGTTTGGCTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcaatctt                                                          <  2:1949328/150‑7 (MQ=255)
                                                                                            aCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGTCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                                                         >  2:2053869/1‑150 (MQ=255)
                                                                                                tGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                                                     >  2:757331/1‑150 (MQ=255)
                                                                                                  aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                   >  2:1498208/1‑150 (MQ=255)
                                                                                                        ttCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                             <  2:795442/149‑1 (MQ=255)
                                                                                                        aTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                             <  2:795432/150‑1 (MQ=255)
                                                                                                        aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCTGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                             >  2:781887/1‑150 (MQ=255)
                                                                                                        aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                             <  2:2177422/150‑1 (MQ=255)
                                                                                                          cGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                           >  1:325324/1‑150 (MQ=255)
                                                                                                          cGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                           >  1:1224825/1‑150 (MQ=255)
                                                                                                            gTCCCCTCGCCCCTGTGGGGTGTGGGTTGGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                         <  1:1984755/150‑1 (MQ=255)
                                                                                                            gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                         >  1:1269263/1‑150 (MQ=255)
                                                                                                             tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                        >  2:534201/1‑150 (MQ=255)
                                                                                                             tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCGATCAACCGCCTTGACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTGTAGCACCCAGGTCGTGTGCCCGAAAGGTACa                                        >  2:1538546/1‑150 (MQ=38)
                                                                                                             tCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                        >  2:749164/1‑150 (MQ=255)
                                                                                                                 cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  1:2149168/1‑150 (MQ=255)
                                                                                                                 cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  2:621001/1‑150 (MQ=255)
                                                                                                                 cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  1:2261697/1‑150 (MQ=255)
                                                                                                                 cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  1:513977/1‑150 (MQ=255)
                                                                                                                 cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  1:1399655/1‑150 (MQ=255)
                                                                                                                 cTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  1:1399620/1‑150 (MQ=255)
                                                                                                                  tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                   >  1:313584/1‑150 (MQ=255)
                                                                                                                  tCGCCCCTCTGGGGAGAGGGATAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                   >  1:1530914/1‑150 (MQ=255)
                                                                                                                     ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                                >  1:2174372/1‑150 (MQ=255)
                                                                                                                        ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa                             >  1:794516/1‑150 (MQ=255)
                                                                                                                         tggggggTGTGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            <  1:482953/147‑1 (MQ=255)
                                                                                                                         tctGGGGCGTGGGTTAGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            <  2:2342985/150‑1 (MQ=255)
                                                                                                                         tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGact                            >  2:1971229/1‑148 (MQ=255)
                                                                                                                         gctGGGGTGTGGGTTTGGTTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            <  1:516375/149‑1 (MQ=255)
                                                                                                                             gggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg                        >  1:1631813/1‑150 (MQ=255)
                                                                                                                             gggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg                        >  2:1572179/1‑150 (MQ=255)
                                                                                                                              ggAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc                       <  2:2410075/150‑1 (MQ=255)
                                                                                                                               gtgtgggTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgcg                      <  2:872892/146‑1 (MQ=255)
                                                                                                                                    ggTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg                 >  1:1370883/1‑150 (MQ=255)
                                                                                                                                    ggTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg                 >  1:184828/1‑150 (MQ=255)
                                                                                                                                      ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg               >  2:310577/1‑150 (MQ=255)
                                                                                                                                      ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg               >  1:1384887/1‑150 (MQ=255)
                                                                                                                                      ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg               >  2:1759971/1‑150 (MQ=255)
                                                                                                                                      ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg               >  2:1245052/1‑150 (MQ=255)
                                                                                                                                         gggTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGcc            <  2:912539/150‑1 (MQ=255)
                                                                                                                                            tGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGa         >  2:1256508/1‑150 (MQ=255)
                                                                                                                                              cggggACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc       <  2:1458695/149‑1 (MQ=255)
                                                                                                                                                gggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgt     >  1:685993/1‑150 (MQ=255)
                                                                                                                                                   aaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAACCCCAGGAATATCGCGGGTAGTGGCCGGATCGTCgg  <  2:1281957/150‑1 (MQ=255)
                                                                                                                                                   aaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCgg  <  1:1348365/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGG  >  NC_000913/3608573‑3608869

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: