Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,081,393 T→G 12.4% intergenic (‑231/+153) sdaC ← / ← amaC serine:H+ symport permease/D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,081,3930TG12.4% 619.8 / 9.4 229intergenic (‑231/+153)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (80/105);  new base G (29/7);  total (116/113)
Fisher's exact test for biased strand distribution p-value = 4.15e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

GGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGC  >  NC_002947/4081245‑4081542
                                                                                                                                                    |                                                                                                                                                      
tcgacagcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCtt                                                                                                                                                       <  1:2572316‑M1/15‑1 (MQ=255)
tcgacagcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCtt                                                                                                                                                       <  2:8080008‑M1/15‑1 (MQ=255)
gggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTT                                                                                                                                                       <  2:2552182‑M1/150‑1 (MQ=255)
gggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTT                                                                                                                                                       <  2:7170079‑M1/150‑1 (MQ=255)
gggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTT                                                                                                                                                       <  1:4898112‑M1/150‑1 (MQ=255)
gggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTT                                                                                                                                                       <  1:326265‑M1/150‑1 (MQ=255)
  gcgcGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGC                                                                                                                                                     <  1:1331135‑M1/150‑1 (MQ=255)
  gacagcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGc                                                                                                                                                     >  1:3606431‑M1/134‑150 (MQ=255)
    gcgATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGT                                                                                                                                                   <  1:1146598‑M1/150‑1 (MQ=255)
     gccggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaccgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                  >  2:820814‑M1/131‑150 (MQ=255)
     agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                  >  1:7441407‑M1/131‑150 (MQ=255)
     agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                  >  1:5855717‑M1/131‑150 (MQ=255)
     agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                  >  1:7285734‑M1/131‑150 (MQ=255)
       cggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                                <  1:9022233‑M1/22‑1 (MQ=255)
       cggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                                >  1:3774287‑M1/129‑150 (MQ=255)
       aTTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGC                                                                                                                                                <  2:5616404‑M1/150‑1 (MQ=255)
              aCGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGCTAggg                                                                                                                                         >  2:2348777/1‑150 (MQ=255)
               cGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                        <  1:8242868/150‑1 (MQ=255)
                gtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGcc                                                                                                                                       <  1:8378191‑M1/31‑1 (MQ=255)
                ggAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGCTAGGGcc                                                                                                                                       >  2:7308671/1‑150 (MQ=255)
                          cACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                             <  2:3175173/150‑1 (MQ=255)
                           aCTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                            >  1:4015743/1‑150 (MQ=255)
                            cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGGTGGTTTCTATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGccc                                                                                                                           <  1:6382789/150‑1 (MQ=255)
                            cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGTGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGccc                                                                                                                           <  1:1424003/150‑1 (MQ=255)
                             ttGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGGGTGGCGGGAGGGCCGCAAAGCGGcccc                                                                                                                          >  1:7837308/1‑150 (MQ=255)
                                cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCCGc                                                                                                                       >  2:6110876/1‑150 (MQ=255)
                                cGAAGCCCTGTTCCGAGGAGCGGTTGCACTTGTCTTTACCCGCGGAAGGCATTGCTGGTGTCGCTCGGTGGAATGGAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCTTTGCGGccctatc                                                                                                                       <  2:7520074/150‑5 (MQ=255)
                                 gAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCatcg                                                                                                                      <  1:211079/150‑4 (MQ=255)
                                  aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGGGGCGCGGGTGGGGCGCCAAGCGGCCCCAGcca                                                                                                                     >  1:8080262/1‑148 (MQ=255)
                                   aGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAt                                                                                                                    <  1:3814357/150‑1 (MQ=255)
                                    gCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAtt                                                                                                                   <  1:8644188/150‑1 (MQ=255)
                                     cccTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAttt                                                                                                                  <  1:8015474/150‑1 (MQ=255)
                                         tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                              <  2:7060173/150‑1 (MQ=255)
                                             cGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCGCAAGTTCAGCGGCa                                                                                                          >  2:7391999/1‑150 (MQ=255)
                                             cGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCa                                                                                                          >  1:4187216/1‑150 (MQ=255)
                                              gaggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCaa                                                                                                         >  1:6981646/1‑150 (MQ=255)
                                                      gTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCTTGCGTCCCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCAGCaa                                                                                                >  2:2799462/1‑150 (MQ=255)
                                                      gTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCaa                                                                                                >  1:8890343/1‑150 (MQ=255)
                                                      gTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGCTGGGGCCGCAAAGCGGCCCCCGCACAGTCAGCGGCcaaccgggca                                                                                                 >  1:630052/1‑143 (MQ=255)
                                                          cACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGGCGCCCTAGGGCCGCCAAGCGGCCCCCGCAAATTCAGCGGCAAAGCGGCCATCgg                                                                                            >  2:5479634/1‑149 (MQ=255)
                                                           aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  2:7607988/150‑1 (MQ=255)
                                                           aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  1:4152675/150‑1 (MQ=255)
                                                             ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCGTCGCGATAGGGCCCCAAAGCGGCCCCACCAATTACACCCGCAAAGCGGCAAACGTAgg                                                                                         >  2:583078/1‑150 (MQ=255)
                                                             ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                         <  1:450392/150‑1 (MQ=255)
                                                                  tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCGTCGCGAt                                                                                                                                             >  2:4243220/1‑94 (MQ=255)
                                                                  tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATATCAGCGGCAAAGCGGCAATCGTAGGGgcgc                                                                                    >  2:3599535/1‑150 (MQ=255)
                                                                  tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAATTCAGCGGCAAAGCTGCAAACGTAGGGgcgc                                                                                    >  2:8972103/1‑150 (MQ=255)
                                                                  tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGTTAGGGCCGCCAAGCCGCCCCCGCAAATTCACCGGCAAAGCTGCCAACCTAGGGgcgc                                                                                    >  2:4768256/1‑150 (MQ=255)
                                                                  tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCCAATTCAGCCGCAAAGCAGCAAACGTATGGgggc                                                                                    >  2:1921183/1‑150 (MQ=38)
                                                                  tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGGTTGGGCCGCAAAGCGGCCCCCCCAAAATCACCGGCAAACCTGCAAACGTAGGGgggc                                                                                    >  2:3320540/1‑150 (MQ=255)
                                                                     aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGGGGCGGGGTGGGGGCGCAAAGCGGCCCCCGCCAGTTAAGCGGCCAAGCTGCAAGACTGGGGGGGCTcc                                                                                 >  1:5223727/1‑150 (MQ=255)
                                                                       aaGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTc                                                                                                                 >  1:3639703/1‑117 (MQ=255)
                                                                            gAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGGGGGGGGCGGGAGGGCCGCCAAGCGGCCCCCCCAAGTTAAGCGGCCAAGCAGGACAAGGGGGGGCGCACCGCAcccg                                                                          >  1:5285840/1‑149 (MQ=255)
                                                                            gAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTCGTCGCGATAGGGCCGCAAAACGGCCCCACCAAATTAAGCGGCAAAACTGCAATAGTAGGGGCGCGCCGCAccca                                                                          >  2:3442109/1‑149 (MQ=38)
                                                                               ggCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCAAGCCACGCGAGATGGCCGCAAAGCGGCCCCCCCCATTAAAGCCGCCAAGCAGCAAGAGTATGGGCGCAccccccccaatt                                                                       >  2:1180318/1‑146 (MQ=255)
                                                                                gCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtc                                                                      >  2:8890186‑M2/1‑148 (MQ=255)
                                                                                 cATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCATGCCTCGCGAAAGGGCCGCAAAACGGCCCCAGCAATTTCAACCGCAAACCTTCAAACGTAGGGGCGCAACGCACCCATTtag                                                                     <  2:5460508‑M2/150‑4 (MQ=255)
                                                                                  aTTTCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgc                                                                    >  2:6626189‑M2/1‑146 (MQ=255)
                                                                                  aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgc                                                                     <  2:6082919‑M1/150‑74 (MQ=255)
                                                                                  aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCa                                                                                                                <  1:4389114/107‑1 (MQ=255)
                                                                                  aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCa                                                                                                                >  2:4389114/1‑107 (MQ=255)
                                                                                    tGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcga                                                                  >  2:2640343‑M2/1‑144 (MQ=255)
                                                                                     gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgac                                                                  <  2:8216180‑M1/150‑77 (MQ=255)
                                                                                     gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTCtcgcgac                                                                 >  2:2253155‑M2/1‑143 (MQ=255)
                                                                                     actcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                 >  2:8378191‑M1/51‑150 (MQ=255)
                                                                                        ggTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCAGCAAATTGAGCGGCAAAGCTGGAATCGGAGGGGCGCTCCCCACCCATTtggcggcgcc                                                              <  2:2274167‑M2/150‑11 (MQ=255)
                                                                                         gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaa                                                             >  1:2758054‑M2/1‑139 (MQ=255)
                                                                                         gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaa                                                             >  1:3328000‑M2/1‑139 (MQ=255)
                                                                                         gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaa                                                             >  1:3394206‑M2/1‑139 (MQ=255)
                                                                                          ccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAg                                                            >  1:469682‑M1/46‑150 (MQ=255)
                                                                                               aTCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtc                                                                      <  1:2641899‑M2/135‑3 (MQ=255)
                                                                                                 cGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTGGTCGCCATAAGGGCCCAAAACCGCCCCACCAAATTAAGCCGCAAAGCCGCCAACCGGAGGGCGCGCCCCACCCCGTTCGCGACCCAAcgccgcga                                                     >  2:77063/1‑148 (MQ=255)
                                                                                                   gTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCCCCAAATTCAGCGGCCAAGCTGCAATCGTTGGGGCGCTACGCCCCCATTtcgcgacgcaaggccgctccc                                                   <  2:5074453‑M2/150‑22 (MQ=255)
                                                                                                        aTGTTGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACaaa                                              <  2:5323724/150‑1 (MQ=255)
                                                                                                        aTGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACaaa                                              <  2:6809258/150‑1 (MQ=255)
                                                                                                                   tGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCCTCGCGATAGGGCCCCAAAGCCGCCCCACCAAAATAACCCGCAAAGCAGCAAAAGTATGGGCCCTCCCCAACCATTTAGCGAACCAATGCCCCACCCCCAAACACCaaacccc                                   >  2:8461627/1‑145 (MQ=255)
                                                                                                                         ccTGGGGGGCGGCGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCAGAAATTACCGCGGAAACGCTGAAATGGTCGGGGCGCTCCGCACCCTTTTCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCcac                             <  2:1631713/150‑1 (MQ=25)
                                                                                                                               gggCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCATGGCCGCTCCCACAAAGTCCAACCGCGCCCCACaggagg                       <  2:4015743/150‑1 (MQ=255)
                                                                                                                               gggCGGCCTGTGGGAGCGGCCATGCGTCGCGATAGGGCCGCAAACCGGCCCCACCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCGCCGCACCCGTTTCGCGTCGCAAGGCCGCGCCCACCAAGGCCAACCCCGCCCCACacgagg                       >  2:1343008/1‑150 (MQ=25)
                                                                                                                                 gcggcCGGTGGGCGCGGCGTTGCGTCGCGATCGGGCCGCAACGCGGCCCCAGCAATGTCAGCGGCAACGCTGCAATCGTCGGGGCGCTCCGCACCCTTTGCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAgggg                     <  2:8669461/150‑1 (MQ=255)
                                                                                                                                  cggcCGGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGAAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTCCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGc                    <  2:6263639/150‑1 (MQ=255)
                                                                                                                                         tgGGAGCGGCCTTGCGTCGCTAaatgggggcggcgcgcccctccgtttgcagctttgccgcttaaattggtggggccgttttggggccctatcgcgacgcacggccgctcccacaggccgccccacaggcgcccaggccgctacattcca              <  2:4749269‑M1/150‑129 (MQ=255)
                                                                                                                                         TGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAACGCTGCAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGtt             >  2:6466626‑M1/1‑150 (MQ=255)
                                                                                                                                          gggCGCGGCGTTGCGTGGCGAactgggtgcggagcgcccctacgtttgccgctttgcggctgaatttgcgggggccgctttgcggccatatcgcgccgcacggccgctcccacaggccgccccacaggcgcccaggccgctacattccac             <  2:4740136‑M1/150‑130 (MQ=255)
                                                                                                                                          gggCGCGGCCTTGCGTCGCTAaatgggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctaccgcgaggcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccac             <  2:8797806‑M1/150‑130 (MQ=255)
                                                                                                                                          gggCGCGGCCTTGCGTCGCGAaatgggtgcggagcgcccctacgattgcagctttgccgcttaaattgctggggccgctttgcggccctaccgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccac             <  2:2186871‑M1/150‑130 (MQ=255)
                                                                                                                                            GAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCACTCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACCGGAGGGGCCGCAGTTGGc          <  2:5758591‑M1/150‑1 (MQ=255)
                                                                                                                                             aGCGGCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccga          >  1:6032526‑M1/1‑18 (MQ=255)
                                                                                                                                              gCGGCCTTGCGTCGCTAaacgggggcggagcgcccctacttttgcaggtttgccgcttaatttgctggggccgctttggggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgat         <  2:4462068‑M1/150‑134 (MQ=255)
                                                                                                                                               cGGCCTTGCGTCGCGAaacgggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatc        <  2:8135124‑M1/150‑135 (MQ=255)
                                                                                                                                                ggCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcg       >  1:2909127‑M1/1‑15 (MQ=255)
                                                                                                                                                 gCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcga      >  2:418931‑M1/1‑14 (MQ=255)
                                                                                                                                                 gCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcga      >  2:5336156‑M1/1‑14 (MQ=255)
                                                                                                                                                  ccTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgac     >  1:2740188‑M1/1‑13 (MQ=255)
                                                                                                                                                  ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgac    <  1:88840‑M2/150‑69 (MQ=255)
                                                                                                                                                   CTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTgg   <  1:5484406‑M1/150‑1 (MQ=255)
                                                                                                                                                    ttGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacac   >  1:2883869‑M1/1‑11 (MQ=255)
                                                                                                                                                    TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAACGCTGCAATCGTAGGGGCGCTCCACACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTGGGCCCGTTGGc  <  1:6115994‑M1/150‑1 (MQ=255)
                                                                                                                                                    TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc  <  1:2063255‑M1/150‑1 (MQ=255)
                                                                                                                                                    TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc  <  1:3865304‑M1/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGC  >  NC_002947/4081245‑4081542

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: