Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,329,773 A→G 16.6% S271G (AGC→GGC)  dacB → D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,329,7730AG16.6% 166.0 / 3.0 69S271G (AGC→GGC) dacBD‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base A (14/42);  new base G (11/2);  total (25/44)
Fisher's exact test for biased strand distribution p-value = 1.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCGATATT  >  NC_000913/3329536‑3330006
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
gTGGAATCTTATTCCCCCGTTACGATTTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAg                                                                                                                                                                                                                                           <  1:657092/239‑1 (MQ=255)
 tGGCATCTTATTCCCCCGTTGCGATGGTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAg                                                                                                                                                                                                                                           <  1:735349/238‑1 (MQ=255)
               cctgTTCCGGTGTTCAGCCAGGTACCCCCCCCCCCCGGTGGTTCTGCCGAAGCGCAATCCTGGGAACGGGATGGGGTGCCAGGCGACCTGACCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTgcg                                                                                                                                                                                                                            <  2:518706/236‑1 (MQ=255)
                 cGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTTCGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAGGATGAGTTAAAACAGGCGGGAGTTACCAGGAACGGGaccaacgtagcca                                                                                                                                                                                                                          >  2:39547/1‑228 (MQ=255)
                  gTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAACCAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGccc                                                                                                                                                                                                                         >  1:477659/1‑238 (MQ=255)
                  gTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATTACCTGGGGCGGGACACTGCTTCGccc                                                                                                                                                                                                                         >  1:809734/1‑238 (MQ=255)
                   ttACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTCGGGCGGGAAActtctcggccca                                                                                                                                                                                                                        >  1:744887/1‑232 (MQ=255)
                   ttACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGGGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATTCCCTTGGGCGGGACAAtgctgccgcca                                                                                                                                                                                                                        >  1:332141/1‑234 (MQ=255)
                   ttACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCAATGCCGGTGCAATTCTGAAAGATGAGTTACAACAAGCGGGTTTCACCCGGACCGGACCCCTGCTGGGcccg                                                                                                                                                                                                                        >  1:666341/1‑237 (MQ=255)
                   ttACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATAACCTGTGGCGGAGCAcagcgcggccag                                                                                                                                                                                                                        >  1:134476/1‑230 (MQ=255)
                    tACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTAGGCAGGGGCAATTCTTAAAGATGACTTAAAACAGTCGGGGAGCACACGGAGCGGAACAcagcgtcggcaaa                                                                                                                                                                                                                       >  2:843198/1‑230 (MQ=255)
                      ccaTGTCCAGCAGGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACt                                                                                                                                                                                                                    <  1:546803/237‑1 (MQ=255)
                         tGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAg                                                                                                                                                                                                                 <  1:572383/239‑1 (MQ=255)
                           ttCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGGGTTAAAACAGGCGGGTTTAACCCGGAGCGGAACACCAGTGCGcccaaaacaagt                                                                                                                                                                                                                >  2:705267/1‑229 (MQ=255)
                                 cgggTCCGCACCCTCCCCCGGGGTTCTGCCGAAGGCCAATAGTGGGAACGGGATTGGGGGCAAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTaa                                                                                                                                                                                                            <  2:195169/234‑1 (MQ=255)
                                    gtccGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGcc                                                                                                                                                                                                                          <  1:100067/217‑1 (MQ=255)
                                     tccccccccTCCCCCGGGGTTGTGCCAAAGGCAAATGCTGGGAACTGGAGGGGGGGCCAGGGGCCCTGAACCGGTTTACCGTGACGGGGTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAAcc                                                                                                                                                                                                      <  2:477659/232‑1 (MQ=255)
                                      gcGCCCCCCCCCCCGTGGGTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTaa                                                                                                                                                                                                            <  1:634584/230‑1 (MQ=255)
                                       cGCCCCCTCCCCCGTTGTTCGGCCGAATCGCAATCCTGGGAACTGGATGGGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTgg                                                                                                                                                                                                   <  2:217940/239‑1 (MQ=255)
                                             cTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATTACCTGGAGCGGAACACTCATCCGCCCAATCCCAGTTAACGAAAATGGAACgtc                                                                                                                                                                                             >  1:127558/1‑237 (MQ=255)
                                                ccccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCCCCCGGGGGGGGACACCTGTTCGCCAGAACACAGGTAACCAACCTGGaaacggtct                                                                                                                                                                                          >  2:634584/1‑232 (MQ=255)
                                                         tCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTaa                                                                                                                                                                                 <  1:802250/239‑1 (MQ=255)
                                                          cTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAAAACTGCTTCGCCCGACTCAGGTTAAAGAAACCGGGACGGTAATTTCCCGTgaa                                                                                                                                                                                >  1:478873/1‑239 (MQ=255)
                                                          cTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGTGCCAGACTCAGGTTAACGAAACTGGGACGGTAGGTGCCAGCaaa                                                                                                                                                                                >  1:743463/1‑239 (MQ=255)
                                                          cTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGAATCACCTGGGGGGGAACCCTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAAGTGCCAGGaaa                                                                                                                                                                                >  1:215407/1‑239 (MQ=255)
                                                                     cAATTCTGTGAACTGGATGTTGTGGCAGGCGACCGGAACCGATTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTGGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAAGTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGccc                                                                                                                                                                      <  1:866345/238‑1 (MQ=255)
                                                                           tGCGAACTGGATGTGGGGCGAGGGGAACTGAACCGCTTTACGCTGGCGGGATGTCTGCCAAAACGTTCTGAGCCGCGCCCGTTGGCTTTTGCCGTGCAGGATGGGGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAACGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCa                                                                                                                                                               <  1:760576/239‑1 (MQ=255)
                                                                             cGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGGCCCCCCGCACg                                                                                                                                                             >  1:292490/1‑239 (MQ=255)
                                                                             cGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTGAACAGTCGGCCCCGCTGCACg                                                                                                                                                             >  1:451335/1‑239 (MQ=255)
                                                                             cGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCCGACCCAGGTTAACGAACCTGGAAACGTAGTTGCCAGGAAACAATCGGgccccccgcccg                                                                                                                                                             >  1:290397/1‑231 (MQ=255)
                                                                             cGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGAGGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGACCACTGCTGCGCCTGACACAGGGGAAAGAAACCGGGACGGTGGTGGGCGGTAAAAAACCGGccccccggccca                                                                                                                                                             >  2:149295/1‑231 (MQ=255)
                                                                               accgtgttttggtgCTAGGCCCCCTGGACCTTTTCACGCTGATGGGTTTGCTGCCACAACGTTCCGGACCGCCCCCGTTGGCTTTTGCACTGCAGGGTGGGGCCAGCTATGCCGGTGCAATTATGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                                                                           <  2:176663/230‑1 (MQ=255)
                                                                               aatggggtttggggCCGGGCGCCCTGACCCCGTTTCCGCTGGCGGGGTCCCTGGCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                                                                           <  2:744887/230‑1 (MQ=255)
                                                                                acgggagtgggtgCCAGGCGACCGGAACCGCTTTAGCCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                                                                           <  1:227317/229‑1 (MQ=255)
                                                                                 cTGGATTGGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCg                                                                                                                                                                    <  1:246685/228‑1 (MQ=255)
                                                                                  tGGGTGTGGTGCCGGGCGACCGGAACCGCTTTACGTGGACGGGATGCCTGCCACAACGTTCTGGGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                                                                           <  2:645689/236‑1 (MQ=255)
                                                                                           tgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGGACACTGCTGCGCCAGAACCAGGGTAACGAACCTGGGAAGGGAGGTGCCAGTAAAAAAACGGCCCCCCTGCACGATCtctggctctta                                                                                                                                               >  2:246685/1‑229 (MQ=255)
                                                                                              cAGGCGGCCTGAACCGCTTTACGGTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATg                                                                                                                                            <  1:576371/239‑1 (MQ=255)
                                                                                               tggCGACCTGAACCGCGTTACGCTGAAGGAATGCCTGCCACAACGTTGTGAGCCGTTCCCGTTGGCTTTTGCCGGGCAGGATGGAGCCATCTATGCCAGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCTGAACACTGCTGCGCCAGGCTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATg                                                                                                                                            <  1:216844/237‑1 (MQ=255)
                                                                                                    cccTGACCCGCTTTTGGTGGGCGGGGTGCCTGCCAAAACGTTCTGGGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGaa                                                                                                                                       <  2:292490/237‑1 (MQ=255)
                                                                                                    aCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAAAACTGCTGCGCCCGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGGCCCCCTGCCCGGACGGCTTAAGATTATGGTGgaa                                                                                                                                      >  1:20548/1‑239 (MQ=255)
                                                                                                             gCTTTACGCGGAGGGGGTTCCTGCCACAACGTTCTGAGCCGCGCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGac                                                                                                                             <  1:737392/239‑1 (MQ=255)
                                                                                                                     tgggcGGGGTGCCTGCCAAAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGGTCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGAt                                                                                                                     <  2:525679/235‑1 (MQ=255)
                                                                                                                        aCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGGACGGTAGGTGCCCGTAAACAGTCGGGCCCGCTGCACGATCTGCTTAAGATTATTCTGAAAAAGCCGGGCAACATGAtag                                                                                                                   >  1:215199/1‑236 (MQ=255)
                                                                                                                                      acaaggtgttGAGCCCTTCCCGTTGGCTTTTGCCGGGCAGGATGGGGCCAGCTGTGCCGGGGAAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACGCTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTcc                                                                                                    <  2:134476/230‑1 (MQ=255)
                                                                                                                                      acaACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTcc                                                                                                    <  1:682311/239‑1 (MQ=255)
                                                                                                                                       caACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCg                                                                                                   <  1:705267/239‑1 (MQ=255)
                                                                                                                                       aaaCGTTCTGGGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTGAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCg                                                                                                   <  1:843198/238‑1 (MQ=255)
                                                                                                                                        aaCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGc                                                                                                  <  1:607937/239‑1 (MQ=255)
                                                                                                                                            ttCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                                                                           <  1:670503/178‑1 (MQ=255)
                                                                                                                                              cTGAGCCGCTCCGGTTGGCTTTTGACGGGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGAt                                                                                             <  1:803200/238‑1 (MQ=255)
                                                                                                                                                      gtCCCGTTGGGTTTTGCCGTGGAGGATGGGGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACCCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATgc                                                                                    <  2:215407/238‑1 (MQ=255)
                                                                                                                                                       ccccGTTGGCTTTTGCCGTGCAGGATGGCGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTGAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACACCATGATCGCCGATACGGTTTTCCGCATGATAGGCCATgcg                                                                                   <  1:686346/238‑1 (MQ=255)
                                                                                                                                                           gTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGAGGATTACCTGGGGCGGGACACTGCTGCGCCAGAATCAGGGTAACGAACCTGGAACCGTAGGTGTCCGTAAAAAATCGGCCCCGCTGCACGATCTGCTTTAGATTATGCTGAAAAAATCGGGCAACAAGATGGCCGATAAGGGTTTCCGCATAAAAAGCCaaggccc                                                                                 >  2:323015/1‑231 (MQ=255)
                                                                                                                                                              ggtttttGCCGTTCAGGATGGGGCCGGTTATGCGTGTGCAATTCTGAAAGATGGTTTAAAAATGGCGGGTATCCCCTTGAGCTGAACGCTGCTGCCCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCa                                                                            <  2:298815/236‑1 (MQ=255)
                                                                                                                                                               gCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                                                                                                                                      <  1:811061/164‑1 (MQ=255)
                                                                                                                                                                 ttttGCCGTGCAGGATGGGGCCAGGTATGCCGGTGCAATTCTGAAAGGTGAGTAAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAt                                                                          <  2:215199/238‑1 (MQ=255)
                                                                                                                                                                 ttttGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAt                                                                          <  1:701533/238‑1 (MQ=255)
                                                                                                                                                                 ttttGCCGGGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGGGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAt                                                                          <  1:864099/238‑1 (MQ=255)
                                                                                                                                                                 ttgtGCCGTGGCGGATGGGTCCAGCGTTGCGGGTTCAGATCTTGAAGAGGAGTAAAAACAGGGTGGTAACCCCTGCGGCGGAACACTGATTCGACGGACTCATGTGAACGAACCTGGAACGGTAGTTGCCATGAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAt                                                                          <  2:893/235‑1 (MQ=255)
                                                                                                                                                                                    gCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACCCTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCggg                                                       <  1:540375/238‑1 (MQ=255)
                                                                                                                                                                                         cTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                                                                           <  1:765787/133‑1 (MQ=255)
                                                                                                                                                                                             gCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGa                                             <  1:478016/239‑1 (MQ=255)
                                                                                                                                                                                                   gaaaTTCCGAAAGAAGAAATAAAACGGGCGGGTCTTACTTGGGGCGGGACACTGGTGCGCCAGACTTGGGTTAAGTAACCTTGAATGGTAGTTTCCAGTAAACAGTCGGCCCCGCTGCCCGATCGTCTGAAGATTAGGCTGAAAAAGTCGGACAACATGATCGCCGATTCGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCg                                        <  2:666341/236‑1 (MQ=255)
                                                                                                                                                                                                   gCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGGCAGACTCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCCCTGCACCATCTGCTTaaaa                                                                                                                                                 >  2:765787/1‑131 (MQ=255)
                                                                                                                                                                                                   gCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGGGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGGCGGGCGcccgg                                      >  1:745037/1‑236 (MQ=255)
                                                                                                                                                                                                                               cGGGTATCACCTGGAGCGGAACCCTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGc               <  1:489240/235‑1 (MQ=255)
                                                                                                                                                                                                                                  gTATCCCGTGGAGGGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGgtgt         <  1:2003/238‑1 (MQ=255)
                                                                                                                                                                                                                                        ccTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAAAATGGCGGGCCGGGGCGGAGGCCGTGCGTCCGAAACCGCGCCAGCAAACCGGGGTGGagaag  >  2:295919/1‑235 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
GTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCGATATT  >  NC_000913/3329536‑3330006

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