Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 538,047 T→A 5.1% I142N (ATC→AAC)  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,0470TA 5.1% 248.1 / 2.4 78I142N (ATC→AAC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base A (2/2);  total (42/36)
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.94e-01

TGACCTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCCTGGCTGGACTTC  >  NC_002947/537903‑538192
                                                                                                                                                |                                                                                                                                                 
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                                                  <  2:277726/149‑1 (MQ=255)
                                                                                             aagTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAgctg                                                  <  5:181810/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                >  5:252788/2‑147 (MQ=255)
                                                                                                                                 tcGGTACCGTTGATAACGAAGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGtt                <  6:252788/146‑1 (MQ=255)
                                                                                                                                 tcGGTACCGTTGATAACGAAGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGt                 <  5:351844/145‑1 (MQ=255)
                                                                                                                                 tcGGTACCGTTGATAACGAAGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGt                 >  6:351844/2‑146 (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)
                                                                                                                                                |                                                                                                                                                 
TGACCTTCGCGGTCCCACTGCGGGTCAAGGTGCGCCTGATCATCTTCGACAAGGAATCGTCGAACAAAGCGATCAAGGACATCAAAGAGCAAGAAGTCTACATGGGTGAAATCCCCCTGATGACTGAGAACGGTACCTTCGTTATCAACGGTACCGAGCGTGTGATCGTTTCCCAGCTGCACCGTTCGCCTGGTGTGTTCTTCGACCACGACCGTGGCAAGACCCACAGCTCGGGCAAGCTGCTGTACTCCGCTCGCATCATCCCTTACCGCGGTTCCTGGCTGGACTTC  >  NC_002947/537903‑538192

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