Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 539,941 T→G 6.9% R773R (CGT→CGG 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,9410TG 6.9% 305.8 / 12.3 104R773R (CGT→CGGrpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (49/47);  new base G (5/2);  total (54/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 = 5.27e-01

Only 100 of 106 total aligned reads displayed.

GCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCC  >  NC_002947/539793‑540089
                                                                                                                                                    |                                                                                                                                                     
gccgcGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGt                                                                                                                                                       <  6:178038/149‑1 (MQ=255)
          tgtgATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGGACGCTGGttg                                                                                                                                             <  7:205796/149‑3 (MQ=255)
            tgATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGGACGCTGGttgat                                                                                                                                           <  4:101818/149‑5 (MQ=255)
             gATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGc                                                                                                                                           >  7:303991/1‑148 (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:185617/1‑148 (MQ=255)
                      ggTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGAc                                                                                                                                  >  1:179553/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)
                                   cGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtgg                                                                                                                    <  6:76695/149‑1 (MQ=255)
                                   cGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtgg                                                                                                                    <  3:113436/149‑1 (MQ=255)
                                    gTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  2:185617/149‑1 (MQ=255)
                                    gTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  2:179553/149‑1 (MQ=255)
                                    gTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  5:380368/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                                                                                                >  5:311209/1‑149 (MQ=255)
                                                      ccGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGc                                                                                                <  4:34806/149‑1 (MQ=255)
                                                           gacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTc                                                                                           >  5:345977/1‑149 (MQ=255)
                                                            acgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTcc                                                                                          <  2:167864/149‑1 (MQ=255)
                                                              gaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   <  7:5884/123‑1 (MQ=255)
                                                              gaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                                   >  8:5884/1‑123 (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)
                                                                                gCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGc                                                                      <  7:83362/149‑1 (MQ=255)
                                                                                   ggTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGATCGCCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGctct                                                                   <  5:17058/149‑1 (MQ=255)
                                                                                     tgtgGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTgg                                                                 <  7:167498/149‑1 (MQ=255)
                                                                                        ggATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGt                                                               >  3:1283/1‑148 (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                                                       >  5:61075/1‑149 (MQ=255)
                                                                                               tACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACAt                                                       >  3:76681/1‑149 (MQ=255)
                                                                                                aCAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATg                                                      >  7:111387/1‑149 (MQ=255)
                                                                                                aCAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATg                                                      >  3:7795/1‑149 (MQ=255)
                                                                                                aCAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATg                                                      >  6:159858/1‑149 (MQ=255)
                                                                                                 cAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgc                                                     <  8:374649/149‑1 (MQ=255)
                                                                                                 cAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgc                                                     <  2:310525/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                           <  7:271677/149‑1 (MQ=255)
                                                                                                                           aCCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGc                           <  1:91956/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)
                                                                                                                                              ccAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACGAACATGGCCGACGGCCCGTCCACCGACATGGGTGCGGTGGCTGGGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTc          >  1:293197/1‑148 (MQ=255)
                                                                                                                                               cAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCAt       >  2:217129/1‑149 (MQ=255)
                                                                                                                                               cAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCAt       >  1:390409/1‑149 (MQ=255)
                                                                                                                                                    tCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGcc  >  5:281086/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCC  >  NC_002947/539793‑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: