Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 374,974 A→G 8.6% intergenic (+93/‑485) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913374,9740AG 8.6% 165.7 / 5.7 58intergenic (+93/‑485)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base G (3/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTA  >  NC_000913/374693‑375256
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                            
gATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTg                                                                                                                                                                                                                                                                                       <  2:310085/289‑1 (MQ=255)
 aTCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGt                                                                                                                                                                                                                                                                                      >  2:54828/1‑289 (MQ=255)
                         aCGCTGGTGTTTACTCGGGCTTCCTGCGGCATTGTGAAACGGCGGGGGCGCGTTATGGCGTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTGTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                                                                              <  1:368103/289‑1 (MQ=255)
                                tgtTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCa                                                                                                                                                                                                                                                       >  2:490318/1‑289 (MQ=255)
                                     aCTCGAGCTTCCTGCGTCACTGTGAAACGGCGGGGGCGCGTTATGGCTGAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                                                                                                                                                                                                  <  2:486124/289‑1 (MQ=255)
                                          aGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCCACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTTTCCGGCCTACGGCTGGCGGGAGAATTTGTAGGCCTGATAACACGCGTTCGCGTCTCATCAGGCATCTgc                                                                                                                                                                                                                                              >  1:10388/1‑288 (MQ=255)
                                               ccTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGa                                                                                                                                                                                                                                        <  1:297483/289‑1 (MQ=255)
                                                       aCTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCCCTGGATCTGCGCAACCACAAATAATGATGACTGCCGAGAACGTGCATTTCGTATGCCCATCCCAGCCAGACTGCAGGATGAGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTCGGTCTGATAAGCTGCGTTAGCGTAGCCTGCGGCATCTCCGCACGACTgtccgc                                                                                                                                                                                                                                >  1:202944/1‑284 (MQ=255)
                                                         ttatcaggcctacacttcgcccgtaagccgtaggcctgataagacgcgcaaagcgtcgcatcaggcatggcgcacgactgccggatgcggcgtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgcatagcgtcgcatcaggcatctgcGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGCTAAGACGCGTTCGGGTCGCATCAGGCCTATGCACACGACTGCCGGATg                                                                                                                                                                                                                              <  1:259263‑M1/117‑1 (MQ=255)
                                                                   cggcggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAc                                                                                                                                                                                                                    >  1:38607/1‑289 (MQ=255)
                                                                                      aaGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGa                                                                                                                                                                                                 <  2:38607/289‑1 (MQ=255)
                                                                                       aGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGa                                                                                                                                                                                                 <  1:490318/288‑1 (MQ=255)
                                                                                            gggggATATTCGCGTTGAGCTGGGAAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGGATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                                                                                            <  2:254315/285‑1 (MQ=255)
                                                                                            ggTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                                                                                            <  2:456421/288‑1 (MQ=255)
                                                                                                     cgcgcaaagcgtcgcatcaggcatggcgcacgactgccggatgcggcgtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgcatagcgtcgcatcaggcatctgcGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGA‑‑taaacgtcttgtccggcctacatttcgcccgtaggcagtca                                                                                                                                                                                >  2:259263‑M1/129‑269 (MQ=255)
                                                                                                          ttGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGGATGGCGCGGGAATTTGTAgg                                                                                                                                                                              >  1:292612/1‑288 (MQ=255)
                                                                                                             aGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCt                                                                                                                                                                           <  2:389944/288‑1 (MQ=255)
                                                                                                                 gggCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACACCAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACGGATGGCGTGAGACTTCGTCGGTCTGATAACACGCGTTAGCGTCGCATCCGGCATCCGCGCACGACTGCCGGATGCTGCCTGATCTCCTTACCCGGCCTCCGGCTGGCGCGGGCAttttgccgcctgatc                                                                                                                                                                       >  1:209487/1‑276 (MQ=255)
                                                                                                                             gATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATTAGACGCG‑TTA‑GCg                                                                                                                                                        >  2:12429/1‑289 (MQ=255)
                                                                                                                              aTGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCg                                                                                                                                                        <  2:292612/288‑1 (MQ=255)
                                                                                                                                  gcagatgcctgatgcgacgctaacgcgtcttatcagtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                                                                                                                                                                                                  <  2:272388‑M1/97‑1 (MQ=255)
                                                                                                                                  gcagatgcctgatgcgacgctaacgcgtcttatcagtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                                                                                                                                                                                                  >  1:272388‑M1/100‑196 (MQ=255)
                                                                                                                                    ggcggcCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGCAAA‑GCGTCGCAt                                                                                                                                                  <  2:331662/289‑1 (MQ=255)
                                                                                                                                                  tatGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCAc                                                                                                                                   <  2:10388/289‑1 (MQ=255)
                                                                                                                                                    ggtgcacgccttatgcggtctaaacgtcgcccataaaccggaggccggataagaggcgttagcgtcgcatcaggcatctgcGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACg                                                                                                                                  >  1:209970‑M1/82‑288 (MQ=255)
                                                                                                                                                          ttGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCg                                                                                                                           >  2:110717/1‑289 (MQ=255)
                                                                                                                                                                  gcgcTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCCTCCGCGCACGACTGCCGGATgcggc                                                                                                                   >  2:274133/1‑289 (MQ=255)
                                                                                                                                                                      tGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTg                                                                                                                <  2:167871/288‑1 (MQ=255)
                                                                                                                                                                               cgcccgtaagccgtaggcctgataagacgcgcaaagcgtcgcatcaggcatggcGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACACTTCGCCC‑ATAAACCGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                      <  2:209970‑M1/234‑1 (MQ=255)
                                                                                                                                                                                 aacaacaaATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                     <  2:351559/288‑1 (MQ=255)
                                                                                                                                                                                   caacaaATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc                                                                                                          >  2:248480/1‑281 (MQ=255)
                                                                                                                                                                                   caacaaATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc                                                                                                          <  1:248480/281‑1 (MQ=255)
                                                                                                                                                                                    aaacgataagGATGACTGCCGAGAACGGGCATTTTGTTGGCCCAGCCCAGCACGACGGCCGGATGCGGCGTGAACGCCTTATCCGGCATACGGATGGCGGGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGGGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCAACTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                  <  1:193600/283‑1 (MQ=14)
                                                                                                                                                                                       aaaGAATGATGACTGCCGAGAACTTGCATTTGGTATGCCCATCCCAGCACGACGGCCGGATGCGGCGTCAACGCCTGATCCGGCTTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCAt                                                                                                                                           <  1:351165/244‑1 (MQ=255)
                                                                                                                                                                                              gatgaCTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGCCTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACGGGTgg                                                                                       >  2:434814/1‑289 (MQ=255)
                                                                                                                                                                                                 gaCTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCCTCAGGCATCTGCGCACgc                                                                                                                                 >  2:351165/1‑243 (MQ=255)
                                                                                                                                                                                                  aCTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGcgcg                                                                                   >  2:487445/1‑289 (MQ=255)
                                                                                                                                                                                                    tGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCAACTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTg                                                                                        >  2:193600/1‑282 (MQ=255)
                                                                                                                                                                                                                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                                                                                             >  2:279761/1‑171 (MQ=40)
                                                                                                                                                                                                                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                                                                                             <  1:194569/171‑1 (MQ=40)
                                                                                                                                                                                                                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                                                                                             <  1:279761/171‑1 (MQ=40)
                                                                                                                                                                                                                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                                                                                             >  2:194569/1‑171 (MQ=40)
                                                                                                                                                                                                                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAg                                                     >  2:53249/1‑289 (MQ=255)
                                                                                                                                                                                                                                    aGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGACTTCTCCGGCCTACGGGTGGCGCGGGCATTTGTAGGCCTGATAAGACGCGTTAGCGTc                                                 >  2:428240/1‑289 (MQ=255)
                                                                                                                                                                                                                                          ACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCAt                                                                                                                                                                                                                                                         <  2:362824‑M1/85‑1 (MQ=255)
                                                                                                                                                                                                                                          ACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCAt                                                                                                                                                                                                                                                         >  1:362824‑M1/1‑85 (MQ=255)
                                                                                                                                                                                                                                          ACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGCCCGT‑AAGCCGTAGGCCTGATAAGACGCG‑CAAAGCGTCGCATc                                                                                                                                                 <  2:343524‑M1/187‑1 (MQ=255)
                                                                                                                                                                                                                                          ACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGCCCGT‑AAGCCGTAGGCCTGATAAGACGCG‑CAAAGCGTCGCATc                                                                                                                                                 >  1:343524‑M1/1‑187 (MQ=255)
                                                                                                                                                                                                                                           CTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                                                                                                                     <  2:76556‑M1/214‑1 (MQ=255)
                                                                                                                                                                                                                                           CTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                                                                                                                     >  1:76556‑M1/1‑214 (MQ=255)
                                                                                                                                                                                                                                                        cgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                             <  2:355642/257‑1 (MQ=40)
                                                                                                                                                                                                                                                        cgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                             >  1:355642/1‑257 (MQ=40)
                                                                                                                                                                                                                                                                                    gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCCTCTGCGCACGACTGCCGGATGCGGAGTTAACGCCTTATCCGGCCCACGGCTGGCCCGAGAATTTGTAGGCCTGGTAAGACGCCTTAGCGTCGCAGCAGGCCTACTCGCACGGCCTCCGCTTCCGTCGTAACCGCCttg  >  1:149046/1‑287 (MQ=11)
                                                                                                                                                                                                                                                                                    gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGAAGCGTTAGCGTCGCATCAGGCTTCCTCGCCCGACCGCCGGGTgcgcgcaaaacgccgtt  >  1:294439/1‑274 (MQ=14)
                                                                                                                                                                                                                                                                                    gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATAAGGCGTTCACGCCGCATCCGGcagtcgt                                                                                                                                                                                                  >  2:368948/1‑92 (MQ=11)
                                                                                                                                                                                                                                                                                    gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATAAGGCGTTCACGCCGCATCCGGcagtcgt                                                                                                                                                                                                  <  1:368948/98‑7 (MQ=11)
                                                                                                                                                                                                                                                                                       tGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAt                   <  1:52066/268‑1 (MQ=40)
                                                                                                                                                                                                                                                                                       tGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAt                   >  2:52066/1‑268 (MQ=40)
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                            
GATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTA‑GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTA  >  NC_000913/374693‑375256

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