Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,696 T→G 5.3% intergenic (+47/+55) 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,6960TG 5.3% 612.2 / 5.9 208intergenic (+47/+55)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (90/99);  new base G (3/11);  total (95/113)
Fisher's exact test for biased strand distribution p-value = 9.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608547‑3608844
                                                                                                                                                     |                                                                                                                                                    
 acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                     >  1:431668/1‑150 (MQ=255)
  cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  2:1452592/150‑1 (MQ=255)
  cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:16787/150‑1 (MQ=255)
   gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCATATACGGGGGCGACGGTGCTGGTGACAGCTAATGCGTTAAGATTGTTGCGCAGTATATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  1:3659820/150‑1 (MQ=255)
   gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  2:3121476/150‑1 (MQ=255)
        tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              >  2:793276/1‑150 (MQ=255)
        tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:1013938/150‑1 (MQ=255)
        tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  1:2811690/150‑1 (MQ=255)
           ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                           >  2:1293731/1‑150 (MQ=255)
                   gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   >  1:1869167/1‑150 (MQ=255)
                    ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGCCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  1:54342/150‑1 (MQ=255)
                       tgtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACAGTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                               >  2:2800191/1‑150 (MQ=255)
                           gctggcCGTGCTGGCAGATCCGGGGGCGACGGCGCTGGTGACAGCGCCTGCGTTAAGCTTGTTGCGCCGGAGACAAGCCACACCGCTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCCCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           <  1:1027385/150‑1 (MQ=37)
                           gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           <  1:3498492/150‑1 (MQ=255)
                             tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                         >  2:3757388/1‑150 (MQ=255)
                                 aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                                                     >  1:2384260/1‑150 (MQ=255)
                                  gTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGAGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                                    >  2:182632/1‑150 (MQ=255)
                                      tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  2:482281/1‑150 (MQ=255)
                                       ggCCGATCCGGGGGCGCGGGTGGTGGTGACAGCGTTTGGTTTAATTTTTTTGCGCTGGATTTACGGCAACCCGATCGCAAAATTGAGCGGGATGGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                               <  2:2499455/150‑1 (MQ=16)
                                       ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                               >  2:2900669/1‑150 (MQ=255)
                                        gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                              >  1:878847/1‑150 (MQ=255)
                                        gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                              >  1:3177541/1‑150 (MQ=255)
                                        gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                              >  1:211160/1‑150 (MQ=255)
                                                   ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATGGGTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGCGCGCGATCGGTCCCCGCGCCCCACTGGGGTGAGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                   <  2:3373858/150‑1 (MQ=255)
                                                    gCGACGGTGCTGGTGACAGCGAATGCTTTAATATTGTTGCGCAGGAGATAAGGCAAACCTATCGCAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                  <  2:603845/150‑1 (MQ=255)
                                                     cGACGTGGCGGGTGACCGCGAATGGGTTAAGTTTGTTGCGCTGGAGATACGGCAACCCTATCGCAACATTTAGCGCGATGGGTCCCCTCGCCCCGCGGGGTAGCGGTTTTGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                 <  2:1647988/150‑1 (MQ=255)
                                                            gcgggTCACAGCTAATGCGTTAATATTGTTGCGCTGGAGACACGGCAAACCGATCGCAACTTTGCGGGCGATCGGTCCCCTCGCCCCTGTGGGGTGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                          <  2:2863877/147‑1 (MQ=255)
                                                            gCTGGTGACAGCGATTGCGTTATGATTGGTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGCGCGCGAGCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                          <  2:2743867/150‑1 (MQ=38)
                                                            gCTGGGGCCAGCTATTGCTTTAAGTTTGTTGCGCAGGAGACAAGGCAACCCGATCGCAACATTGCGCGCGAGGGGTCCCCTCGCCCCTCGGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATTGAAGCCAATCAGCCCCTATCAAcc                                                                                          <  2:3269202/150‑1 (MQ=255)
                                                              tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        <  2:642052/150‑1 (MQ=255)
                                                                 tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                     >  1:2337294/1‑150 (MQ=255)
                                                                  gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                    >  2:1946605/1‑150 (MQ=255)
                                                                   aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                   >  1:2538432/1‑150 (MQ=255)
                                                                   aCAGCGAATGCGATAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                   >  1:290119/1‑150 (MQ=255)
                                                                    ccgCGAATGCGTTAAGTTTGTTGCGCAGGATATAAGGCAAACCGATCGCAACATTGCGCGCCAGCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTCGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                  <  2:3908969/148‑1 (MQ=38)
                                                                    cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                  >  1:3998198/1‑150 (MQ=255)
                                                                     aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                 >  1:2481891/1‑150 (MQ=255)
                                                                       cgtttGCTTTAATATTTTTGCGCAGGATATATGGCAAACCGATCGCAACTTTGTGCGCGATCGGCCCCCTCGCCCCTGTGGGTAGTGGTTTCGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               <  2:3834703/146‑1 (MQ=255)
                                                                       cgtaTGCGTTAAGATGGTGGCGCAGGAGATCAGGCAAACCGATCGCAACATTGCGCGGGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTTCGaa                                                                               <  2:3753052/147‑1 (MQ=38)
                                                                        gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              >  1:1501574/1‑150 (MQ=255)
                                                                          aTGCGTTAAGATTGTGGCGAAGGACATAAGGCAAACCGACCGCAACATTGCGCGCGATCGGGCCCGTCGCCCCTCTGGGGAGTGGGTTTGGTTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCa                                                                            <  1:4018166/150‑1 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:2491969/1‑150 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:802669/1‑150 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:1331461/1‑150 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:3951057/1‑150 (MQ=255)
                                                                               ttatgtttgttgCGCTGGTGTTCCGGAAACCCGATCGCAAATTTGGGCGGGTGGGGGCCCCGCGCCCGTGTGGGTTGTGGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       <  2:393236/144‑1 (MQ=255)
                                                                               ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       >  1:1120273/1‑150 (MQ=255)
                                                                               ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       >  1:845512/1‑150 (MQ=255)
                                                                                  aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                    >  1:354930/1‑150 (MQ=255)
                                                                                    aTTGTTGCGCAGGAGATACGGCAACCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                  <  2:3108215/150‑1 (MQ=255)
                                                                                    aTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                  >  2:663409/1‑150 (MQ=255)
                                                                                     ttgttgCGCAGGAGATATGGCAACCCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 <  2:3998709/150‑1 (MQ=255)
                                                                                     ttgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCATCGCCCCTCTGGGTAGAGGTTTTGGTTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 <  2:2563457/150‑1 (MQ=255)
                                                                                     ttggtgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 <  2:4125043/150‑1 (MQ=255)
                                                                                     ggggtgCGCCGGAGATACGGCAAACCGACCGCAACATTGACCGCGATCGGGCCCCTCCCCCCTCGGGGGCGTGGGTCAGGGGGAGGGGAAAAGGCGGCATCGACGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 <  2:3055189/146‑1 (MQ=255)
                                                                                        ttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                              <  2:688874/150‑1 (MQ=255)
                                                                                         tgCGCTGGTGATACGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                             <  2:2550843/150‑1 (MQ=39)
                                                                                         tgCGCTGGAGTTTCGGCAACCCTATCGCAACTTTGTGCGCGTTCGGGCCCCTCGCCCGTGTGGGGTGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                             <  1:2261358/150‑1 (MQ=255)
                                                                                         tgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                             >  1:1557565/1‑150 (MQ=255)
                                                                                              aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                        >  1:597185/1‑150 (MQ=255)
                                                                                                   aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                   >  1:3465003/1‑150 (MQ=255)
                                                                                                     ccggCACACCTATCGAACCTTTGCGCGCGCTCGGGCCCCGCGCCCGTGGGGGGTGTGGGTGTGGGGGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 <  1:2899631/148‑1 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  1:1218370/1‑150 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGGCCATCCGTCTCtt                                              >  2:1142187/1‑150 (MQ=255)
                                                                                                         caccccGATCTCAACTTTGCGCGCTTTGGGCCCCCCCCCCCTTGTGGGTTGTGGGTTTGGTGGGGGGCACAAGGCGGTATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             <  1:687316/146‑1 (MQ=255)
                                                                                                         cAACCCGATCGAAACTTTGCGCGCGATGGGTCCCCGCGCCCGGGGGGGGTGGGGTTTCGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             <  1:1155632/150‑1 (MQ=255)
                                                                                                          aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCCTCCGTCTCTTTa                                            >  2:204846/1‑150 (MQ=255)
                                                                                                          aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  1:3117418/1‑150 (MQ=255)
                                                                                                            aCCGATCGAACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                          <  2:2096430/150‑1 (MQ=255)
                                                                                                              ctaTCGCAACTTTGCGCGCGTTGGGGCCCCCCGCCCCTGTGGGTTGTGGGTTGGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                        <  2:3529054/148‑1 (MQ=255)
                                                                                                              cGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                        >  1:2740652/1‑150 (MQ=255)
                                                                                                               tctCGCAACATTGCGCGCGATCGGTCCCCGCGCCCCTCTGGGGTGAGGTTTTGGTGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       <  1:4126787/148‑1 (MQ=255)
                                                                                                               gATCGCAACATTGCGCGCGAGCGGCCCCCTCGCCCCTCGGGGGCGCGGGTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       <  2:3607943/150‑1 (MQ=255)
                                                                                                                aTCGCAACATTCAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGATAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGCCAGTCCATCCGTCACTTTAGCAAcc                                      >  2:36743/1‑150 (MQ=255)
                                                                                                                    caccATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  2:2456748/4‑150 (MQ=255)
                                                                                                                    cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  2:2714457/1‑150 (MQ=255)
                                                                                                                    cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  1:1816601/1‑150 (MQ=255)
                                                                                                                      aCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGTAGTGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCACCAGTTCgt                               <  1:3092817/150‑1 (MQ=255)
                                                                                                                           gAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                           <  2:2318650/150‑1 (MQ=255)
                                                                                                                                  ttggttCCCCTCGCCCCTGTGGGGTGTGGTTTTGGGTGGGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  2:1869167/145‑1 (MQ=255)
                                                                                                                                  aTCGGTCTCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  2:2079531/150‑1 (MQ=255)
                                                                                                                                  aTCGGCCCCCTCGCCCCTCTGGGTAGTGGTTTTGGGTGTGGTGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  1:263015/150‑1 (MQ=255)
                                                                                                                                   tggggccccGCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                   <  2:431668/145‑1 (MQ=255)
                                                                                                                                   tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                   >  1:2206157/1‑150 (MQ=255)
                                                                                                                                    cGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                  <  2:2011803/150‑1 (MQ=255)
                                                                                                                                          ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa            <  1:3757367/150‑1 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc           >  1:1345883/1‑150 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc           >  1:1592990/1‑150 (MQ=255)
                                                                                                                                              gCCCCTCTGGGGAGAGGGTTCGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  2:1664904/150‑1 (MQ=255)
                                                                                                                                               ccccTCTGGGTTGTGGTTTTGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       <  1:1979140/150‑1 (MQ=255)
                                                                                                                                               ccccTCTGGGTTGAGGGTTTGGTTTTGGTGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       <  2:542013/150‑1 (MQ=255)
                                                                                                                                               ccccTCTGGGTAGTGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       <  2:3795062/150‑1 (MQ=255)
                                                                                                                                                 ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa     >  2:2384016/1‑150 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa    >  1:2918693/1‑150 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa    >  1:2985657/1‑150 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat   <  1:1946605/150‑1 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat   >  1:2580070/1‑150 (MQ=255)
                                                                                                                                                    ctGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  <  2:537456/150‑1 (MQ=255)
                                                                                                                                                    ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCACCCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  <  2:3755114/150‑1 (MQ=255)
                                                                                                                                                    ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  >  1:3294921/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608547‑3608844

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: