Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,034,775 T→C 9.9% intergenic (+61/+75) dppA‑IV → / ← PP_5284 dipeptide transport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,034,7750TC 9.9% 316.0 / 18.6 115intergenic (+61/+75)dppA‑IV/PP_5284dipeptide transport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/61);  new base C (6/5);  total (48/67)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 117 total aligned reads displayed.

GGCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAG  >  NC_002947/6034627‑6034918
                                                                                                                                                      |                                                                                                                                               
ggcgACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAATCCCGCTCCTACAGGGGTAGCGCCAGGCGGGGAGtt                                                                                                                                                >  2:5141745/1‑150 (MQ=255)
ggcgACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGtt                                                                                                                                                >  1:489440/1‑150 (MQ=255)
   gACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGAc                                                                                                                                             <  2:919022/150‑1 (MQ=255)
   gACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGAc                                                                                                                                             >  2:382216/1‑150 (MQ=255)
   gACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGAc                                                                                                                                             <  2:2658325/150‑1 (MQ=255)
       ttGCTGCGCCAAGGCGTCGATGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                         <  2:1172169/150‑1 (MQ=255)
       ttGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                         <  1:1499214/150‑1 (MQ=255)
       ttGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                         <  1:5170446/150‑1 (MQ=255)
        tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                        <  1:4504380/150‑1 (MQ=255)
        tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                        <  1:4082315/150‑1 (MQ=255)
         gctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCACCAACAGGGGAAGAGCCAGCCTGGCAGTTTAAGCTAcc                                                                                                                                       >  1:2209814/1‑150 (MQ=39)
          ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                      <  1:1514963/150‑1 (MQ=255)
          ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                      <  1:1514966/150‑1 (MQ=255)
          ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                      <  2:3475941/150‑1 (MQ=255)
          ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                      <  1:339707/150‑1 (MQ=255)
          ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                      <  1:339411/150‑1 (MQ=255)
          ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                      <  1:339405/150‑1 (MQ=255)
          ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                      <  2:2473966/150‑1 (MQ=255)
            gcgcCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCt                                                                                                                                    <  2:5330438/150‑1 (MQ=255)
              gcCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGt                                                                                                                                  <  1:468119/150‑1 (MQ=255)
               ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTa                                                                                                                                 <  2:3231708/150‑1 (MQ=25)
               ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGACTGGCGCTACCCCTGTa                                                                                                                                 <  2:5279961/150‑1 (MQ=25)
               ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTa                                                                                                                                 <  2:489440/150‑1 (MQ=255)
               ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTa                                                                                                                                 <  1:3534444/150‑1 (MQ=255)
               ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTa                                                                                                                                 <  1:1992552/150‑1 (MQ=255)
               ccAAGGCGTCGAGGGGGACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCGTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTa                                                                                                                                 <  2:3854067/150‑1 (MQ=25)
                cAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAg                                                                                                                                <  1:2034097/150‑1 (MQ=21)
                   ggCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCGTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAg                                                                                                                             <  2:3958820/150‑1 (MQ=255)
                      gTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCgg                                                                                                                          <  2:2405907/150‑1 (MQ=255)
                      gTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCgg                                                                                                                          <  1:3318141/150‑1 (MQ=255)
                       tCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCggg                                                                                                                         <  1:4384904/150‑1 (MQ=255)
                       tCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCggg                                                                                                                         <  1:4384868/150‑1 (MQ=255)
                       tCGAGGGGGCCCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTTAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCGGTAGGAGCggg                                                                                                                         <  2:3330136/150‑1 (MQ=255)
                       tCGAGGGGGACCAGGTGAGCCCGTTCGGGCGGCGGGATTTCAGCCCGGTGACGGTGACCAAGTACGCCTGTACCGGCCGCTGCGCGGGTAACCCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGGAGGAGCggg                                                                                                                         <  2:3330141/150‑1 (MQ=255)
                                aCCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCcgc                                                                                                                <  2:4786491/150‑1 (MQ=255)
                                 ccAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCGCTTCGCGGGTAACCCCGCTCCTACCGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCcgcg                                                                                                               <  2:4807240/150‑1 (MQ=255)
                                    ggTGAGCCCGTTCGGGCGGCTGTATTTCCGCATGGTGACGGTGACCACGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCTGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                            <  1:774942/150‑1 (MQ=38)
                                      tGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaagag                                                                                                          <  2:2623222/150‑4 (MQ=255)
                                        cgCCGGTTGGGGCGGCTGGTTTTCAGCACGGTGACGGTGACCACGTAAGCCTGTACCGGCCTCTTCGCGGGTAACCCCGCTCCTACCGGGGTAGCGCCAGCCGGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                        <  2:1219093/149‑1 (MQ=37)
                                            cGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGACCTCTTCGCGGGTTAAACCCCGCCAACAGGGGTAGCGCCCCCCGGGGCGGTGTCCCGACCCCCGGAGGGGCGGGGGTTCCCCCGCAACGGccccg                                                                                                    >  1:4159112/1‑148 (MQ=255)
                                                cGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTTAAAACCCGCCGAAAGGGGGAGCGCCCCCCGGGCCGGTGTCGCGACCCCGGGGTGAGGGGGGTTAACCCCGCAAGGGCCCGGAcag                                                                                                >  1:4133241/1‑148 (MQ=255)
                                                       cTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCCGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGATGACGCAACCCCAGTAGGAGCCGGGTTACCCGCCACAGGGCCCCTACTGAcccaaa                                                                                         >  2:3106981/1‑147 (MQ=37)
                                                        tGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGa                                                                                                                                              >  2:4881579/1‑96 (MQ=255)
                                                        tGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGa                                                                                                                                              <  1:4881579/96‑1 (MQ=255)
                                                        tGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGt                                                                                                                                  >  1:4447045/1‑108 (MQ=255)
                                                        tGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGt                                                                                                                                  <  2:4447045/108‑1 (MQ=255)
                                                         ggATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATaa                                                                                       <  2:3229658/150‑1 (MQ=255)
                                                               cAGCAAGGTGCCGGTCACAAAGAACGCCTGCACCGGCCTCTTCGCGGTAAAACCCGCTCCTACCGGGGTCGCGCCAGCCGGGCAGTTGCCGCAACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGctct                                                                                 <  1:3857646/150‑1 (MQ=37)
                                                               cAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCGACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCCTACTGACCCGATAAGGctcg                                                                                 >  2:2085323/1‑149 (MQ=255)
                                                                 gctaGGTGACGGTGACCATGTAAGCATGTACGGGCCTCTTCGCGGGTTATCCCGCTCCTCCAGGGGTAGCGCCAGCCTGGCATTTGAGGCTACCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCa                                                                               <  2:2209814/147‑1 (MQ=255)
                                                                  cAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAg                                                                              <  2:3626488/150‑1 (MQ=255)
                                                                     ggtgacggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACg                                                                                                                                            <  1:4534137/85‑1 (MQ=255)
                                                                     ggtgacggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACg                                                                                                                                            >  2:4534137/1‑85 (MQ=255)
                                                                           ggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCGCCGACAGGGGTCGCGCCAGCCGGGGAGGTGACGCGACCCCGGGCGGTGCGGGTTTTGCCGCGACACGGCCCGTGATGGCCCGATATGGGTCTGCGTGCTcccgc                                                                     >  2:1842838/1‑148 (MQ=255)
                                                                            gtgacCATGTAAGCCTTTACCGGCCTGTTCGCGGGTAAACCCGCTCCTACAGGGGTTGCGCCAGCCTGGCTGTTTACGCTACCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCa                                                                    <  2:1909373/150‑1 (MQ=38)
                                                                              gacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCAGGAGGGGCGGGTTTACCCGCGAAAGGGCCCGTACAGACCCTATAAGGCGCTCAGTACTCCCTCAGa                                                                  >  2:4137735/1‑150 (MQ=38)
                                                                                     tAAGCCTGTACCGGCCTCTTCGCGGGTAAACCGGCTCCTACCGGGGTAGCGCCAGCCGGGCCTTTGACGCTACCCGTGTTGGCGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAAt                                                           <  2:1381477/150‑1 (MQ=38)
                                                                                      acgCCTGTCCGGGCCTCTTGGCGGTAAAACCCGCTCCTACCGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGCGGGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATc                                                          <  1:3725928/148‑1 (MQ=255)
                                                                                          cTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCCGGTGACGCTACCCCGGTAGGAGCGGGGTTACCCGCCCAAGGGCCCCTACTGACCCAATAAGGCTCGCAGTACTCCCGCAGACCACACGCCAcc                                                      >  2:1880697/1‑150 (MQ=25)
                                                                                           tGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGTTGGGCCAGCGTGGCAGTTGACGCTACCCGTGTTGGAGGGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAccc                                                     <  2:1923213/150‑1 (MQ=38)
                                                                                            gTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCa                                                    <  2:2565542/150‑1 (MQ=255)
                                                                                                cGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCt                                                >  2:5337068/1‑150 (MQ=255)
                                                                                                   ccTCTTCGCGGGTAAACCCGCTCCTACAGGGGTTGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGa                                             <  2:847594/150‑1 (MQ=255)
                                                                                                        gcgcgGGTAAACCCGCTCCTCCAGGGGAAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATc                                        <  2:4380780/149‑1 (MQ=255)
                                                                                                          gcgGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGa                                      >  1:2423361/1‑150 (MQ=255)
                                                                                                               tAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACg                                 >  2:4240060/1‑150 (MQ=21)
                                                                                                                 aaCCCGCTCCTACAGGGTTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTTGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGc                               <  1:4918448/150‑1 (MQ=255)
                                                                                                                 aaCCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGc                               >  1:4407430/1‑150 (MQ=255)
                                                                                                                  aCCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCt                              >  1:3036211/1‑150 (MQ=21)
                                                                                                                   cccGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTc                             >  2:4408820/1‑150 (MQ=21)
                                                                                                                    ccGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTAc                                                                                                  <  1:1740273/80‑1 (MQ=255)
                                                                                                                    ccGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTAc                                                                                                  >  2:1740273/1‑80 (MQ=255)
                                                                                                                    ccGCTCCCCCAGGGGTCGCGCCAGCCTGCCAGATGACGCTACCCCCGACGGAGCGGGATTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcc                            <  2:4946348/150‑1 (MQ=255)
                                                                                                                     cGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTccc                           <  1:382216/150‑1 (MQ=255)
                                                                                                                      gCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcccc                          >  1:3693775/1‑150 (MQ=255)
                                                                                                                         ccTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATc                       >  2:2130283/1‑150 (MQ=25)
                                                                                                                          cTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCa                      >  2:691560/1‑150 (MQ=255)
                                                                                                                           tACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                     >  1:672595/1‑150 (MQ=25)
                                                                                                                           tACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCGGTACTCCCTCAGACCCCCATCCACCCATGCTCGACCATCGACAACGGCTCCCCCTCAc                     >  2:624532/1‑150 (MQ=255)
                                                                                                                           tACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                     >  2:624649/1‑150 (MQ=255)
                                                                                                                           tACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                     >  2:2138673/1‑150 (MQ=255)
                                                                                                                           tACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                     >  2:3286122/1‑150 (MQ=255)
                                                                                                                            aCAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAcc                    >  1:4589418/1‑150 (MQ=255)
                                                                                                                             cAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGTCCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCa                   >  2:5376947/1‑150 (MQ=25)
                                                                                                                              aGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCaa                  >  1:3261185/1‑150 (MQ=255)
                                                                                                                              aGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCaa                  >  1:598505/1‑150 (MQ=255)
                                                                                                                                gggTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatg                >  1:1770935/1‑150 (MQ=255)
                                                                                                                                  gTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgat              >  1:2152273/1‑150 (MQ=255)
                                                                                                                                  gTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgat              >  1:5227996/1‑150 (MQ=255)
                                                                                                                                   tAGCGCCACCCTGGCAGTTGACGCTACCCCTGGAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACAGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAaagagg             >  2:2907363/1‑148 (MQ=255)
                                                                                                                                    aGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACACTCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatga            <  2:4407430/150‑1 (MQ=255)
                                                                                                                                     gcgcCAGCCTGGCCGTGCCCCCAACCCCTTTTGGCGCGGGTTAACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat           <  1:4593697/150‑1 (MQ=255)
                                                                                                                                     gcgcCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat           <  1:3448677/150‑1 (MQ=255)
                                                                                                                                        ccAGCCTGGCAGTTGACGCTACCCCTGTAGGAGGGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGAtgtg        <  1:3783312/150‑1 (MQ=255)
                                                                                                                                             cTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCa                                                                               >  2:4489975/1‑74 (MQ=255)
                                                                                                                                             cTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCa                                                                               <  1:4489975/74‑1 (MQ=255)
                                                                                                                                             cTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCa   >  2:203168/1‑150 (MQ=255)
                                                                                                                                             cTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGACGTGATCCa   >  2:203157/1‑150 (MQ=255)
                                                                                                                                             cTGGCAGTTGACGCTACCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCa   >  1:1083850/1‑150 (MQ=255)
                                                                                                                                              tGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGAGGTGATCCAg  >  2:5014150/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GGCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAG  >  NC_002947/6034627‑6034918

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