Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,721 A→C 9.8% 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.8% 668.7 / 13.4 237intergenic (+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 (135/73);  new base C (2/27);  total (137/100)
Fisher's exact test for biased strand distribution p-value = 1.33e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCG  >  NC_000913/3608572‑3608868
                                                                                                                                                      |                                                                                                                                                   
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                     >  1:2588492/1‑150 (MQ=255)
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                     >  1:2918342/1‑150 (MQ=255)
        aGTGATGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGAGTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGACAAggcggc                                                                                                                                             >  2:1305492/1‑150 (MQ=255)
             tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                        >  1:1175971/1‑150 (MQ=255)
               gctgATAGGGGGGCGACGGTGCTGGTGACAGCGATTGCGTTAATTTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAg                                                                                                                                      <  2:2745673/147‑1 (MQ=38)
                    tACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAAt                                                                                                                                 >  1:910233/1‑150 (MQ=255)
                    tACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAAt                                                                                                                                 >  2:2060507/1‑150 (MQ=255)
                     aCGGGGGCGCCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGTAGATAAGGCAAACCTCTCGCACCTTTGCGCGGGATGGGTCCCCGCGCCCCTCTGGGGAGAGGGTTCGGGGGCGGGGACAAGGCGGCATCGCAGCCAATc                                                                                                                                <  2:1995191/150‑1 (MQ=255)
                      cGGGGGCGCCGGTGCTGGTGACAGCGAATGCTTTATGTTTGGTGCGCAGGAGTACCCGCAAACCGATCCCAACTTTGAGCGCGATCGGGCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                                               <  2:2355355/150‑1 (MQ=25)
                      cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                                               >  2:1240956/1‑150 (MQ=255)
                       gggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAg                                                                                                                              >  1:2996177/1‑150 (MQ=255)
                        ggggggAGGGTGCTTGTGACAGCTATTGCTTTAAGTTTTTTGGGCAGTTTATAAGGAAAACCGTTCGCAACATTGAGCGCTATCGGTCCCCGCGCCCCTCGGGGTTGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                                             <  2:1909315/150‑1 (MQ=255)
                          ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                                          >  2:1000073/1‑150 (MQ=255)
                            cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                         >  2:2466269/1‑150 (MQ=255)
                            cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                         >  2:1311129/1‑150 (MQ=255)
                             gACGGTGGTGGTGACAGCTATTGCGTTATGTTTTTTGCCCTGTAGAAATGGAAAACCTATCGCAAATTTGAGCGCTATCGGTCCCCTCGCCCCTGTGGGTATTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                        <  2:1927944/150‑1 (MQ=255)
                                  tGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                   >  2:551114/1‑150 (MQ=255)
                                        tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                             >  1:2649780/1‑150 (MQ=255)
                                        tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                             >  1:575900/1‑150 (MQ=255)
                                           ccgCGAATGCGTTCAGATTGTTGCCCAGGATATACGGCAACCCGCTCGCACCTTGGAGCGCGTTGGGGCCCCTCGCCCCCCGGGGGAGAGGGTTGGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                          <  2:2815374/148‑1 (MQ=255)
                                           cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                          >  1:3649979/1‑150 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                         >  2:508412/1‑150 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                         >  1:117379/1‑150 (MQ=255)
                                              cgtattggTTATGTTTTTTGCGCTGCAGATAAGGCAACCGCTGCTCAACATTGAGCGCGTTCGGTCCCATCGCCACTCTGGGGTTTGGGTTAGGGTTAGGGTAAAAGGCGGCATCTAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                       <  2:2401132/143‑1 (MQ=255)
                                              cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                       >  1:1764575/1‑150 (MQ=255)
                                                     ggtAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                >  1:240331/3‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                >  1:2559324/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                >  1:1315088/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                >  1:834776/1‑150 (MQ=255)
                                                      ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                               >  1:1358248/1‑150 (MQ=255)
                                                        aaGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACg                                                                                             >  2:3807157/1‑150 (MQ=255)
                                                             tgttgCGCAGTAGATACGGCAACCCGATCGCACCCTTGAGCGCGACCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                                        <  2:4020226/150‑1 (MQ=38)
                                                             tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                                        >  1:22358/1‑150 (MQ=255)
                                                              gttgCGCCGTATATACGCCAACCCCATCCCAACTTTGAGCGCCTTCGGCCCCCTCGCCCCTGTGGGTAGCGGGTTTGGGTGTGGGGAAACGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                                       <  1:3256915/150‑1 (MQ=255)
                                                                tgCGCTGGATATAAGGCAAACCGATCGCAACATTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGTGAGGTTTAGGTTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                     <  2:2147329/150‑1 (MQ=255)
                                                                   tctgtttatatGGAAAACCGATCGATACTTTGCGCGGGAGCGGTCCCCGCGCCGGGGTGGGTCGTGTGTTCGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc                                                                                  <  1:624955/143‑1 (MQ=255)
                                                                     aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                                                >  2:64609/1‑150 (MQ=255)
                                                                        tgaTAAGGCAAACCGATCGCAACTTTGAGCGCGATGGGTCCCCTCGCCCCTGTGGGTAGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                                             <  1:3062478/149‑1 (MQ=255)
                                                                        agaTACGCCAAACCGAGCGCAACATTGCGCGCGATGGGTCCCCCCCCCCGTCGGGGGTGCGGGTTTGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                                             <  1:227247/150‑1 (MQ=255)
                                                                          aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                                           >  2:723409/1‑150 (MQ=255)
                                                                              ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                       >  2:3320707/1‑150 (MQ=255)
                                                                               gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGTAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                      <  2:672536/150‑1 (MQ=255)
                                                                               gCAAACCGATCGCAACATGGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                      <  2:3686778/150‑1 (MQ=255)
                                                                                cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                     <  1:2730988/150‑1 (MQ=255)
                                                                                 aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGTGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                    <  2:2638516/150‑1 (MQ=255)
                                                                                 aaaCCGAGCGCAACTTGGCGGGCGTGCGGTCCCCGCGCCCCTGTGGGGTGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                    <  1:2522898/150‑1 (MQ=255)
                                                                                 aaaCCGACCGCACCACGGCGCGCGAGCGGACCCCCCCCCCCTCTGGGTTGAGGGTTAGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                    <  1:1647769/150‑1 (MQ=255)
                                                                                  acccgaccgaAACTTGGCGCGCTCTCGGCCCCCTCGCCCGTGTGGGGTGTGGGTTTGGGGCCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                                                   <  1:1435968/148‑1 (MQ=255)
                                                                                      gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                               >  2:2766683/1‑150 (MQ=255)
                                                                                       gcccGCAACTTGGCGCGCGATCGGCCCCCGCGCCCCTGTGGGGTGTGGGTTGGGGGGTGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                               <  1:716982/147‑1 (MQ=255)
                                                                                        tCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                                             <  2:2044872/150‑1 (MQ=255)
                                                                                        tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                                             >  2:1719395/1‑150 (MQ=255)
                                                                                           gCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcccgg                                                           >  2:31113/1‑147 (MQ=255)
                                                                                            cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                          >  2:608699/1‑150 (MQ=255)
                                                                                            cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                          >  2:3511956/1‑150 (MQ=255)
                                                                                             aaatttGCGCGCTATGGGTCCCCTCGCCCCGCTGGGGAGTGGGTTCGGGTGTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCATCTCTTTAGCAACCAGTTc                                                         <  2:1772186/146‑1 (MQ=255)
                                                                                             aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                         >  1:2133487/1‑150 (MQ=255)
                                                                                             aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                         >  1:2970421/1‑150 (MQ=255)
                                                                                             aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                         >  2:609909/1‑150 (MQ=255)
                                                                                             aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACAAGTTc                                                         >  1:1514136/1‑150 (MQ=255)
                                                                                               cctggcgcgcGATGGGGCCCCACGCCCCTGTGGGGTGTGGGTTGGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                                                       <  2:3155911/144‑1 (MQ=255)
                                                                                               cATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                                                       >  2:4090900/1‑150 (MQ=255)
                                                                                               cATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                                                       >  1:3442948/1‑150 (MQ=255)
                                                                                                 ttGCGCGCGATCGGTCCCCTCGCCCGTGGGGGGTGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                                     <  1:2181948/150‑1 (MQ=255)
                                                                                                 ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTATACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTATAGCAACCAGTTCgagt                                                     >  1:3520257/1‑150 (MQ=255)
                                                                                                  tGACCGCGATCGGACCCCTCCCCCCTCTGGGGATTGGTTTAGGGTGAGGCCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                                                    <  1:528174/150‑1 (MQ=255)
                                                                                                   gcgcgcGATGGGCCCCCCCCCCCCTCGGGGGAGAGGTTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                                                   <  2:1092062/148‑1 (MQ=255)
                                                                                                   gAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                  <  1:1000073/150‑1 (MQ=255)
                                                                                                    cgcgcgcTGGGGCCCCGCCCCCCTCTGGGGAGGGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                  <  2:1901577/149‑1 (MQ=255)
                                                                                                    aGCGCGATCGGTCCCCTCGCCCCTATGGGGAGAGGGTTAGGGTGAGGGGAAAGGGCGGCATCGAAGCAAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGACTCTTTTGCACCCAGTTCGTGTTcc                                                  >  2:442848/1‑150 (MQ=38)
                                                                                                      ggcgttgggTCCCCTCGCCCCTATGGGGAGTGGGTTAGGGTGAGGTGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATACCGATAAGGCAGTCCATCCGACTCTTTAGCAACCAGTTCGTGTTCCAt                                                <  2:514063/143‑1 (MQ=255)
                                                                                                      cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                                                >  1:146003/1‑150 (MQ=255)
                                                                                                       gcgATCGGTCCCCTCGCCCCTCGGGGGTGCGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                                               <  2:575900/150‑1 (MQ=255)
                                                                                                         gATCGGGCCCCTCGCCCCTGTGGGGTGTGGTTTCGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                             <  2:173960/150‑1 (MQ=255)
                                                                                                             ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGTGAAAATGATGCATCGAAGCCAATCAGCACCTATCAACCGCCTTAACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACAAGTTCGTGTTCCATAAAGGTa                                         >  2:4083000/1‑150 (MQ=38)
                                                                                                               tCCCCTCGCCCCTCTGGGGAGAGGTTAAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCATTAAGAATCAAATAACGATAAGGCAGTCCATCAGTATCTTTAGCAACCAGTTCGAGTTACATAAAGGTACa                                       >  1:964283/1‑150 (MQ=38)
                                                                                                               tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                       >  2:2900006/1‑150 (MQ=255)
                                                                                                                acccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                                      >  1:3887793/2‑150 (MQ=255)
                                                                                                                  ccTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa                                    <  1:2934211/150‑1 (MQ=255)
                                                                                                                    tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                  >  2:2427488/1‑150 (MQ=255)
                                                                                                                      gCCCCTCTGGGTTGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                <  1:2594926/150‑1 (MQ=255)
                                                                                                                      cccccTCTGGGGAGAGGGTTCGGTGGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                <  1:3911228/149‑1 (MQ=255)
                                                                                                                           tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                           >  2:17235/1‑150 (MQ=255)
                                                                                                                           tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAACAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAACGTACAAAACCCAGGaaa                           >  2:2659389/1‑149 (MQ=255)
                                                                                                                             tGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAatat                         >  1:1043977/1‑150 (MQ=255)
                                                                                                                              ggggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATc                        >  1:4119154/1‑150 (MQ=255)
                                                                                                                                 gtgtgtgTTTGGGGGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgcg                     <  2:2257160/144‑1 (MQ=255)
                                                                                                                                   gagGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCggg                   >  1:1622301/1‑150 (MQ=255)
                                                                                                                                     gggTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTa                 <  1:3558296/150‑1 (MQ=255)
                                                                                                                                     gggTTCGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTa                 <  1:1816442/150‑1 (MQ=255)
                                                                                                                                           gggTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGcc           >  2:3064016/1‑150 (MQ=255)
                                                                                                                                            ggTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCg          <  1:1736338/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGGAAAAGGCGGCATCGAAGACAATCAGCCCCTATCAACCGCCTTTACGAATCATATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCg          <  2:3375233/150‑1 (MQ=255)
                                                                                                                                                aGGGGACCAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc      <  2:3754504/150‑1 (MQ=255)
                                                                                                                                                   ggACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc   <  1:2067622/150‑1 (MQ=255)
                                                                                                                                                   ggACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc   <  2:2996177/150‑1 (MQ=255)
                                                                                                                                                   ggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc   >  2:747305/1‑150 (MQ=255)
                                                                                                                                                   ggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGAAGTCCATCCGTCTCTTTAGCAACCAGTTAGTGTTCCATAAAGGTACAAAAACCATGAATATCGCGGGTAGTGGCCGGAtcgtc   >  2:3554042/1‑150 (MQ=255)
                                                                                                                                                    gAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtcg  >  2:3613344/1‑150 (MQ=255)
                                                                                                                                                    gAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtcg  >  1:1018557/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCG  >  NC_000913/3608572‑3608868

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