Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,070,658 G→C 7.8% R124R (CGC→CGG) ‡ mdtC ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,070,6580GC 7.8% 541.6 / 26.9 166R124R (CGC→CGG) ‡mdtCmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (85/68);  new base C (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 = 9.99e-01

Only 100 of 172 total aligned reads displayed.

TCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCG  >  NC_002947/4070509‑4070805
                                                                                                                                                     |                                                                                                                                                   
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)
           cGACACCTGAGCCAGGCTCTGCGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACt                                                                                                                                          >  2:1442968/1‑150 (MQ=255)
            gACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTg                                                                                                                                         <  2:2436586/150‑1 (MQ=255)
              cacCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGCTCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTggg                                                                                                                                       <  2:1221992/150‑1 (MQ=255)
               acCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgc                                                                                                                                      <  2:4655363/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                                                                                                                                    >  1:530402/1‑150 (MQ=255)
                 cTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcag                                                                                                                                    >  1:1865256/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)
                           ctctGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCgg                                                                                                                          <  1:1941159/150‑1 (MQ=255)
                              tgtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTa                                                                                                                       >  2:1127030/1‑150 (MQ=255)
                              tgtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTa                                                                                                                       <  1:4757999/150‑1 (MQ=255)
                               gtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAg                                                                                                                      <  1:1328236/150‑1 (MQ=255)
                                tgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGc                                                                                                                     <  2:4523083/150‑1 (MQ=255)
                                  aaaGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGt                                                                                                                   >  1:3043645/1‑150 (MQ=255)
                                     ggATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGa                                                                                                                <  1:3124526/150‑1 (MQ=255)
                                       aTGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGCATGCGCAGCATGcccaccta                                                                                                                                       <  1:399125/124‑7 (MQ=255)
                                       aTGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGCATGCGCAGCATGcccaccta                                                                                                                                       >  2:399125/1‑118 (MQ=255)
                                             tCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCt                                                                                                        >  1:504094/1‑150 (MQ=255)
                                               ggCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTg                                                                                                      >  1:2690609/1‑150 (MQ=255)
                                                gCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGa                                                                                                     >  1:1035374/1‑150 (MQ=255)
                                                 ccAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  1:1039919/150‑1 (MQ=255)
                                                 ccAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  2:1887744/150‑1 (MQ=255)
                                                     gTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagca                                                                                                                                  <  2:3088950/116‑1 (MQ=255)
                                                     gTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagca                                                                                                                                  >  1:3088950/1‑116 (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)
                                                                    tttCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATc                                                                                 <  1:4480098/150‑1 (MQ=255)
                                                                       cTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGAt                                                                              <  2:1974747/150‑1 (MQ=255)
                                                                          cAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTc                                                                           >  1:45001/1‑150 (MQ=255)
                                                                              acaTCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACgg                                                                       <  1:353787/150‑1 (MQ=255)
                                                                                tcgtggtgTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGcc                                                                     >  2:4158034/8‑150 (MQ=255)
                                                                                tcgtggtgTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGcc                                                                     >  2:4158045/8‑150 (MQ=255)
                                                                                  ccGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCa                                                                   >  2:2019078/1‑150 (MQ=255)
                                                                                   cGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCaa                                                                  <  2:3630100/150‑1 (MQ=255)
                                                                                    gAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAAc                                                                 >  1:1948296/1‑150 (MQ=255)
                                                                                        tCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    >  2:5274477/1‑111 (MQ=255)
                                                                                        tCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGACCCAACTCGa                                                             >  2:822612/1‑150 (MQ=255)
                                                                                         cAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            <  1:77240/150‑1 (MQ=255)
                                                                                          aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGa                                                             <  2:3067246/148‑1 (MQ=255)
                                                                                          aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAg                                                           >  2:261776/1‑150 (MQ=255)
                                                                                          aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAg                                                           <  2:4412488/150‑1 (MQ=255)
                                                                                               cAGCACCCTGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  2:3258583/150‑1 (MQ=255)
                                                                                               cAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  1:945324/150‑1 (MQ=255)
                                                                                               cAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGAt                                                      <  2:2902162/150‑1 (MQ=255)
                                                                                                   aCCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACt                                                  >  1:265486/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)
                                                                                                     cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  1:482369/1‑150 (MQ=255)
                                                                                                     cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  2:520495/1‑150 (MQ=255)
                                                                                                     cATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTAc                                                >  1:420004/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                                                                                                    >  1:2494749/1‑93 (MQ=255)
                                                                                                          tCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGaa                                                                                                    <  2:2494749/93‑1 (MQ=255)
                                                                                                             gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGAc                                                   <  2:3787055/139‑1 (MQ=255)
                                                                                                             gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGAc                                                   >  1:3787055/1‑139 (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)
                                                                                                                          gTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           <  1:5245174/150‑1 (MQ=255)
                                                                                                                               tcttctTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGa                           >  2:2058740/1‑145 (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:368552/150‑1 (MQ=255)
                                                                                                                                  tcttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTgg                   >  1:4087704/1‑150 (MQ=255)
                                                                                                                                   cttGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGt                  <  2:482369/150‑1 (MQ=255)
                                                                                                                                     tGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGa                <  1:5281922/150‑1 (MQ=255)
                                                                                                                                     tGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGa                <  1:4789313/150‑1 (MQ=255)
                                                                                                                                      gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATcgccgcaac                                                                                                      <  1:166636/63‑4 (MQ=255)
                                                                                                                                      gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATcgccgcaac                                                                                                      >  2:166636/1‑60 (MQ=255)
                                                                                                                                      gTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCt                                                                                         <  1:3934418/76‑1 (MQ=255)
                                                                                                                                        ccagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            <  1:5394973/98‑1 (MQ=255)
                                                                                                                                        ccagtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                            >  2:5394973/6‑103 (MQ=255)
                                                                                                                                         cagtggCATGCGCAGCATGCCGCTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt            >  1:841608/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:530402/150‑1 (MQ=255)
                                                                                                                                              gCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCCGCTGGTGAGCGTGGTGa       >  2:4471013/1‑150 (MQ=255)
                                                                                                                                               cATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGAc      >  1:4426703/1‑150 (MQ=255)
                                                                                                                                                 tgctgcGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGc    <  1:317699/150‑1 (MQ=255)
                                                                                                                                                   ctgcGCATGCCACTGGGCAACAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCg  >  1:1564938/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGCCG  >  NC_002947/4070509‑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: