Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,070,657 T→G 8.6% S125R (AGC→CGC)  mdtC ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,070,6570TG 8.6% 537.6 / 27.3 168S125R (AGC→CGC) mdtCmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base T (87/66);  new base G (10/5);  total (97/71)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

CTCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCG  >  NC_002947/4070508‑4070805
                                                                                                                                                     |                                                                                                                                                    
cTCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCAtgct                                                                                                                                                      <  1:5072700/150‑1 (MQ=255)
 tCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCAtgctg                                                                                                                                                     >  1:3020619/1‑150 (MQ=255)
   cccTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTgcg                                                                                                                                                   >  1:2534778/1‑150 (MQ=255)
        cccccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGc                                                                                                                                              <  2:880539/150‑1 (MQ=255)
        cccccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGc                                                                                                                                              <  2:3882980/150‑1 (MQ=255)
         ccccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGcc                                                                                                                                             >  2:620145/1‑150 (MQ=255)
           ccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCAc                                                                                                                                           >  1:4312807/1‑150 (MQ=255)
             gACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTg                                                                                                                                         <  2:2023884/150‑1 (MQ=255)
               cacCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGCTCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTggg                                                                                                                                       <  2:1221992/150‑1 (MQ=255)
                acCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGTGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgc                                                                                                                                      >  2:2100493/1‑150 (MQ=255)
                acCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgc                                                                                                                                      <  1:5021619/150‑1 (MQ=255)
                acCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgc                                                                                                                                      <  2:4655363/150‑1 (MQ=255)
                 ccTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgca                                                                                                                                     <  2:2931431/150‑1 (MQ=255)
                  cTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcag                                                                                                                                    >  1:1865256/1‑150 (MQ=255)
                  cTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcag                                                                                                                                    >  1:530402/1‑150 (MQ=255)
                   tGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagc                                                                                                                                   >  1:3630100/1‑150 (MQ=255)
                       ccAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGt                                                                                                                               >  1:3057555/1‑150 (MQ=255)
                          ggCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGc                                                                                                                            >  1:523292/1‑150 (MQ=255)
                          ggCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGc                                                                                                                            >  1:523409/1‑150 (MQ=255)
                            ctctGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCgg                                                                                                                          <  1:1941159/150‑1 (MQ=255)
                             tctGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCggg                                                                                                                         >  2:1944064/1‑150 (MQ=255)
                               tgtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTa                                                                                                                       <  1:4757999/150‑1 (MQ=255)
                               tgtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTa                                                                                                                       >  2:1127030/1‑150 (MQ=255)
                                gtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAg                                                                                                                      <  1:1328236/150‑1 (MQ=255)
                                 tgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGc                                                                                                                     >  1:4523083/1‑150 (MQ=255)
                                 tgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGc                                                                                                                     <  2:4523083/150‑1 (MQ=255)
                                   aaaGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGt                                                                                                                   >  1:3044094/1‑150 (MQ=255)
                                      ggATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGa                                                                                                                <  1:3124526/150‑1 (MQ=255)
                                        aTGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGCATGCGCAGCATGcccaccta                                                                                                                                       >  2:399125/1‑118 (MQ=255)
                                             gtCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATcgccgc                                                                                                         >  2:1319514/1‑150 (MQ=255)
                                                ggCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTg                                                                                                      >  1:2690609/1‑150 (MQ=255)
                                                ggCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTg                                                                                                      <  2:3868940/150‑1 (MQ=255)
                                                 gCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGa                                                                                                     >  1:1035374/1‑150 (MQ=255)
                                                  ccAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  2:1887744/150‑1 (MQ=255)
                                                  ccAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  1:1039919/150‑1 (MQ=255)
                                                      gTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagca                                                                                                                                  <  2:3088950/116‑1 (MQ=255)
                                                      gTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagca                                                                                                                                  >  1:3088950/1‑116 (MQ=255)
                                                         aTACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACg                                                                                             >  2:2965124/1‑150 (MQ=255)
                                                         aTACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACg                                                                                             >  2:2965070/1‑150 (MQ=255)
                                                              gCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTg                                                                                        >  2:4021467/1‑150 (MQ=255)
                                                                  gCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGcc                                                                                    >  1:4193167/1‑150 (MQ=255)
                                                                        cTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGAt                                                                              <  2:412604/150‑1 (MQ=255)
                                                                        cTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGAt                                                                              <  2:3494378/150‑1 (MQ=255)
                                                                        cTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGAt                                                                              <  2:1974747/150‑1 (MQ=255)
                                                                         tGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATg                                                                             <  1:4492063/150‑1 (MQ=255)
                                                                               acaTCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACgg                                                                       <  1:353787/150‑1 (MQ=255)
                                                                               acaTCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCGCACGGGCTGCGCCATCGATGTCACgg                                                                       >  2:1655079/1‑150 (MQ=255)
                                                                                 tcgtggtgTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGcc                                                                     >  2:4158045/8‑150 (MQ=255)
                                                                                  tCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGccc                                                                    >  2:2570202/1‑150 (MQ=255)
                                                                                   ccGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagc                                                                                                                                   <  1:3256992/86‑1 (MQ=255)
                                                                                   ccGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCa                                                                   >  2:2019078/1‑150 (MQ=255)
                                                                                    cGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGAATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCaa                                                                  >  1:910873/1‑150 (MQ=255)
                                                                                     gAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAAc                                                                 >  1:1948296/1‑150 (MQ=255)
                                                                                         tCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    >  2:5274477/1‑111 (MQ=255)
                                                                                         tCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  1:5274477/111‑1 (MQ=255)
                                                                                         tCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGACCCAACTCGa                                                             >  2:822612/1‑150 (MQ=255)
                                                                                           aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAg                                                           <  2:4412488/150‑1 (MQ=255)
                                                                                                cAGCACCCTGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  2:3258583/150‑1 (MQ=255)
                                                                                                cAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  2:2902162/150‑1 (MQ=255)
                                                                                                    aCCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACt                                                  >  1:265486/1‑150 (MQ=255)
                                                                                                      cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  1:482369/1‑150 (MQ=255)
                                                                                                      cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  2:520495/1‑150 (MQ=255)
                                                                                                      cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                <  2:4178213/150‑1 (MQ=255)
                                                                                                      cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  1:2944628/1‑150 (MQ=255)
                                                                                                      cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  1:849308/1‑150 (MQ=255)
                                                                                                       aTGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAcc                                               >  1:3281976/1‑150 (MQ=255)
                                                                                                       aTGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAcc                                               >  1:2129148/1‑150 (MQ=255)
                                                                                                        tGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAccc                                              <  1:454611/150‑1 (MQ=255)
                                                                                                           tCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  2:2494749/93‑1 (MQ=255)
                                                                                                              gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGAc                                                   >  1:3787055/1‑139 (MQ=255)
                                                                                                              gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGAc                                                   <  2:3787055/139‑1 (MQ=255)
                                                                                                              gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCg                                        >  2:368552/1‑150 (MQ=255)
                                                                                                                 ccTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAAc                                     >  2:1721873/1‑150 (MQ=255)
                                                                                                                   tGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAAccc                                   <  1:2430802/150‑1 (MQ=255)
                                                                                                                        cGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGtt                              <  2:2542764/150‑1 (MQ=255)
                                                                                                                           gTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           >  2:519500/1‑150 (MQ=255)
                                                                                                                                tcttctTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAgctgc                      >  1:4364134/1‑150 (MQ=255)
                                                                                                                                tcttctTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAgctgc                      >  2:3252003/1‑150 (MQ=255)
                                                                                                                                tcgtcttGTAGGTGGGCATGCGGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           <  1:2058740/142‑1 (MQ=255)
                                                                                                                                 cttcttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAgctgct                     >  1:4160328/1‑150 (MQ=255)
                                                                                                                                   tcttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTgg                   >  1:4087704/1‑150 (MQ=255)
                                                                                                                                   tcttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTgg                   <  1:368552/150‑1 (MQ=255)
                                                                                                                                    cttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGt                  <  2:482369/150‑1 (MQ=255)
                                                                                                                                     ttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTg                 >  2:101034/1‑150 (MQ=255)
                                                                                                                                      tGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCTCCATCGATGTCCCGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGa                <  1:4789394/150‑1 (MQ=255)
                                                                                                                                       gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATcgccgcaac                                                                                                      <  1:166636/63‑4 (MQ=255)
                                                                                                                                       gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCt                                                                                         <  1:3934418/76‑1 (MQ=255)
                                                                                                                                         ccagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            >  2:5394973/6‑103 (MQ=255)
                                                                                                                                         ccagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            <  1:5394973/98‑1 (MQ=255)
                                                                                                                                          ggTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  1:4070808/1‑150 (MQ=255)
                                                                                                                                          cagtggCATGCGCAGCATGCCGCTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  1:841344/5‑150 (MQ=255)
                                                                                                                                          cagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  2:4604100/5‑150 (MQ=255)
                                                                                                                                             tggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         >  2:4608467/2‑150 (MQ=255)
                                                                                                                                             tggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         >  1:2591005/2‑150 (MQ=255)
                                                                                                                                             gggCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         <  2:530402/150‑1 (MQ=255)
                                                                                                                                             gggCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgtggt         <  1:728792/150‑1 (MQ=255)
                                                                                                                                               gCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCCGCTGGTGAGCGTGGTGa       >  2:4471013/1‑150 (MQ=255)
                                                                                                                                                cATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGAc      >  1:2520980/1‑150 (MQ=255)
                                                                                                                                                cATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGAc      >  1:4426703/1‑150 (MQ=255)
                                                                                                                                                    ctgcGCATGCCACTGGGCAACAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCg  >  1:1564938/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCG  >  NC_002947/4070508‑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: