Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 539,942 C→A 6.7% P774T (CCG→ACG)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947539,9420CA 6.7% 339.5 / 9.8 105P774T (CCG→ACG) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (50/48);  new base A (5/2);  total (55/50)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 105 total aligned reads displayed.

TGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCC  >  NC_002947/539803‑540089
                                                                                                                                           |                                                                                                                                                    
tgtgATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGGACGCTGGttg                                                                                                                                             <  7:205796/149‑3 (MQ=255)
  tgATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGGACGCTGGttgat                                                                                                                                           <  4:101818/149‑5 (MQ=255)
   gATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGc                                                                                                                                           >  7:303991/1‑148 (MQ=255)
     tCGCCTCGGTCGACGCCAGCCTTTTGGTTGTTCGGGTTGCCGACTACTATGTTTAAACCGGCTATGGGGGTGTGTATATCTACAACCTGACAATGTACACCCGTTCTAACCATAACACCTGCATCAACCAGCGTCCGCTGGt                                                                                                                                               <  8:178477/142‑1 (MQ=255)
     tCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGATCGCCTGCATCAACCAGCGTCCGCTGGTGAGCAAg                                                                                                                                        >  6:17058/1‑149 (MQ=255)
     tCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTgggatcg                                                                                                                                        >  7:178477/1‑143 (MQ=255)
     tCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCaa                                                                                                                                         <  4:195178/148‑1 (MQ=255)
     tCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCACCCTCCGCTGGAGAGCAAg                                                                                                                                        >  8:7024/1‑149 (MQ=255)
       taCTCGGTCGACGCAAGCCGTTTCGTTGTTCGCTTTGCCGACGACGAAGTGGAAACCGGCTAAGGGGGTGTGTAATTCTACAACGTGACAATGTACACCGGTTCGAACCATAACACGTGCATCAACCAGCGCCCGCTGGTTAGAAAgg                                                                                                                                       <  8:303991/147‑1 (MQ=255)
       gACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACCCCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAggg                                                                                                                                      >  4:231701/1‑149 (MQ=255)
       gACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAgg                                                                                                                                       <  7:245564/148‑1 (MQ=255)
       gACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGCGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAgg                                                                                                                                       <  1:376715/148‑1 (MQ=255)
        aCTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGt                                                                                                                                     >  8:289462/1‑149 (MQ=255)
        aCTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGt                                                                                                                                     >  6:116911/1‑149 (MQ=255)
         cTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGt                                                                                                                                     <  2:168821/148‑1 (MQ=255)
         cTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGt                                                                                                                                     >  8:167498/1‑148 (MQ=255)
            ggTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGa                                                                                                                                   <  6:208095/147‑1 (MQ=255)
            ggTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGAc                                                                                                                                  <  4:316866/148‑1 (MQ=255)
            ggTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGAc                                                                                                                                  >  7:238522/1‑148 (MQ=255)
            ggTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGAc                                                                                                                                  >  1:179553/1‑148 (MQ=255)
            ggTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGAc                                                                                                                                  >  1:185617/1‑148 (MQ=255)
              tCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAg                                                                                                                               <  1:197870/149‑1 (MQ=255)
                   gCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCa                                                                                                                          <  4:47907/149‑1 (MQ=255)
                     cAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGc                                                                                                                        <  3:156623/149‑1 (MQ=255)
                      aGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCg                                                                                                                       >  5:92399/1‑149 (MQ=255)
                        ccGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtg                                                                                                                     <  8:238522/149‑1 (MQ=255)
                          gTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  5:380368/149‑1 (MQ=255)
                          gTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  2:185617/149‑1 (MQ=255)
                          gTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  2:179553/149‑1 (MQ=255)
                             tCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGa                                                                                                                 <  7:289462/148‑1 (MQ=255)
                               gttgttCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGAca                                                                                                              >  8:184845/1‑148 (MQ=255)
                                 tgttCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCATGGGTGACAAGGTTCAGCGTGGTGAcatca                                                                                                           >  5:178520/1‑149 (MQ=255)
                                  gttCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCaa                                                                                                                                         <  7:269499/119‑1 (MQ=255)
                                  gttCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCaa                                                                                                                                         >  8:269499/1‑119 (MQ=255)
                                  gttCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGAcatcat                                                                                                          >  8:299984/1‑149 (MQ=255)
                                      gcgTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGcc                                                                                                      >  5:190590/1‑149 (MQ=255)
                                        gTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCgt                                                                                                    <  7:299984/149‑2 (MQ=255)
                                         ttGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGAc                                                                                                   >  8:83362/1‑149 (MQ=255)
                                            ccGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGc                                                                                                <  4:34806/149‑1 (MQ=255)
                                             cgacgacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGc                                                                                                <  7:184845/148‑1 (MQ=255)
                                                 gacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTc                                                                                           >  5:345977/1‑149 (MQ=255)
                                                  acgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTcc                                                                                          <  2:167864/149‑1 (MQ=255)
                                                    gaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   >  8:5884/1‑123 (MQ=255)
                                                    gaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  7:5884/123‑1 (MQ=255)
                                                    gaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCAc                                                                                        >  7:374649/1‑149 (MQ=255)
                                                    gaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCAc                                                                                        >  8:31307/1‑149 (MQ=255)
                                                           aaaCCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATg                                                                                 <  5:116911/149‑1 (MQ=255)
                                                            aaCCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATgg                                                                                <  6:311209/149‑1 (MQ=255)
                                                                 gCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAg                                                                           >  6:328624/1‑149 (MQ=255)
                                                                   gAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCt                                                                         >  2:91956/1‑149 (MQ=255)
                                                                         ggTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGATCGCCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGctct                                                                   <  5:17058/149‑1 (MQ=255)
                                                                           tgtgGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTgg                                                                 <  7:167498/149‑1 (MQ=255)
                                                                              ggATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGt                                                               >  3:1283/1‑148 (MQ=255)
                                                                               gATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGt                                                               >  2:32272/1‑147 (MQ=255)
                                                                               gATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCa                                                             >  1:273729/1‑149 (MQ=255)
                                                                                atatCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGt                                                               >  6:296623/1‑146 (MQ=255)
                                                                                atatCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGt                                                               >  1:310525/1‑146 (MQ=255)
                                                                                    cTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACa                                                        >  6:1747/1‑149 (MQ=255)
                                                                                     tACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACAt                                                       >  3:76681/1‑149 (MQ=255)
                                                                                     tACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACAt                                                       >  5:61075/1‑149 (MQ=255)
                                                                                      aCAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATg                                                      >  7:111387/1‑149 (MQ=255)
                                                                                      aCAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATg                                                      >  6:159858/1‑149 (MQ=255)
                                                                                      aCAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATg                                                      >  3:7795/1‑149 (MQ=255)
                                                                                       cAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgc                                                     <  2:310525/149‑1 (MQ=255)
                                                                                       cAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgc                                                     <  8:374649/149‑1 (MQ=255)
                                                                                        aaCCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgcg                                                    <  6:92399/149‑1 (MQ=255)
                                                                                        aaCCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgcg                                                    >  2:262202/1‑149 (MQ=255)
                                                                                            tGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATc                                                >  8:271677/1‑149 (MQ=255)
                                                                                              aCCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATc                                                <  5:1747/147‑1 (MQ=255)
                                                                                                 cgatgtacCCGTTGGAACCAGATCACCTGCATCAACCAGCGTCGGCTGGTGAGCAAGTGTGACATGGTTCAGCGTGGTGACATCTTGGCCGACGGCCCTTCCAGCGACATTGGTGAGCTGGCTCTGTGTCAGAACATGCGCATCGCGtt                                           <  8:25735/143‑1 (MQ=255)
                                                                                                   gaacacCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCa                                         >  7:112101/3‑149 (MQ=255)
                                                                                                      cacCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATg                                       <  6:345977/148‑1 (MQ=255)
                                                                                                        cccGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGtt                                           >  7:25735/1‑142 (MQ=255)
                                                                                                         ccGTTCGAACCAGAACACCTGCATCAACCAGCTTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCGCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgc                                                     <  5:418113/131‑1 (MQ=255)
                                                                                                         ccGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTTACAAGGTTCAGCGTGGTGACATCAAGGCCGACGGCCCGTCCACCGACATGGGTCAGCTGGCTATGCGTCAGAACATgc                                                     >  6:418113/1‑131 (MQ=255)
                                                                                                         ccGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGt                                   >  4:107546/1‑149 (MQ=255)
                                                                                                              cGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGaa                               <  3:107546/148‑1 (MQ=255)
                                                                                                                 aCCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGc                           <  1:91956/149‑1 (MQ=255)
                                                                                                                 aCCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGc                           <  7:271677/149‑1 (MQ=255)
                                                                                                                  ccAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGc                           <  7:92895/148‑1 (MQ=255)
                                                                                                                       acacCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAAc                     <  7:31307/149‑1 (MQ=255)
                                                                                                                              cATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa               <  4:7795/148‑1 (MQ=255)
                                                                                                                              cATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGa                <  8:111387/147‑1 (MQ=255)
                                                                                                                              cATCAACCAGCCTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa               <  8:186963/148‑1 (MQ=255)
                                                                                                                               aTCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGa             >  6:147854/1‑149 (MQ=255)
                                                                                                                                tCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGAc            <  3:231701/149‑1 (MQ=255)
                                                                                                                                 tcaCCAGCGGATGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACt           >  8:63575/3‑149 (MQ=255)
                                                                                                                                 tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACt           >  1:151370/3‑149 (MQ=255)
                                                                                                                                 tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACt           >  6:87280/3‑149 (MQ=255)
                                                                                                                                 tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTc          >  6:240871/3‑149 (MQ=255)
                                                                                                                                 tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCAGTCCACCGACATGGGTGAGCTGGCTATGGGTCAGCACATGCGCATCGCGTTCATGGCGTTGAACGGCTTCAACTTCGAAGaat           >  6:111873/3‑147 (MQ=255)
                                                                                                                                   aCCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCt                          >  7:81997/1‑132 (MQ=255)
                                                                                                                                   aCCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCt                          <  8:81997/132‑1 (MQ=255)
                                                                                                                                   aCCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa               >  6:34008/1‑143 (MQ=255)
                                                                                                                                   aCCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa               <  5:34008/143‑1 (MQ=255)
                                                                                                                                   aCCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTcc         <  5:147854/149‑1 (MQ=255)
                                                                                                                                    ccAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACGAACATGGCCGACGGCCCGTCCACCGACATGGGTGCGGTGGCTGGGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTc          >  1:293197/1‑148 (MQ=255)
                                                                                                                                     cAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCAt       >  1:390409/1‑149 (MQ=255)
                                                                                                                                     cAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCAt       >  2:217129/1‑149 (MQ=255)
                                                                                                                                          tCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGcc  >  5:281086/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
TGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCC  >  NC_002947/539803‑540089

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: