Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,115,770 C→T 6.3% A485T (GCA→ACA)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,115,7700CT 6.3% 187.7 / 4.2 66A485T (GCA→ACA) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (38/23);  new base T (2/2);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

CGCCGCGTCCTAACCATTTGTAGCCAGCCAGCGCCGCTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGC  >  NC_000913/4115498‑4116047
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                     
cgcctcGTCCTAACCCGTTGTAGACGGCCAGCGCCGCTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGt                                                                                                                                                                                                                                                                                     <  4:121692/275‑1 (MQ=255)
cgccgcGTCCTAACCATTTGTAGCCAGCCAGCGCCGCTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGt                                                                                                                                                                                                                                                                                     <  2:336096/275‑1 (MQ=255)
   cgcgTCCTAACCATTTGTAGCCAGCCAGCGCCGCTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTGCCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAAc                                                                                                                                                                                                                                                                                  >  2:431664/1‑275 (MQ=255)
    tgggcctaccattttGTGGCCGGCCAGGGACGGTGATTCGGGGAAGGGGGGAAATTCGAGGGCAAGGGCACGTTTCATAGCAAACTGTAAGCAGTAAGGAAGGGCGCGAAGTGTAGCAAATGGGAAGGGAGGGATTATGGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCAGTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAAc                                                                                                                                                                                                                                                                                  <  4:456348/262‑1 (MQ=255)
                tttGTAGCCAGCCAGCGCCGCTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGCGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCCTGTTCTTCCCACGCCATCCCTCGTTTAACCGCTTTTTTCCACCCTTCTTACCGGTAAtttccct                                                                                                                                                                                                                                                                      >  4:74935/1‑269 (MQ=255)
                  tGTAGCCAGCCAGCGCCGCTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAg                                                                                                                                                                                                                                                                   >  2:373251/1‑275 (MQ=255)
                              gagccgcGGTTTCGGGGACGCGGGAAAATTCGATGGCAAGTTCTCGTATAATAGGAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGc                                                                                                                                                                                                                                                       <  4:373250/273‑1 (MQ=255)
                                   gcTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGGCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGCGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGa                                                                                                                                                                                                                                                  <  2:74935/275‑1 (MQ=255)
                                                 gcgGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGCCGGAACTCGCGCt                                                                                                                                                                                                                                     >  2:356863/1‑274 (MQ=255)
                                                 gcgGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAAt                                                                                                                                                                                                                                 >  1:90949/1‑278 (MQ=255)
                                                        aaTTGGATGGCAAGTTCGCGGTTCATAGCAAACTCTAAGCAGTAAGGAAGGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCAc                                                                                                                                                                                                                              <  2:380566/274‑1 (MQ=255)
                                                                     ggtgtgggttaaTAGCAAATTCTAGGCAGTAAGGAAGGGGGGGAAGGGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCa                                                                                                                                                                                                                <  4:74974/264‑1 (MQ=255)
                                                                                  tgCAAACTTGAAGGAGTAAGGAAGGGCGCGAAGTGAAGCACAGGGGAAGGTAGGGATTAGTGATGTGTGCGGGGTTGCCCCCGCACTTTTCGGGACGACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGGGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCCGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGtt                                                                                                                                                                                                   <  4:399162/274‑1 (MQ=255)
                                                                                       aCTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCAc                                                                                                                                                                                                                              >  2:297219/1‑243 (MQ=255)
                                                                                                  tAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATCCGGCATTTGCATGATTCGTCGTGTTCTTCCCACGCCATCGCGCGTATAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCCCAGTGGTTTCGATGCCTGTACGGAACCCGCGCCCAATCCCCGCTTTCGCTTACAGCACGGCGGCGTTACGGCAGAAGCCaaa                                                                                                                                                                                    >  4:248358/1‑273 (MQ=255)
                                                                                                   aaGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCCGTGGTTTCGATGCCTGGACGGAACTCGCCCCCAATCCCCGCTTTCCATTTCAGCCCGTCGAGGTTCTGCCAGCAGCCCAccc                                                                                                                                                                                  >  4:660112/1‑274 (MQ=255)
                                                                                                       aaTGGCGCGAAGTGGAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAg                                                                                                                                                                              <  2:357530/275‑1 (MQ=255)
                                                                                                         tGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTTGTTTCGATGCCTGGTCGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCCAGGTTCCGCCAGAAGCCAACCGCCAgc                                                                                                                                                                             >  2:217000/1‑273 (MQ=255)
                                                                                                          ggCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGc                                                                                                                                                                            >  2:27795/1‑274 (MQ=255)
                                                                                                            cgcgAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCCGGCCt                                                                                                                                                                          >  4:242484/1‑274 (MQ=255)
                                                                                                            cgcgAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCt                                                                                                                                                                          >  2:503080/1‑274 (MQ=255)
                                                                                                            cgcgAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCCCTCAATAACCGCTTTCTGTTTCAGATCCTCGCGGTTCTGCCCGCAACCAACCACCAGCCCt                                                                                                                                                                          >  4:573878/1‑274 (MQ=255)
                                                                                                                     tAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCCAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGCTAgg                                                                                                                                                                >  2:277536/1‑275 (MQ=255)
                                                                                                                             gggAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTAc                                                                                                                                                                                                                                                                         >  2:148645/1‑162 (MQ=255)
                                                                                                                             gggAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCt                                                                                                                                                                                                                          <  1:183466/209‑1 (MQ=255)
                                                                                                                             gggAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCt                                                                                                                                                                                                                          >  3:183466/1‑209 (MQ=255)
                                                                                                                                 agggagggATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTc                                                                                                                                                                                                                                                                     <  4:148644/162‑1 (MQ=255)
                                                                                                                                                  gtgtgCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCGCCGCTTTCGCTTTCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCCCCATGCCTGCGAGATAGGCCGCACCCAATAGGGTGACCTCGGCCACt                                                                                                                                   >  4:493719/1‑275 (MQ=255)
                                                                                                                                                  gtgtgCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCCGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCCCGCAc                                                                                                                                    >  2:191770/1‑274 (MQ=255)
                                                                                                                                                   tgtgCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGCACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGCTAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCg                                                                                                                               >  1:130771/1‑278 (MQ=255)
                                                                                                                                                    gtgCGGGGTTGCCCCCGCACGTTTCCGGACCNNGGGATGCTGAATAAACGCTTCATTGGGCATNTACATTATTCGGCGTGTTCTTCCCACGCAATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTANGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCg                                                                                                                               <  3:90949/277‑1 (MQ=255)
                                                                                                                                                     tgCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCACCAGTCCTGCGAGATAGGCCGCACCCACTGCGGTGACTTCGCGCCCTTcc                                                                                                                                >  2:90811/1‑275 (MQ=255)
                                                                                                                                                      gCGNGTNTGCCCCCTCACGTTTCTGGACTACCGGATGCGGCATAAGCGCTTCAGNCGGCATGTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGa                                                                                                                                                                                           <  4:727934/215‑1 (MQ=255)
                                                                                                                                                          ggTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGCAGCCACCCGCCAGGCCTCCGAGCTGGGCCGACCCCAATCCGGTGACTTCGCCAACTCCCGGgc                                                                                                                            >  4:256268/1‑274 (MQ=255)
                                                                                                                                                               cccccGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGCGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCaa                                                                                                                      >  2:121319/1‑275 (MQ=255)
                                                                                                                                                               cccccGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCaa                                                                                                                      >  2:569726/1‑275 (MQ=255)
                                                                                                                                                               cccccGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTCTCCATCCTGCGTAACGGTAATCACGCTCAGTGGTTTCGATGCCTGGACGGAACACGCGCGCAATAACCGCTTTCGCTTGAAGCTCGTCCAGGTTCTGCCAGAAGCCAACCGCCATGCCTGCAAGCAAGGCCGCACCCACTGCGGTCACTACGCGCCCTCCAGGGCGCCCaa                                                                                                                      >  4:67817/1‑275 (MQ=255)
                                                                                                                                                                    gCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGCCTTCGCGCCCTTCCGGGCCCTCAACGCgg                                                                                                                 >  2:372150/1‑275 (MQ=255)
                                                                                                                                                                                  tACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGc                                                                                                                                                            >  2:273488/1‑218 (MQ=255)
                                                                                                                                                                                    ccGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGCTAGGCCGCACCCAATCCGGTGACTTCGCGCACTTCCGGGCCCTCCACCCCGGTCCCGAACATATcc                                                                                                  >  2:117350/1‑273 (MQ=255)
                                                                                                                                                                                    ccGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACCCCGGTGCCGACAatgacgg                                                                                                 >  4:572466/1‑270 (MQ=255)
                                                                                                                                                                                    ccGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCCAGATCGGCCGCACCAAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACCCGGGTGCCGAGaaatgcg                                                                                                  >  4:551210/1‑269 (MQ=255)
                                                                                                                                                                                      ggATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCAccc                                                                                                                                                        <  4:273487/218‑1 (MQ=255)
                                                                                                                                                                                              cATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGt                                                                                                                                                >  2:200550/1‑218 (MQ=255)
                                                                                                                                                                                                  aaCGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACt                                                                                                                                            <  4:200549/218‑1 (MQ=255)
                                                                                                                                                                                                        tCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGTGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGAt                                                                                                                                                                   >  4:388931/1‑189 (MQ=255)
                                                                                                                                                                                                        gcagtgggCATTTTCCTTAGTCGTCTTGGTTTTACCACGCCAACGCGCGTTTAATCTCTTTTTTTCAGCCTGCGTAACGGTAATGGCGCTCAGTGGTGTCGATGCCTTGCCTGAACTCGCGCTCGATCACCGCTTTCTCTTGCAGCTCGTCGCGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTTCCGAGAATATCGGACTGGAACTGCATCAGGaaa                                                                             <  4:290408/269‑1 (MQ=255)
                                                                                                                                                                                                            tCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGTGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGc                                                                                                                                                               <  2:388932/189‑1 (MQ=255)
                                                                                                                                                                                                                 aTTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTGTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCt                                                                    <  2:493720/275‑1 (MQ=255)
                                                                                                                                                                                                                    tACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGTGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACt                                                                 <  2:548188/275‑1 (MQ=255)
                                                                                                                                                                                                                     aCATTTTTCGTGGTGTGCTTCCCACGCCATCGCGCGTTTAACCGCTGTTGTCCAGCCTGCGTAACGGGAATTACGCTCTGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAAGAACCGCTTGTTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTg                                                                <  4:356862/275‑1 (MQ=255)
                                                                                                                                                                                                                          attCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGTGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTTACTTCCCCGCCTTCCGGGCGCTCAACGCGGGTGCCGACAATATCGGACTGGAACTGCATCCGGAAATTTTTTGCTACTGccccc                                                           >  4:548187/1‑274 (MQ=255)
                                                                                                                                                                                                                          attCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCg                                                                                                                                                                                                                                                           >  2:722116/1‑83 (MQ=255)
                                                                                                                                                                                                                          attCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGGTAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGGAAATCGGAATGGAACTGCATCAGGAAATTTTTTGCTTCTGCCcca                                                           >  4:478287/1‑275 (MQ=255)
                                                                                                                                                                                                                          attCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCCCGCCCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCCCTTCCGGGCGCCCCACCCGGGTGCCGCGCATCCCGGCATGGAACTGCCGCCGGACATTGTTTCCCACtacccc                                                            >  4:74891/1‑269 (MQ=255)
                                                                                                                                                                                                                              gtcgtGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGc                                                                                                                                                                                                                                                       <  4:722115/83‑1 (MQ=255)
                                                                                                                                                                                                                              gtcgtGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTa                                                                   <  2:478288/263‑1 (MQ=255)
                                                                                                                                                                                                                                cgtGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCCCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCCCACCCAATTCGGGCACCTCCGGCACTTCCGGGCGCCCAACGCGGGTGCCGAGAATATAGGACGGGAACTGCATCAGGAACCTGTTTGCTACTCAGCCCTCATcc                                                     >  4:4870/1‑275 (MQ=255)
                                                                                                                                                                                                                                   gtattggcccgggcATCGCGCGTTTAACCGCCTGTTGCCAGCCTGCGTAACGTGAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAAGCACCGCTTTCTCTTGGAGCTCGTGGAGGTTGTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCAcg                                                  <  4:191769/262‑1 (MQ=255)
                                                                                                                                                                                                                                     ttttcCCAGGCCATGTGGCGTTGAACCGCATTTNGCAGGCCAGCGAAACGGGAATTACGNTCAGTGGTTTCGGTGCCGGGACGGAACTCGCGCTCAATCACCGCTTTCTCTGGCCGCGCGGCGAGGTTCGGCCAGAAGCCNACCGCCAGGCCGGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAgg                                            <  3:130771/277‑1 (MQ=255)
                                                                                                                                                                                                                                     tcttcCCACGCCATCGGGCGTTTAACCGCTTTTTTCCGGCCTGGGTAACGGGAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCAcgcg                                                <  4:569725/275‑1 (MQ=255)
                                                                                                                                                                                                                                         cccACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCCCCATCCACGCCCAgg                                            >  2:533184/1‑275 (MQ=255)
                                                                                                                                                                                                                                         cccACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAgg                                            >  2:310122/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   cgcgcgTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCCGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAACATATCGGACTGGAACTGCCTCAGGAAATTGTTTGCTCATGCGCCCCCATCCACGCGCAGGGCGTGCCGCCGGac                              >  1:139419/1‑278 (MQ=255)
                                                                                                                                                                                                                                                         tttAACCCCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATAc                            <  2:4870/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAgg                                                                                                                                                                >  2:487826/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                                cGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGc                                                                                                                                                            <  4:487825/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACGGCGCCACCATCCACGCGCCGGGCGTGCAGACGGATACCAGAGTCGGACTGCATCGCTTCCatc  >  1:79041/1‑276 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                     
CGCCGCGTCCTAACCATTTGTAGCCAGCCAGCGCCGCTGATTCGGTGACGCGGGAAAATTCGATGGCAAGTTCTCGTCTCATAGCAAACTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGC  >  NC_000913/4115498‑4116047

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