Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,671 A→C 6.3% intergenic (+22/+80) 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,6710AC 6.3% 1847.4 / 3.7 601intergenic (+22/+80)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (310/220);  new base C (2/56);  total (314/287)
Fisher's exact test for biased strand distribution p-value = 1.57e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 608 total aligned reads displayed.

GGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCC  >  NC_000913/3608522‑3608820
                                                                                                                                                     |                                                                                                                                                     
 gCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAg                                                                                                                                                      >  2:1576875/1‑150 (MQ=255)
 gCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAg                                                                                                                                                      >  2:6278129/1‑150 (MQ=255)
  cTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAgc                                                                                                                                                     >  1:2045483/1‑150 (MQ=255)
   tGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAgcg                                                                                                                                                    >  1:11103350/1‑150 (MQ=255)
    gAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAgcgc                                                                                                                                                   <  2:6899244/150‑1 (MQ=255)
       aGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGAt                                                                                                                                                <  2:9100122/150‑1 (MQ=255)
       aGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGCTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGAt                                                                                                                                                >  1:8412290/1‑150 (MQ=255)
         ggATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCg                                                                                                                                              <  1:5184200/150‑1 (MQ=255)
         ggATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCg                                                                                                                                              <  1:7771747/150‑1 (MQ=255)
           aTCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGt                                                                                                                                            >  2:8756009/1‑150 (MQ=255)
             cTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGttt                                                                                                                                          <  2:5863123/150‑3 (MQ=255)
               tCCTCGTCACCTCGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGGGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTcccc                                                                                                                                        <  2:5412355/150‑1 (MQ=255)
               tCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTcccc                                                                                                                                        >  1:5843754/1‑150 (MQ=255)
                  tCGTCACCCCGCTGTTTGGGATGACCGGGTTGTGGCTGGCAGGGCTGGCAGATACGGGGGCGAGGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCCTCg                                                                                                                                     <  1:10632832/150‑1 (MQ=255)
                  tCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCg                                                                                                                                     <  2:198807/150‑1 (MQ=255)
                  tCGTCACCACGCTGTTAGGGAGGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGTCAGCGAATGCGGTCAGATGGTTGCGCAGGAGATACGGCAATCCGATCGCACCATTGAGCGCGATCGGTCCCCTCg                                                                                                                                     >  1:4963935/1‑150 (MQ=255)
                   cGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGc                                                                                                                                    <  1:7432943/150‑1 (MQ=255)
                   cGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGc                                                                                                                                    >  1:1582940/1‑150 (MQ=255)
                    gTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGcc                                                                                                                                   >  2:7961150/1‑150 (MQ=255)
                     tCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGccc                                                                                                                                  <  2:1726653/150‑1 (MQ=255)
                        ccacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                               <  2:7073235/150‑1 (MQ=255)
                         cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                              >  1:4106155/1‑150 (MQ=255)
                         cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                              >  1:3490499/1‑150 (MQ=255)
                          acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                             <  2:3490499/150‑1 (MQ=255)
                            gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                           <  1:10836158/150‑1 (MQ=255)
                            gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                           <  2:6390538/150‑1 (MQ=255)
                            gCTGTTAGGGATCACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                           >  2:8491283/1‑150 (MQ=255)
                             cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                          <  1:2979529/150‑1 (MQ=255)
                                ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                       <  2:8829338/150‑1 (MQ=255)
                                ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                       <  1:1168024/150‑1 (MQ=255)
                                 tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                      <  2:5615623/150‑1 (MQ=255)
                                  aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                     >  1:9317652/1‑150 (MQ=255)
                                  aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                     <  2:375090/150‑1 (MQ=255)
                                     gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                  >  2:8358061/1‑150 (MQ=255)
                                        gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                               <  1:4441479/150‑1 (MQ=255)
                                        gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                               <  2:7807366/150‑1 (MQ=255)
                                        gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                               >  1:468362/1‑150 (MQ=255)
                                          ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATTCGTTAAGATTGTTGCGCAGGAGATAAGTCAAACCGATCGCAACATTGAGCGCGATCGGTCCCATCGCCCCTCAGTTGAGAGGGTTAGGGt                                                                                                             >  1:6377167/1‑150 (MQ=255)
                                          ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                             >  1:6902636/1‑150 (MQ=255)
                                           cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                            <  1:124463/150‑1 (MQ=255)
                                                    gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCACACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                   <  2:962848/150‑1 (MQ=255)
                                                     ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                  >  1:817020/1‑150 (MQ=255)
                                                     ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                  >  1:4239961/1‑150 (MQ=255)
                                                       ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                >  1:549578/1‑150 (MQ=255)
                                                        gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                               >  2:4238045/1‑150 (MQ=255)
                                                              cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCg                                                                                         >  1:2777041/1‑150 (MQ=255)
                                                               tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                        >  1:8542360/1‑150 (MQ=255)
                                                                ggCATATCCGGGGGCGACGGTGCTGGTGCACGCGAATGCTTTATTATTGTTGCGCAGTAGATATGGCAAACGTATCGCAACTTTGTGCGCTATCGGTCCCCTCGCCCCTCTGGGGTTAGGTTTTGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                       <  2:4378736/150‑1 (MQ=255)
                                                                 gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                      >  1:8098576/1‑150 (MQ=255)
                                                                  ctgATACGGGGGCGACGGTGCTGGTGACAGCGAATGCTTTATGATTGTTGGGATGTTGATATGGAAAACCTATCGCAACATTTAGGGCGATGGGTCCCCTCGCCCCTCTGGGGTGAGGTTTAGGGTGAGGGTAAAAGGCGGCATCGAAGc                                                                                     <  2:1211425/148‑1 (MQ=255)
                                                                  cAGATAGGGGGGCGAGGGTGATGTTGACAGGGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGAAAACGGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGc                                                                                     <  2:2606571/150‑1 (MQ=38)
                                                                   aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                    >  1:10363500/1‑150 (MQ=255)
                                                                   aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGAGTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                    >  2:6070415/1‑150 (MQ=255)
                                                                         gggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAg                                                                              >  2:10927710/1‑150 (MQ=255)
                                                                          ggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                             >  2:8138360/1‑150 (MQ=255)
                                                                            ggcgcgggTGCGGGTGACAGCGATGGCGTTAAGATTGGTGCGCAGGAGATAAGGCAAACCGCTCGCAACTTTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                           <  2:10138523/144‑1 (MQ=37)
                                                                              cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                         >  2:3081039/1‑150 (MQ=255)
                                                                              cGACGGTGCTGGTGACAGCGAATGCGTTAAGAGTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGGTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                         >  1:10932783/1‑150 (MQ=255)
                                                                               gcgggTGCTGGTGACAGGGAATGCGTTATTTTTGTTGCGCAGGATATAAGGCAACCCTTTCGAAACTTTGCGCGCGCTCGGTCCCCGCGCCCCTGGGGGGTGTGGTTTTGGGGGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                        <  2:2930157/147‑1 (MQ=255)
                                                                               gagggTGCTGGTGACAGCGACTGCGTTATGTTTGTTGCGCTGGAGATACGGCAAACCGATCGCATCCTTGAGCGCGAGGGGTCCCCTCGCCCGTGGGGGGAGTGGTTTCGGGTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                        <  2:6574284/147‑1 (MQ=255)
                                                                               gACGGTGCTGGTGACAGCGAATGCGTTATGTTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                        <  2:7387828/150‑1 (MQ=255)
                                                                               gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGAGAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCACGCAGCCCCTa                                                                        >  2:2458849/1‑150 (MQ=39)
                                                                                 gggTGCGGGTGAAAGGTAATGTTTTATTTTTGTTGCGCTGGATTTACGGAAAACCGATCGAAACTTTGAGCGGTCTGGGCCCCCGCGCCCCTCTGGGTAGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                      <  2:1259067/149‑1 (MQ=255)
                                                                                 cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                      <  2:10232682/150‑1 (MQ=255)
                                                                                 cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                      >  1:9999646/1‑150 (MQ=255)
                                                                                       tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                >  2:5107293/1‑150 (MQ=255)
                                                                                        ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                               >  1:2641265/1‑150 (MQ=255)
                                                                                         gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGCTAAGGCAAACCGATCGCAACATTGAGCGCGAGCGGTCCCCTCGCCCCTCTGGGGCGAGGGGTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTCTCAACCGCCt                                                              >  2:7319798/1‑150 (MQ=39)
                                                                                          tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                             >  2:9292930/1‑150 (MQ=255)
                                                                                             cAGCGTTTGCTTTAATTTTGTTGGGCTGGTTATAAGCAAACCCGATCGAAACTTTGCGCGCGCTGGGCCCCCCCGCCCCTGTGGGTTGTGGGTTTGGGTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                          <  1:7912049/150‑1 (MQ=255)
                                                                                              aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                         >  2:11062620/1‑150 (MQ=255)
                                                                                              aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                         >  1:3563443/1‑150 (MQ=255)
                                                                                                ctattgggTTAAGTTGTTTGCGCAGTATATAAGGCAACCCCAACGCAAATTGGACCGCTCGCGGTCCCCCCGCCCCCCTGGGGTGCTGTTTCGTGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                       <  1:5726444/143‑1 (MQ=255)
                                                                                                cgtaTGGGTTAAGTTTGTTGCGCAGGATATCCGGCACCCCTATCGCAACTGTGAGCGCTAGCGGGCCCCTCGCCCCTCTGGGTCGGGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                       <  2:979255/147‑1 (MQ=255)
                                                                                                cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                       >  1:7310921/1‑150 (MQ=255)
                                                                                                 gAATGCTTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGTGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                      <  2:1796447/150‑1 (MQ=255)
                                                                                                    tGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaa                                                   >  2:10588893/1‑150 (MQ=255)
                                                                                                         tAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                              >  2:8676451/1‑150 (MQ=255)
                                                                                                           tttttgttgCGCAGGTGATACGGCAACCCGATCGCAACTTTGGGCGGGATCGGCCCCCTCGCCCGTGGGGGTTGTGGGTTTGGGTGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                            <  2:4840895/147‑1 (MQ=255)
                                                                                                            gATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACTAt                                           >  1:10000428/1‑150 (MQ=255)
                                                                                                                        ggAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                               >  1:7131615/1‑150 (MQ=255)
                                                                                                                        ggAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                               >  2:3232609/1‑150 (MQ=255)
                                                                                                                           gaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                            >  1:980926/1‑150 (MQ=255)
                                                                                                                            ttAAGGCACACCGAGCGCAACATTGCGCGCGATCGGTCCCCGCGCCCCTGTGGGGTGTGGGTTCGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                           <  2:10404423/149‑1 (MQ=255)
                                                                                                                             aaaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                          >  1:8969885/2‑150 (MQ=255)
                                                                                                                              aaGCCAACCCGCGCGCAACATTGCGCGCGCGCGGCCCCCCCCCCCCTGGGGGGAGAGGTTTTGGGTGCGGGGCAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtca                         <  1:3885741/150‑2 (MQ=255)
                                                                                                                               aGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                        >  1:3273944/1‑150 (MQ=255)
                                                                                                                                ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                       >  2:3941723/1‑150 (MQ=255)
                                                                                                                                ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                       >  1:1355468/1‑150 (MQ=255)
                                                                                                                                 gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAATGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                      >  2:9350049/1‑150 (MQ=255)
                                                                                                                                 gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                      >  2:7242222/1‑150 (MQ=255)
                                                                                                                                 gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                      >  2:3364158/1‑150 (MQ=255)
                                                                                                                                 gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                      >  1:87581/1‑150 (MQ=255)
                                                                                                                                   aaaCCCACCGAAACTTGGCGCGCGTTGGGTCCCCCCGCCCCTGTGGGTTGTGTGTTCGGGTGTGGTGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                    <  1:8491327/150‑1 (MQ=255)
                                                                                                                                      ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                 >  2:2316212/1‑150 (MQ=255)
                                                                                                                                         aTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc              >  1:2511311/1‑150 (MQ=255)
                                                                                                                                              aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGAGGAGAGGGTTAGGGTTAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc         <  1:10178713/150‑1 (MQ=255)
                                                                                                                                                cATTGCGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt       <  1:6796429/150‑1 (MQ=39)
                                                                                                                                                 tttGCGCGCTATCGGTCCCCTCGCCCCTCGGGGGTGTGGTTTTGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg      <  2:405182/149‑1 (MQ=255)
                                                                                                                                                  ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt     >  1:6146933/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCC  >  NC_000913/3608522‑3608820

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: