Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,934,7140GC30.0% 63.0 / 12.2 40A142P (GCC→CCC) aarL‑alanine/L‑glutamate racemase
Reads supporting (aligned to +/- strand):  ref base G (22/6);  new base C (0/12);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 3.32e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACATTCGTGAAACCAGGTGATCACGTACTTGCTGGAGATACCATCGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTG  >  KT2440/4934594‑4934848
                                                                                                                        |                                                                                                                                      
acaTTCGTGAAACCAGGTGATCACGTACTTGCTGGAGATACCATCGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTaaaaaa                                                                                                                        >  6:104950/1‑137 (MQ=255)
  aTTCTTCAAACCCGGTGATCACTTAGTTGCGGTATATACCATCTTTAAACCCACCCGTTGGCTGCTTAATCATTTTTAAAATATTGGTTAAAAAATTCAAACCTACACACCACGCCACCCCGCGACTCTAAAAAATa                                                                                                                      <  4:159141/137‑1 (MQ=255)
    tCGTGAAACCAGGTGATCACGTACTTGCTGGAGATACCATCGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATaaa                                                                                                                    >  5:2991/1‑137 (MQ=255)
                   aTCACGTACTTGCTGGAGATACCATCGTTTACTACACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACa                                                                                                     >  4:124759/1‑137 (MQ=255)
                            ttGCTGGAGATACCATCGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAAta                                                                                            >  2:274667/1‑136 (MQ=255)
                                 ggataaaCCATAGTTAACTCCACACGATGGCTGCTCAATCACTTTTAAAATATTGGTTTAAAAATTCAAAACTTCACACCTCCCCAACCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATAtt                                                                                       <  5:103502/131‑1 (MQ=255)
                                   agaTACCATCGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATAtttt                                                                                     >  6:13621/1‑136 (MQ=255)
                                       aCCATCGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGaaaa                                                                                 >  2:30548/1‑137 (MQ=255)
                                            cGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAg                                                                            >  6:349028/1‑137 (MQ=255)
                                              tttACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGgc                                                                          >  1:395728/1‑137 (MQ=255)
                                                aaCTCCACCCGAGGGGTGGTTAAAAATTTTTAAATTTTTGTTTAAAAAATTTAATACTAACCACCTAGCGACCCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGgcgc                                                                        <  6:2991/136‑1 (MQ=255)
                                                aaCTCCACACGATGGCTGCTTAATCACTTTGAAAATATTGGTGAAAAAATTTAAAAATTCACACCTAGCGAACCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGgcgc                                                                        <  4:300010/136‑1 (MQ=255)
                                                      acacGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCAc                                                                  >  3:359520/1‑137 (MQ=255)
                                                       aactaTGGCTTCTTAATCAATTTTAAAATATTTGTTTAAAAATTTAATACTTCACACATAGCCAACCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACg                                                                 <  5:77384/133‑1 (MQ=255)
                                                        acGATGGCTTCTCAATCACTTTTAAAATTTTGGTTTAAAATTTTAATACTTCACACAAAGCCAACCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGt                                                                <  3:430484/137‑1 (MQ=255)
                                                          tattcctgctgAATCACTTTGAAAATATTGTTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTAc                                                              >  5:350888/6‑137 (MQ=255)
                                                          gATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTAc                                                              >  6:338302/1‑137 (MQ=255)
                                                           ttGGCTGCTAAATCACTTTTAAAATATTGGTTAAAAAATTGAATACAGCACACCTAGCCACGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACg                                                             <  6:146044/136‑1 (MQ=255)
                                                             gggtgcttAATCACTTTAAAAATATTGTTTTAAAAATTAAATACTTCACACCTCGCGACCCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAg                                                           <  7:61226/134‑1 (MQ=255)
                                                             ggCTTCTAAAACAATTTAAAAATATTGTTTAAAAAATATAATAATACACACCTAGCCACCCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAg                                                           <  3:51727/137‑1 (MQ=255)
                                                               tttattaaTCATTTTAAAAATTTTGGTTAAAAAATTTAAAACTCCACACATCGCCAACCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGt                                                         <  2:232039/131‑1 (MQ=255)
                                                               gtgcttAACCAATTTTAAAATTTTGGTAAAAAAATACAATACTTCACACCTCGCCACCCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGt                                                         <  4:284257/136‑1 (MQ=255)
                                                               gtgctgAAATATTTTTAAAATATTGTTGAAAAAATTTAAAATTACACACCTTGCACCCCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGt                                                         <  2:236811/136‑1 (MQ=255)
                                                               ctgcttAATCACTTTCAAAATATTGGTTTAAAAATTGAATACTTCACACCTAGCCAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGt                                                         <  8:85676/137‑1 (MQ=255)
                                                               aagcttAATCACTTTTAAAATATTGGTTTAAAAATTGAATACTACACACCTAGCCAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGt                                                         <  8:152660/135‑1 (MQ=255)
                                                                tgctgAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTc                                                        >  2:406252/1‑137 (MQ=255)
                                                                  ctgAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCAt                                                      >  7:141420/1‑137 (MQ=255)
                                                                       tCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGata                                                 >  1:232039/1‑137 (MQ=255)
                                                                             tGAAAATATTGGTGTAAAAAATGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTAGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTg                                           >  2:56924/1‑137 (MQ=255)
                                                                                aaaTATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGctg                                        >  7:152660/1‑135 (MQ=255)
                                                                                 aaTATTTGTTTAAAAATTTAATAATTAAAAACACGCTAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGAt                                       <  3:77373/137‑1 (MQ=255)
                                                                                 aaTATTGTTTAAAAAATTTAAAACTTCACACCTAGCTAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGAt                                       <  5:104950/137‑1 (MQ=255)
                                                                                                aTTTAATACTTAACACGTAGCGACCCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGaca                        <  5:334397/137‑1 (MQ=255)
                                                                                                aTTTAATACTACACACCTAGCAACGCCGCTACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTCCGAGGTAATGTATATTAGTGCGATCGTCGAATCAAGaca                        <  6:350888/137‑1 (MQ=255)
                                                                                                            acacCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAgc            >  4:33587/1‑137 (MQ=255)
                                                                                                              acCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAgcgc          >  4:18034/1‑137 (MQ=255)
                                                                                                                     gAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATtgt   >  2:250804/1‑137 (MQ=255)
                                                                                                                     gAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATtgt   >  7:278731/1‑137 (MQ=255)
                                                                                                                      aaGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATtgtg  >  3:479118/1‑137 (MQ=255)
                                                                                                                      aaGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACACAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACAACCCAGCGCTGATtttg  >  8:172187/1‑137 (MQ=255)
                                                                                                                        |                                                                                                                                      
ACATTCGTGAAACCAGGTGATCACGTACTTGCTGGAGATACCATCGTTTACTCCACACGATGGCTGCTGAATCACTTTGAAAATATTGGTGTAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTG  >  KT2440/4934594‑4934848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑