Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,070,659 C→A 7.6% R124L (CGC→CTC) ‡ mdtC ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,070,6590CA 7.6% 543.0 / 23.1 166R124L (CGC→CTC) ‡mdtCmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base C (85/68);  new base A (10/3);  total (95/71)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCG  >  NC_002947/4070511‑4070805
                                                                                                                                                    |                                                                                                                                                  
cccTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTgcg                                                                                                                                                   >  1:2534778/1‑150 (MQ=255)
  caaCCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCa                                                                                                                                                 >  1:3355740/3‑150 (MQ=255)
  cTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCa                                                                                                                                                 <  1:3615156/150‑1 (MQ=255)
     cccccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGc                                                                                                                                              <  2:880539/150‑1 (MQ=255)
      ccccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGcc                                                                                                                                             >  2:620145/1‑150 (MQ=255)
        ccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCAc                                                                                                                                           >  1:4312807/1‑150 (MQ=255)
         cGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACt                                                                                                                                          <  2:1167563/150‑1 (MQ=255)
         cGACACCTGAGCCAGGCTCTGCGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACt                                                                                                                                          >  2:1442968/1‑150 (MQ=255)
          gACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTg                                                                                                                                         <  2:2023884/150‑1 (MQ=255)
          gACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTg                                                                                                                                         <  2:2436586/150‑1 (MQ=255)
             acCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgc                                                                                                                                      <  2:4655363/150‑1 (MQ=255)
             acCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgc                                                                                                                                      <  1:5021433/150‑1 (MQ=255)
             acCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgc                                                                                                                                      <  1:5021619/150‑1 (MQ=255)
              ccTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgca                                                                                                                                     <  2:1099680/150‑1 (MQ=255)
              ccTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgca                                                                                                                                     <  2:2931431/150‑1 (MQ=255)
               cTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcag                                                                                                                                    >  2:5224856/1‑150 (MQ=255)
               cTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcag                                                                                                                                    >  1:1865256/1‑150 (MQ=255)
                tGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagc                                                                                                                                   >  1:3630100/1‑150 (MQ=255)
                    ccAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGt                                                                                                                               >  1:2622023/1‑150 (MQ=255)
                       ggCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGc                                                                                                                            >  1:523409/1‑150 (MQ=255)
                         ctctGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCgg                                                                                                                          <  1:1941159/150‑1 (MQ=255)
                            tgtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTa                                                                                                                       <  1:4757999/150‑1 (MQ=255)
                             gtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAg                                                                                                                      <  1:1328236/150‑1 (MQ=255)
                              tgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGc                                                                                                                     <  2:4523083/150‑1 (MQ=255)
                              tgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGc                                                                                                                     >  1:4523083/1‑150 (MQ=255)
                                aaaGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGt                                                                                                                   >  1:3044094/1‑150 (MQ=255)
                                aaaGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGt                                                                                                                   >  1:3043645/1‑150 (MQ=255)
                                 aaGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGtt                                                                                                                  <  1:3903340/150‑1 (MQ=255)
                                   ggATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGa                                                                                                                <  1:3124526/150‑1 (MQ=255)
                                    gATGGTGTGGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGAt                                                                                                               >  2:4492063/1‑150 (MQ=255)
                                     aTGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGCATGCGCAGCATGcccaccta                                                                                                                                       >  2:399125/1‑118 (MQ=255)
                                          gtCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATcgccgc                                                                                                         >  2:1319514/1‑150 (MQ=255)
                                             ggCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTg                                                                                                      >  1:2690609/1‑150 (MQ=255)
                                              gCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGa                                                                                                     >  1:5388024/1‑150 (MQ=255)
                                              gCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGa                                                                                                     >  1:1035374/1‑150 (MQ=255)
                                               ccAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  2:1887744/150‑1 (MQ=255)
                                               ccAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  1:1039919/150‑1 (MQ=255)
                                                      aTACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACg                                                                                             >  2:2965070/1‑150 (MQ=255)
                                                      aTACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACg                                                                                             >  2:197191/1‑150 (MQ=255)
                                                           gCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTg                                                                                        >  2:4021467/1‑150 (MQ=255)
                                                               gCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGcc                                                                                    >  1:4193167/1‑150 (MQ=255)
                                                                 cTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCAt                                                                                  >  1:3093012/1‑150 (MQ=255)
                                                                      tGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATg                                                                             <  1:4492063/150‑1 (MQ=255)
                                                                              tcgtggtgTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGcc                                                                     >  2:4158045/8‑150 (MQ=255)
                                                                              tcgtggtgTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGcc                                                                     >  2:4158034/8‑150 (MQ=255)
                                                                                ccGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagc                                                                                                                                   <  1:3256992/86‑1 (MQ=255)
                                                                                 cGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCaa                                                                  <  2:3630100/150‑1 (MQ=255)
                                                                                 cGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGAATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCaa                                                                  >  1:910873/1‑150 (MQ=255)
                                                                                   aaGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACt                                                                <  2:2929054/150‑1 (MQ=255)
                                                                                      tCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    >  2:5274477/1‑111 (MQ=255)
                                                                                      tCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGACCCAACTCGa                                                             >  2:822612/1‑150 (MQ=255)
                                                                                       taGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            <  2:198823/149‑1 (MQ=255)
                                                                                       cAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            <  1:77240/150‑1 (MQ=255)
                                                                                        aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGa                                                             >  1:3067246/1‑148 (MQ=255)
                                                                                        aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGa                                                             <  2:3067246/148‑1 (MQ=255)
                                                                                        aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAg                                                           <  2:4412488/150‑1 (MQ=255)
                                                                                        aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAg                                                           >  2:261776/1‑150 (MQ=255)
                                                                                             cAGCACCCTGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  2:3258583/150‑1 (MQ=255)
                                                                                             cAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  1:3687383/150‑1 (MQ=255)
                                                                                             cAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  1:945324/150‑1 (MQ=255)
                                                                                                   cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                <  2:4178213/150‑1 (MQ=255)
                                                                                                   cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  1:2944628/1‑150 (MQ=255)
                                                                                                    aTGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAcc                                               >  1:2129148/1‑150 (MQ=255)
                                                                                                     tGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAccc                                              <  1:454611/150‑1 (MQ=255)
                                                                                                        tCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    >  1:2494749/1‑93 (MQ=255)
                                                                                                           gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGAc                                                   <  2:3787055/139‑1 (MQ=255)
                                                                                                           gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCg                                        >  2:368552/1‑150 (MQ=255)
                                                                                                                tGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAAccc                                   <  1:2430802/150‑1 (MQ=255)
                                                                                                                   gACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTgg                                <  2:3153470/150‑1 (MQ=255)
                                                                                                                     cGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGtt                              <  2:2542764/150‑1 (MQ=255)
                                                                                                                        gTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           >  2:519500/1‑150 (MQ=255)
                                                                                                                        gTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           <  1:5245174/150‑1 (MQ=255)
                                                                                                                             tcttctTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           >  2:2058740/1‑145 (MQ=255)
                                                                                                                             tcttctTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAgctgc                      >  2:3252003/1‑150 (MQ=255)
                                                                                                                             tcgtcttGTAGGTGGGCATGCGGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           <  1:2058740/142‑1 (MQ=255)
                                                                                                                              cttcttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAgctgct                     >  1:4160174/1‑150 (MQ=255)
                                                                                                                                 cttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGt                  <  2:482369/150‑1 (MQ=255)
                                                                                                                                   tGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCTCCATCGATGTCCCGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGa                <  1:4789394/150‑1 (MQ=255)
                                                                                                                                   tGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGa                <  1:5281922/150‑1 (MQ=255)
                                                                                                                                    gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATcgccgcaac                                                                                                      >  2:166636/1‑60 (MQ=255)
                                                                                                                                    gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATcgccgcaac                                                                                                      <  1:166636/63‑4 (MQ=255)
                                                                                                                                    gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCt                                                                                         <  1:3934418/76‑1 (MQ=255)
                                                                                                                                    gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCt                                                                                         >  2:3934418/1‑76 (MQ=255)
                                                                                                                                      ccagtggCATGCGCAGCATGCCCCTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            <  1:4161449/98‑1 (MQ=255)
                                                                                                                                      ccagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            >  2:5394973/6‑103 (MQ=255)
                                                                                                                                       ggTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  1:4070808/1‑150 (MQ=255)
                                                                                                                                       cagtggCATGCGCAGCATGCCGCTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  1:841608/5‑150 (MQ=255)
                                                                                                                                       cagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  1:4381836/5‑150 (MQ=255)
                                                                                                                                       cagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  2:4604100/5‑150 (MQ=255)
                                                                                                                                        agtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtg           >  2:3873025/4‑150 (MQ=255)
                                                                                                                                          tggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         >  1:2591005/2‑150 (MQ=255)
                                                                                                                                          gggCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         <  2:4645289/150‑1 (MQ=255)
                                                                                                                                          gggCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         <  2:1865256/150‑1 (MQ=255)
                                                                                                                                          gggCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         <  1:728792/150‑1 (MQ=255)
                                                                                                                                          gggCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         <  2:530402/150‑1 (MQ=255)
                                                                                                                                           ggCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggtg        >  2:2525165/1‑150 (MQ=255)
                                                                                                                                             cATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGAc      >  1:1942695/1‑150 (MQ=255)
                                                                                                                                             cATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGAc      >  1:5157359/1‑150 (MQ=255)
                                                                                                                                             cATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGAc      >  1:4426703/1‑150 (MQ=255)
                                                                                                                                                 ctgcGCATGCCACTGGGCAACAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCg  >  1:1564938/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCG  >  NC_002947/4070511‑4070805

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