Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 538,049 A→G 7.5% N143D (AAC→GAC)  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,0490AG 7.5% 230.8 / 6.0 80N143D (AAC→GAC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base G (3/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TGACCTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCCTGGCTGGACTTCGAGTT  >  NC_002947/537903‑538197
                                                                                                                                                  |                                                                                                                                                    
tgACCTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAAc                                                                                                                                                    <  2:41456/149‑1 (MQ=255)
  aCCTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACg                                                                                                                                                   <  5:16369/148‑1 (MQ=255)
    cTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTa                                                                                                                                                >  5:103730/1‑149 (MQ=255)
    cTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTa                                                                                                                                                >  8:192581/1‑149 (MQ=255)
       cgcgGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTAcc                                                                                                                                              <  8:175656/148‑1 (MQ=255)
        gcgGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGa                                                                                                                                            <  7:192581/149‑1 (MQ=255)
              ccACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCgtgt                                                                                                                                      <  3:66730/149‑1 (MQ=255)
              ccACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCgtgt                                                                                                                                      >  1:53564/1‑149 (MQ=255)
               cACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGtgtg                                                                                                                                     <  7:150538/149‑1 (MQ=255)
                aCTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGa                                                                                                                                    >  2:303870/1‑149 (MQ=255)
                aCTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGa                                                                                                                                    >  1:97632/1‑149 (MQ=255)
                    cGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGt                                                                                                                                <  2:97632/149‑1 (MQ=255)
                       gTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTc                                                                                                                             >  1:271667/1‑149 (MQ=255)
                          aaGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTccc                                                                                                                           <  2:53564/148‑1 (MQ=255)
                          aaGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTcc                                                                                                                            >  7:175353/1‑147 (MQ=255)
                          aaGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTcc                                                                                                                            >  8:86609/1‑147 (MQ=255)
                             gTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCt                                                                                                                       <  6:202480/149‑1 (MQ=255)
                                 gcCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATAAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCAc                                                                                                                   >  5:378967/1‑149 (MQ=255)
                                     gATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGtt                                                                                                               >  5:350395/1‑149 (MQ=255)
                                      atcatcTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGt                                                                                                                <  6:103730/147‑1 (MQ=255)
                                      atcatcTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTc                                                                                                              >  6:265613/1‑149 (MQ=255)
                                      atcatcTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTc                                                                                                              <  8:280838/149‑1 (MQ=255)
                                         atcTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCCAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGAt                                                                                                                                   <  1:69722/125‑1 (MQ=255)
                                         atcTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCCAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGAt                                                                                                                                   >  2:69722/1‑125 (MQ=255)
                                              cGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGgt                                                                                                       <  5:265613/148‑1 (MQ=255)
                                              cGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTCgt                                                                                                       <  5:289243/148‑1 (MQ=255)
                                                   aGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCGACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGttc                                                                                                 >  3:261377/1‑149 (MQ=255)
                                                     gAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGttctt                                                                                               >  1:144492/1‑149 (MQ=255)
                                                     gAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGttctt                                                                                               >  7:87188/1‑149 (MQ=255)
                                                      aaTCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGttcttc                                                                                              <  4:281390/149‑1 (MQ=255)
                                                        tcgtcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGgt                                                                                                       <  3:189566/138‑1 (MQ=255)
                                                        tcgtcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGgt                                                                                                       >  4:189566/1‑138 (MQ=255)
                                                        tcgtcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGtt                                                                                                  <  6:354538/143‑1 (MQ=255)
                                                        tcgtcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGtt                                                                                                  >  5:354538/1‑143 (MQ=255)
                                                          gtcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCa                                                                                                                          <  5:360262/117‑1 (MQ=255)
                                                          gtcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCa                                                                                                                          >  6:360262/1‑117 (MQ=255)
                                                          gtcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGgtg                                                                                                      >  1:293474/1‑137 (MQ=255)
                                                           tcgAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCa                                                                                         >  3:11054/1‑149 (MQ=255)
                                                                 aaaGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCg                                                                                   >  3:185010/1‑149 (MQ=255)
                                                                 aaaGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCg                                                                                   >  1:277726/1‑149 (MQ=255)
                                                                  aaGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGt                                                                                  <  4:11054/149‑1 (MQ=255)
                                                                          aaGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGAc                                                                          >  6:181810/1‑149 (MQ=255)
                                                                                     agagCAAGAAGTCTACATGGGTTAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCgg                                                               <  6:350395/149‑1 (MQ=255)
                                                                                      gagCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCggg                                                              >  1:23428/1‑149 (MQ=255)
                                                                                             aagTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGAcc                                                                                    <  6:176660/120‑1 (MQ=255)
                                                                                             aagTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGAcc                                                                                    >  5:176660/1‑120 (MQ=255)
                                                                                             aagTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAgctg                                                       <  5:181810/149‑1 (MQ=255)
                                                                                             aagTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAgctg                                                       <  2:277726/149‑1 (MQ=255)
                                                                                              agTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCGACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAgctg                                                       <  4:261377/148‑1 (MQ=255)
                                                                                                 cTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCCCAGCTCGGGCAAGctgctg                                                    >  7:199780/1‑148 (MQ=255)
                                                                                                  tACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGt                                                   >  3:332707/1‑148 (MQ=255)
                                                                                                    cATGGGGGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCGGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACt                                                <  6:105704/149‑1 (MQ=255)
                                                                                                            aaaTCCCCCTGATGACTGAGAATGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCg                                         >  3:49587/1‑148 (MQ=255)
                                                                                                                   ccTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAgctgct                                                     <  3:178689/129‑1 (MQ=255)
                                                                                                                   ccTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAgctgct                                                     >  4:178689/1‑129 (MQ=255)
                                                                                                                   ccTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATcc                                 <  4:185010/149‑1 (MQ=255)
                                                                                                                    cTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATccc                                >  8:279884/1‑149 (MQ=255)
                                                                                                                     tgatgaCTGAGAACGGTACCTTCGTTATCAACGGTACCGAGGGTGTGATGGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCcac                                                                      <  6:30886/110‑1 (MQ=255)
                                                                                                                     tgatgaCTGAGAACGGTACCATCGTTATCAACGGTACCGAGCGTGAGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCcac                                                                      >  5:30886/1‑110 (MQ=255)
                                                                                                                         gaCTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGAATCATCCCTTAcc                           >  5:203398/1‑149 (MQ=255)
                                                                                                                          aCTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTAcc                           <  7:86609/148‑1 (MQ=255)
                                                                                                                              agaATGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCgg                       <  4:49587/148‑1 (MQ=255)
                                                                                                                              agaACGGTACGTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCgg                       <  2:271667/148‑1 (MQ=255)
                                                                                                                                 tcGGTACCGTTGATAACGAAGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGtt                     <  6:252788/146‑1 (MQ=255)
                                                                                                                                 tcGGTACCGTTGATAACGAAGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGtt                     >  5:252788/2‑147 (MQ=255)
                                                                                                                                 tcGGTACCGTTGATAACGAAGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGt                      >  6:351844/2‑146 (MQ=255)
                                                                                                                                 tcGGTACCGTTGATAACGAAGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGt                      <  5:351844/145‑1 (MQ=255)
                                                                                                                                  cGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCct                  <  2:23428/149‑1 (MQ=255)
                                                                                                                                  cGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCct                  >  4:214819/1‑149 (MQ=255)
                                                                                                                                   ggTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCctg                 >  7:167353/1‑149 (MQ=255)
                                                                                                                                      aCCTTCGTTATCAACGGTACCGCGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCAGCCTCCCTGACCGCGGTTCctggct              >  6:9165/1‑149 (MQ=255)
                                                                                                                                      aCCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCctggct              >  6:6724/1‑149 (MQ=255)
                                                                                                                                        cTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCctggctg             >  8:132237/1‑148 (MQ=255)
                                                                                                                                        cTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCctggctg             >  1:341413/1‑148 (MQ=255)
                                                                                                                                          tCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTc                                          <  6:352458/117‑1 (MQ=255)
                                                                                                                                          tCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTc                                          >  5:352458/1‑117 (MQ=255)
                                                                                                                                             ttATCAACGGGACCGAGCCGGTGATCGTTTCCCAGCTGCACCGTGCGCCTGGTGTGGTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCCTGGCTGGACTTc       <  8:175353/149‑1 (MQ=255)
                                                                                                                                               aTCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCCTGGCTGGACTTc       >  2:337749/1‑147 (MQ=255)
                                                                                                                                                tCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCct                  <  7:212353/135‑1 (MQ=255)
                                                                                                                                                tCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCct                  >  8:212353/1‑135 (MQ=255)
                                                                                                                                                  aaCGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACAGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCCTGGCTGGACTTCGAGtt  >  6:140708/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGACCTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCCTGGCTGGACTTCGAGTT  >  NC_002947/537903‑538197

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: ‑