Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,554,0000TG52.2% 3.5 / 14.4 23Y557S (TAC→TCC) maeAmalate dehydrogenase (oxaloacetate‑decarboxylating)
Reads supporting (aligned to +/- strand):  ref base T (1/10);  new base G (12/0);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 9.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGGCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCG  >  NC_000913/1553755‑1554221
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                             
cagGGCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTTGCGGGATTCg                                                                                                                                                                                                                          >  1:169205/1‑251 (MQ=255)
            aTGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTACCAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTCGTCGCGGTAATCGGGTTTCCAgca                                                                                                                                                                                                               >  2:119446/1‑248 (MQ=255)
                      gAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCa                                                                                                                                                                                                     <  2:104343/250‑1 (MQ=255)
                           cAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATAGTCAATGGcc                                                                                                                                                                                               >  1:280554/1‑251 (MQ=255)
                              tgtTTTATAAGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGCCGCGGGATTCGGGTTGCCCGGACATATAGTCAATGGCCTgg                                                                                                                                                                                            >  2:143083/1‑250 (MQ=255)
                                   tataCGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCa                                                                                                                                                                                        <  1:238735/250‑1 (MQ=255)
                                           aCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGCTGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGTGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAATTATCGTCAATGGCCTGTTGGAGGGCTTCgg                                                                                                                                                                               >  1:145213/1‑251 (MQ=255)
                                                                  aacaacTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAATTATCGTCAATGGCCTGGTGCAGGGCTTCGGCCGAGGGTTTCAccgccccgccc                                                                                                                                                        >  1:171434/1‑250 (MQ=255)
                                                                    caacTGCAGCCTTCATAGATACTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAACTGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCCGGTTGGCAGAAAATATTGGCAATTGCCTGTTGCCGGGATTCGGCCGAGGGTTTCACCCCCACGccctg                                                                                                                                                      >  1:99627/1‑248 (MQ=255)
                                                                      acTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAACTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCCATGGCCCGTTGCCGGGCTTCCGCCGAGGTTTTTACCGCCAcccccttct                                                                                                                                                    >  2:296999/1‑243 (MQ=255)
                                                                                    aGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCCttttt                                                                                                                                      >  1:180326/1‑250 (MQ=255)
                                                                                    aGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCCGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCCttttt                                                                                                                                      >  1:39309/1‑250 (MQ=255)
                                                                                                                    cgggggAGGAAAATGAGATTTATTCGCCACTCCCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTCCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAGGGCCTGTGGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACttt                                                                                                      <  2:209961/246‑1 (MQ=255)
                                                                                                                                                                 agaTCTGAAAAAGGGAGAGGGAAAGAGCCCGGTAGCCTTCACTCCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGTTATTCGGCTGGCCAGAAATTACCGTCAATGGCCGTTTGCAGGGCTTCGGAAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGTAGTACCATACCTTCGCCGTTCAg                                                         <  2:106246/251‑1 (MQ=255)
                                                                                                                                                                         aaaaGGGGGAGGAAAATTGCCCGGTAGCCTTCACTCCCGGGCGCAGGCTGAGATGGAGGTACGGCGGAAGTCGCGGTGTTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGTCTTCGGGAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTGTGCCAACCGCAAACGCAATTGCGCGGGATACTTTCTGAATAGCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATg                                                 <  1:41872/251‑1 (MQ=255)
                                                                                                                                                                         aaaaGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCCAACGCAATTGCGCGGGGGACTTTCTGAATATCTTTCAGTTCCGGCAGGACCATACCTTCGCCGTTCAGCACCAATg                                                 >  1:152637/1‑251 (MQ=255)
                                                                                                                                                                              ggAGGGGGAAATAGCCGGGTAGCCTCCACTCCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGGGGTATTCGGCTTGCAAGAAATTATCGTCAAGGGCCTGGTGCGGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGaa                                            <  2:180326/251‑1 (MQ=255)
                                                                                                                                                                              ggAGAGGGAAATAGCCGGGTGCCCTTCCCTCCGGGGCGCAGGTTAAGATGGGGGTAGGGCGGTATTCGGGGTATTTGGTTTGCCAGAAATTATCTTCAATGGCCTGTTGAAGGGTTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTGCAGTTCCGGCAGTGCCATACCTTCGCCGTTCAGCACCAATGGTGaa                                            <  1:119446/251‑1 (MQ=255)
                                                                                                                                                                                 gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCa                                                                                                                                                                                        <  2:321263/108‑1 (MQ=255)
                                                                                                                                                                                 gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATTGTTAATGGCCTGGTGcc                                                                                                                                                                                        >  1:321263/1‑107 (MQ=255)
                                                                                                                                                                                                  aGCCTTCACTCCCGGGCGCAGGCTTAGATGGGGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCGTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACt                        <  1:196630/251‑1 (MQ=255)
                                                                                                                                                                                                      ttCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTTACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTTCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCCTACCTTCGCCGTTCAGCACCAATGGTGAATACCGCGCCAGCGTTTCACTTGcc                    >  1:191285/1‑251 (MQ=255)
                                                                                                                                                                                                                          tAGAGGGAGGAACGGGGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAGGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCg  <  1:298482/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                             
CAGGGCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCG  >  NC_000913/1553755‑1554221

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