Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 538,397 A→G 6.1% I259V (ATC→GTC)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947538,3970AG 6.1% 285.2 / 6.6 98I259V (ATC→GTC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (49/43);  new base G (2/4);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGT  >  NC_002947/538250‑538542
                                                                                                                                                   |                                                                                                                                                 
gCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACAt                                                                                                                                                  >  4:405025/1‑149 (MQ=255)
       tgctgctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGAc                                                                                                                                           >  5:47526/1‑149 (MQ=255)
         ctgctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGa                                                                                                                                         <  4:207119/149‑1 (MQ=255)
           gctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgca                                                                                                                                       <  4:305506/149‑5 (MQ=255)
            ctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgc                                                                                                                                        <  5:2922/147‑4 (MQ=255)
            ctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACAAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAAc                                                                                                                                      >  3:379876/1‑149 (MQ=255)
               cgcgcCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCgg                                                                                                                                   >  8:323398/1‑149 (MQ=255)
               cgcgcCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACAAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAAcc                                                                                                                                     <  4:379876/147‑1 (MQ=255)
                   cccTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCaaa                                                                                                                               <  4:60023/149‑1 (MQ=255)
                   cccTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCaaa                                                                                                                               <  3:416311/149‑1 (MQ=255)
                    ccTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAg                                                                                                                              <  1:54093/149‑1 (MQ=255)
                        ggTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCa                                                                                                                           <  6:365884/148‑1 (MQ=255)
                         gTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtc                                                                                                                         <  5:3854/149‑1 (MQ=255)
                             cAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtcgtcg                                                                                                                     <  6:418994/149‑1 (MQ=255)
                             cAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtcgtc                                                                                                                      <  7:175352/148‑1 (MQ=255)
                             cAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtcgtc                                                                                                                      <  6:94568/148‑1 (MQ=255)
                                cACTGAAGAAGTGTTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAg                                                                                                                   >  1:157313/1‑148 (MQ=255)
                                 aCTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTACAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGc                                                                                                                  >  8:399098/1‑148 (MQ=255)
                                  cTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGCCGACACCGGCAAAGTCATCGTGGAGc                                                                                                                  >  4:342203/1‑147 (MQ=255)
                                  cTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGc                                                                                                                  >  4:17569/1‑147 (MQ=255)
                                       gaagTGCTGAACACCTTCTACACAACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGgccg                                                                                                           >  4:354378/1‑149 (MQ=255)
                                        aagTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgc                                                                                                                                        <  5:359185/119‑4 (MQ=255)
                                        aagTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgc                                                                                                                                        >  6:359185/1‑116 (MQ=255)
                                          gTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGccgccg                                                                                                        >  2:225091/1‑149 (MQ=255)
                                           tGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAAc                                                                                                                                      >  7:210797/1‑118 (MQ=255)
                                           tGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAAc                                                                                                                                      <  8:210797/118‑1 (MQ=255)
                                            gCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGa                                                                                                                                            <  4:368762/111‑1 (MQ=255)
                                            gCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGa                                                                                                                                            >  3:368762/1‑111 (MQ=255)
                                             cTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTAt                                                                                                     <  3:239193/149‑1 (MQ=255)
                                             cTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTAt                                                                                                     >  6:146513/1‑149 (MQ=255)
                                                aaCACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTAc                                                                                                  >  8:202517/1‑149 (MQ=255)
                                                  cacCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTa                                                                                                   >  2:112963/1‑146 (MQ=255)
                                                  cacCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTa                                                                                                   <  1:112963/146‑1 (MQ=255)
                                                   acCTTCTACACCACCAACGTGTTCCACATTTCCGTCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCAGCCGTATTACcgc                                                                                               >  5:62607/1‑149 (MQ=255)
                                                     cTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgc                                                                                             >  6:311711/1‑149 (MQ=255)
                                                     cTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgc                                                                                             >  3:300605/1‑149 (MQ=255)
                                                        cTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCgcgcgc                                                                                           >  7:111146/1‑148 (MQ=255)
                                                        cTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCgcgcgc                                                                                           <  5:311711/148‑1 (MQ=255)
                                                        cTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCgcgcgc                                                                                           <  7:84199/148‑1 (MQ=255)
                                                        cTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCgcgcgc                                                                                           <  4:300605/148‑1 (MQ=255)
                                                            accaccAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACAt                                                                                      >  7:70903/1‑149 (MQ=255)
                                                              caccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCgcgcgc                                                                                           <  6:322536/142‑1 (MQ=255)
                                                              caccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCgcgcgc                                                                                           >  5:322536/1‑142 (MQ=255)
                                                              caccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATc                                                                                     >  8:332172/1‑148 (MQ=255)
                                                                ccaACGTGTTCCACATTTCCGGCGTAAAACTCAGCCTGTAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAAc                                                                                  <  6:62607/149‑1 (MQ=255)
                                                                 caACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAAcc                                                                                 <  8:111146/149‑1 (MQ=255)
                                                                 caACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAAc                                                                                  >  6:135470/1‑148 (MQ=255)
                                                                        ttCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGa                                                                                                                    >  6:345740/1‑107 (MQ=255)
                                                                          ccACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAATGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCgag                                                                         >  8:74054/1‑148 (MQ=255)
                                                                          ccACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCgag                                                                         >  3:97147/1‑148 (MQ=255)
                                                                          ccACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCgag                                                                         >  6:218965/1‑148 (MQ=255)
                                                                           cacaTTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTa                                                                                                   >  2:151365/1‑121 (MQ=255)
                                                                           cacaTTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGaa                                                                       <  5:135470/149‑1 (MQ=255)
                                                                              aTTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgc                                                                                               <  2:152662/122‑1 (MQ=255)
                                                                              aTTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgc                                                                                               >  1:152662/1‑122 (MQ=255)
                                                                                 tCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCg                                                                  <  8:294875/148‑1 (MQ=255)
                                                                                 tCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGGCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCg                                                                  <  4:158342/148‑1 (MQ=255)
                                                                                   cGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGc                                                                <  5:105767/148‑1 (MQ=255)
                                                                                       gAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCaa                                                           <  7:326222/149‑1 (MQ=255)
                                                                                       gAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCaa                                                           <  7:323398/149‑1 (MQ=255)
                                                                                        aaaaaCTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAg                                                          >  3:175874/1‑149 (MQ=255)
                                                                                            aCTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGAAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAAc                                                      >  7:288775/1‑149 (MQ=255)
                                                                                                  ccTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACg                                                >  3:60978/1‑149 (MQ=255)
                                                                                                   cTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGt                                               <  3:354378/149‑1 (MQ=255)
                                                                                                   cTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACAAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGt                                               >  7:155322/1‑149 (MQ=255)
                                                                                                      gAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAATGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTc                                             <  7:74054/148‑1 (MQ=255)
                                                                                                      gAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTc                                             <  3:342203/148‑1 (MQ=255)
                                                                                                             tACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGaa                                     <  8:155322/149‑1 (MQ=255)
                                                                                                             tACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGaa                                     <  6:47526/149‑1 (MQ=255)
                                                                                                             tACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGaa                                     <  7:347571/149‑1 (MQ=255)
                                                                                                               ccGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATa                                   <  4:97147/149‑1 (MQ=255)
                                                                                                                cGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATAc                                  >  6:2974/1‑149 (MQ=255)
                                                                                                                cGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATAc                                  <  7:332172/149‑1 (MQ=255)
                                                                                                                 gCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACg                                 >  4:312354/1‑149 (MQ=255)
                                                                                                                         tGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCCCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGc                         >  4:331180/1‑149 (MQ=255)
                                                                                                                          gCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGcc                        >  8:305349/1‑149 (MQ=255)
                                                                                                                           cGTGGTGAAGTTGCAGTCATGGACATCCATGACTAAACCGGCAAAGTCAtcg                                                                                                                        <  3:23178/52‑1 (MQ=255)
                                                                                                                           cGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtcg                                                                                                                        >  4:23178/1‑52 (MQ=255)
                                                                                                                           cGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCg                       <  8:70903/149‑1 (MQ=255)
                                                                                                                                 gAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTAc                 >  2:137806/1‑149 (MQ=255)
                                                                                                                                  aaGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTAcc                >  8:122911/1‑149 (MQ=255)
                                                                                                                                  aaGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTAcc                >  6:116858/1‑149 (MQ=255)
                                                                                                                                   aGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCTCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCg               <  5:116858/149‑1 (MQ=255)
                                                                                                                                    gTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGc              >  4:41121/1‑149 (MQ=255)
                                                                                                                                       gCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGcc             >  2:345899/1‑147 (MQ=255)
                                                                                                                                        cAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGGCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAg          <  3:17569/149‑1 (MQ=255)
                                                                                                                                         aGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCTTAATAACCCGCGCCACATAAACCAGCTCGAGAACGCCGGCGTCAAGCAACTGGACGTTACAATGGAAAACGTCCTGGGCCCCACTAACGCCAAgg         >  1:115049/1‑149 (MQ=255)
                                                                                                                                         aGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGCATACGTCCTGGGCCGCACTACCGCCAAgg         <  1:32981/149‑1 (MQ=255)
                                                                                                                                         aGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAgg         >  2:397722/1‑149 (MQ=255)
                                                                                                                                          gTCATGGATGTCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCa                     >  8:118928/1‑136 (MQ=255)
                                                                                                                                          gTCATGGATGTCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCa                     <  7:118928/136‑1 (MQ=255)
                                                                                                                                          gTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGc        <  5:2974/149‑1 (MQ=255)
                                                                                                                                           tCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGc        >  7:168388/1‑148 (MQ=255)
                                                                                                                                           tCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGc        >  5:192983/1‑148 (MQ=255)
                                                                                                                                             aTGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTAt     >  3:326949/1‑149 (MQ=255)
                                                                                                                                               ggACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCg   >  3:358640/1‑149 (MQ=255)
                                                                                                                                                gACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGTCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCg   <  5:146513/148‑1 (MQ=255)
                                                                                                                                                gACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGt  <  3:312354/149‑1 (MQ=255)
                                                                                                                                                   aTCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTc                                             <  4:167886/103‑1 (MQ=255)
                                                                                                                                                   aTCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTc                                             >  3:167886/1‑103 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGT  >  NC_002947/538250‑538542

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: