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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950CG14.8% 541.0 / 28.4 183intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (89/64);  new base G (2/27);  total (91/92)
Fisher's exact test for biased strand distribution p-value = 2.26e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAAT  >  NC_000913/3608546‑3608843
                                                                                                                                                     |                                                                                                                                                    
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCTCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                     >  2:2912284/1‑150 (MQ=255)
   cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                   <  2:2166796/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                  <  1:930582/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                  <  2:2764591/150‑1 (MQ=255)
     cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                 <  1:152531/150‑1 (MQ=255)
       gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                               <  2:2202540/150‑1 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                              <  2:2176709/150‑1 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                              <  1:1537105/150‑1 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  1:2588021/150‑1 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  1:2272350/150‑1 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCGGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  1:1701865/150‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            >  2:466374/1‑150 (MQ=255)
            ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgggg                                                                                                                                          >  2:1366376/1‑149 (MQ=255)
             gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgggga                                                                                                                                         <  2:988399/150‑3 (MQ=255)
                gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                      >  1:1990586/1‑150 (MQ=255)
                gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                      <  1:28989/150‑1 (MQ=255)
                 aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                     <  1:358368/150‑1 (MQ=255)
                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                   >  1:2031005/1‑150 (MQ=255)
                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                   >  1:1164215/1‑150 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACTGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                 <  2:231289/150‑1 (MQ=255)
                           ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                           >  1:2051129/1‑150 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                          >  1:750447/1‑150 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGaaaa                                                                                                                          >  2:2872972/1‑150 (MQ=255)
                               ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTGAAGATGGTGGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                       >  2:90890/1‑150 (MQ=255)
                                  aGTGCTGGCAGATACGGGGGCGCCGGTGCGGGTGACCGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAAggcggc                                                                                                                    <  2:3024287/150‑1 (MQ=255)
                                  aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                                                    >  2:445657/1‑150 (MQ=255)
                                    tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                  >  1:174907/1‑150 (MQ=255)
                                    tGCTGGCAGAAACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCGCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                  >  1:1108095/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTCGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                               >  2:792988/1‑150 (MQ=255)
                                        ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGATTGCTTTAAGATTGTTGCGCAGGAGATACGCCAACCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGTGTGGTTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                              <  2:2846126/150‑1 (MQ=255)
                                        ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGCGCGCGATCGGTCCCCTCGCCCCTGGGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                              <  2:1550725/150‑1 (MQ=255)
                                        ggCAGATACGGGGGCGACGGTGCTGGGGACTGCTATTGCGTTAAGATTGTTGCGCAGGATATACGGCAAACCGATCGAAACATTGAGGGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                              <  2:1807489/150‑1 (MQ=38)
                                         gCATATCCGGGGGCGCGGGTGCTTGTGACAGCGAAGGCTTTATTTTTGTGGCGCGGGAGATCAGGCAAACCGTTCGCAACTTTGAGGGCTATCGGTCCCCTCGCCCCTGTGGGGTGCGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAg                                                                                                             <  2:2451675/150‑1 (MQ=255)
                                         gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                             >  1:1672830/1‑150 (MQ=255)
                                         gCAGATACGGGGGCGACGGTGATGTACACAGCGAATGCGTTAAGATTGGTGCGCAGGAGATAACGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGCAGAGGGTTAGGGTGAGGGGAAAAGGCGTCATCGAAg                                                                                                             >  2:2684068/1‑150 (MQ=38)
                                              taggggggCGCGGGTGGTGGTGACCGCGATTGCTTTAAGATTGTGGCGCCGGATATAAGCCAACCCGATCGCAACTTTGCGCGCGATGGGTCCCCTCGCCCCTCTGGGTTTAGGGTTTGGGTGCGGGTAAAAGGCGGCATCGAAGCCAAt                                                                                                        <  2:2414685/147‑1 (MQ=255)
                                                   gggCGCCGGGGCTGGTGACCGCGATGGCGTTAATTTTGGGGCGCCGGAGATACGGCAACCCGATCGCAACATGCCGCGCTAAGGGCCCCCACGCCCCCCTGGGGAGGGGGTTCGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                   <  2:1188077/150‑1 (MQ=255)
                                                       gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                               >  1:2205168/1‑150 (MQ=255)
                                                        ccGGTGCGGGTTACCGCGAATGCTTTATTTTTTTTGCGCTGCATATACGCAAAACCTATCGCAACTTTGCGCGGGATCGGACCCCCCGCCCCTCTGGGGCGTGGGTTAGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTAt                                                                                              <  2:874714/149‑1 (MQ=255)
                                                          ggTGCTGTTGACACCGAATGCTTTAAGATTGTTGCGCAGGATATAAGGAAAACCGATCGCAACATGTAGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                            <  2:475995/150‑1 (MQ=38)
                                                          ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                            >  2:2997869/1‑150 (MQ=255)
                                                           gtggtggtgACCGCTAATGCTTTTAGTTGGTTGCGCTGGAGATACGGCACCCCGATCGCAACTTTGCGCGCGATGGGCCCCCTCGCCCAGCGGGGTAGTGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCaa                                                                                           <  2:2400335/150‑1 (MQ=255)
                                                               tGGTGACCGCTATTGCTTTAAGTTTTTTGCGCTGGATATACGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGTTGAGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                       <  2:112712/150‑1 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAGCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                       >  1:51028/1‑150 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                       >  2:3064252/1‑150 (MQ=255)
                                                                  tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                    >  1:2485151/1‑150 (MQ=255)
                                                                   gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                   >  1:574926/1‑150 (MQ=255)
                                                                    aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                  >  1:1340312/1‑150 (MQ=255)
                                                                     cAGGGATTGCGTTAAGTTTGTTGCGCTGGAGATAAGGCAACCCGCTCGCAACATTGTGCGCGTTCGGTCCCCGCGCCCCTCTGGGGTGAGGGTTTGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                 <  2:1265547/150‑1 (MQ=255)
                                                                       gctattGCTTTATGTTTGTTGCGCAGGTGACACGGCAAACCGATCGCAACATTGCGGGCGAGCGGTCCCCTCGCCCCTCGGGGGTGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                               <  2:278050/145‑1 (MQ=255)
                                                                       gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                               >  2:2134017/1‑150 (MQ=255)
                                                                       gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCGTTACGa                                                                               >  2:608785/1‑150 (MQ=255)
                                                                        cGATTGCTTTATGATTGTTGCGCTGGATATACGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                              <  2:971814/150‑1 (MQ=38)
                                                                               gTTATGTTTGTTGCGCAGGATATAAGGCAACCCGATGGCAACATTGAGCGCGATCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                       <  2:1977346/150‑1 (MQ=38)
                                                                                ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                      >  2:265635/1‑150 (MQ=255)
                                                                                     aTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                 >  2:940578/1‑150 (MQ=255)
                                                                                       tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                               >  2:601687/1‑150 (MQ=255)
                                                                                        gttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                              >  1:2696788/1‑150 (MQ=255)
                                                                                          tgCGCTGTATATAAGGAAAACCGTTCGAAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGTTTAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                            <  2:362034/150‑1 (MQ=255)
                                                                                               ttgagaTACGGAAAACCTCTCTCAACTTTGTGCGCGTTGGCTCCCCTCGCCCGTGGGGGTTTTGGGTGGGGGGGTGGGAAACAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                       <  1:1114755/148‑1 (MQ=255)
                                                                                                ggATATAAGGCAACCCTATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGCGGTTTCGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                                                      <  2:1595086/150‑1 (MQ=255)
                                                                                                  agaTAAGGCAACCCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                    <  1:3040376/150‑1 (MQ=255)
                                                                                                     tacgcaaaaCCGATCGCAACATTGAGCGCTATCGGCCCCCTCGCCCCTGTGGGGAGTGGGTTCGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                 <  1:774350/144‑1 (MQ=255)
                                                                                                     tAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                 >  1:3006131/1‑150 (MQ=255)
                                                                                                      aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                >  1:1791177/1‑150 (MQ=255)
                                                                                                        ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                              >  2:1107314/1‑150 (MQ=255)
                                                                                                        ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCAATCCGTctct                                              >  1:1638698/1‑150 (MQ=255)
                                                                                                          cAACCCGATCGAAACTTTGCGCGGGATGGGCCCCCGCGCCCCTCTGGGGTGTGGTTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                            <  2:116825/150‑1 (MQ=255)
                                                                                                          cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTCTCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                            >  1:1358537/1‑150 (MQ=255)
                                                                                                           aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                           <  2:2063712/150‑1 (MQ=255)
                                                                                                           aaaCCGACCGAAACTTTGCGCGCGCTGGGTCCCCTCGCCCCTGGGGGGTGTGGGTTAGGGGGGGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                           <  2:17467/150‑1 (MQ=255)
                                                                                                               cGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                       >  1:2339644/1‑150 (MQ=255)
                                                                                                                gATCGCAACATTGAGCGCGTTCGGTCCCCTCGCCCCTCTGGGGAGAGTTTTTGGGTGAGGTTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                      <  2:1687894/150‑1 (MQ=255)
                                                                                                                gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                      >  2:1497525/1‑150 (MQ=255)
                                                                                                                gATCGCAACATTGAGCGCGACCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                      <  2:1687912/150‑1 (MQ=255)
                                                                                                                   cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                   >  2:2020512/1‑150 (MQ=255)
                                                                                                                      aaCATTGCGCGCGATCGGCCCCCGCGCCCCTGTGGGGAGAGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                <  1:466374/150‑1 (MQ=255)
                                                                                                                          ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                            >  1:1627131/1‑150 (MQ=255)
                                                                                                                           cgcccgcgACCGGCCCCCGCGCCCCGCGGGGGAGCGGGTGAGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                           <  2:2967856/146‑1 (MQ=255)
                                                                                                                             aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                         >  2:2791873/1‑150 (MQ=255)
                                                                                                                             aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                         >  2:2705545/1‑150 (MQ=255)
                                                                                                                              gcgcGACCGGCCCCCCCGCCCGTGGGGGGCGCGGGTTCGGTTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                        <  1:265285/150‑1 (MQ=255)
                                                                                                                                gcgTTCGGCCCCCCCGCCCCTGTGGGGTTTGGTTTTGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                      <  1:2129290/150‑1 (MQ=255)
                                                                                                                                  gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                    >  1:2519569/1‑150 (MQ=255)
                                                                                                                                   accctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                   >  1:2822786/8‑150 (MQ=255)
                                                                                                                                    tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTTTTTAGCAACCAGTTCGTGTTCCATAAAgg                  >  2:1226609/1‑150 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                 >  1:2953522/6‑150 (MQ=255)
                                                                                                                                       gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc               >  1:45853/1‑150 (MQ=255)
                                                                                                                                           ccgcGCCCGTGTGGGTTGTGGGTTGGGGTGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa           <  2:267535/147‑1 (MQ=255)
                                                                                                                                           ccTCGCCCGTGGGGGGTGTGGGTTTGGGTGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa           <  2:590500/150‑1 (MQ=255)
                                                                                                                                            cgccccccACTGGGGTGCGGGTTCGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCACCCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc          <  1:1598529/146‑1 (MQ=255)
                                                                                                                                               gCCCCTGTGGGGTGTGGGTTCGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa       <  1:2392364/150‑1 (MQ=255)
                                                                                                                                               gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa       >  1:2961849/1‑150 (MQ=255)
                                                                                                                                                ccccTGTGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg      <  1:357958/150‑1 (MQ=255)
                                                                                                                                                   ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa   >  1:647721/1‑150 (MQ=255)
                                                                                                                                                   ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa   <  2:2493288/150‑1 (MQ=255)
                                                                                                                                                   ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa   >  1:1478545/1‑150 (MQ=255)
                                                                                                                                                   ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa   <  2:937294/150‑1 (MQ=255)
                                                                                                                                                   ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa   >  1:298401/1‑150 (MQ=255)
                                                                                                                                                    tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat  >  2:1383329/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAAT  >  NC_000913/3608546‑3608843

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: