Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,737,4430TC36.1% 46.1 / 17.0 38I229I (ATT→ATCpaoBAldehyde oxidoreductase FAD‑binding subunitPaoB
Reads supporting (aligned to +/- strand):  ref base T (19/4);  new base C (0/13);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.03e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGGCATTCCATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTT  >  ACB143/3737316‑3737572
                                                                                                                               |                                                                                                                                 
ttGCGGCATACCATCGCCCCCCCCAGGACACCCCGGAGCTACAAAATCAGCTGGAAATCTATCAAGTGGTGGGCGCAATCGAGCTGCCGCCCCCGGTGTTCCAAAAACATTACGTTTATATTAAAATTCGCGAcc                                                                                                                            <  4:235607/135‑1 (MQ=255)
  gCGGCATTCCATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACcgcg                                                                                                                         >  3:357047/1‑136 (MQ=255)
          ccATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTAc                                                                                                                 >  8:328721/1‑136 (MQ=255)
                      cGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTg                                                                                                     >  1:72337/1‑136 (MQ=255)
                                ccGGACCGATACATTCCGCTGCCATTCGATTAACTGGTGGGCGCAATGGAGGGGCCGCCCCCGTGGTCCAAAAAACAATACCTTAACCTTAAAACTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTcg                                                                                           <  8:1312/136‑1 (MQ=255)
                                  ggAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTcgcc                                                                                         >  3:29969/1‑136 (MQ=255)
                                                 gCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTc                                                                          >  6:309486/1‑136 (MQ=255)
                                                     ccAATCGATGAAGTGGTGGGCGCCATCGAGCTGCCGCCCCGGGTGTTCTAAAAAAATTACCTTTATCTTAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATg                                                                      <  5:99860/135‑1 (MQ=255)
                                                      cAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTACTACGCGTTTGCGCTGTTATCTGTCGCCGCAGCCTTGAGCTTCGATgg                                                                     >  2:97010/1‑136 (MQ=255)
                                                           gATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAgg                                                                >  7:15017/1‑136 (MQ=255)
                                                             tGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGAGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTAGACCTTCGATGGCCAGGt                                                               >  8:151237/1‑135 (MQ=255)
                                                               aaCTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGAt                                                            >  1:240486/1‑136 (MQ=255)
                                                                          gCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCACGCTc                                                 >  2:267353/1‑136 (MQ=255)
                                                                          gCCAGCGTGCTGCCCCCCCCGGTTTAATAAAAAAATAACCTTTAACTAAAATTACGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTc                                                 <  7:279741/136‑1 (MQ=255)
                                                                            cATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGc                                               >  6:346778/1‑136 (MQ=255)
                                                                            cATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGc                                               >  4:356166/1‑136 (MQ=255)
                                                                               cgcgcggcccccCCCGGCGTACCAAAAACAAACCCGTAACCAAAAAACCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTg                                            <  3:270239/130‑1 (MQ=255)
                                                                               cGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTg                                            >  5:297112/1‑136 (MQ=255)
                                                                                         ccccgGGTTTACAAAAAACATAACCTATATCTAAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGggcg                                  <  6:21612/136‑1 (MQ=255)
                                                                                          cccgGGTGTCCTAAAAACACACCCGTTCTCTTACATCCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGggcgg                                 <  3:221152/136‑1 (MQ=255)
                                                                                          ccccGGTGTTCCAAAAACATAACCGTTATCTCAAAAACCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGggcgg                                 <  6:126169/136‑1 (MQ=255)
                                                                                                  ttCCAAAAAAATACCCTTAATCTAAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCt                          <  4:357047/135‑1 (MQ=255)
                                                                                                      aaaaaCATACCCTTAATCTAAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACaa                     <  3:66686/136‑1 (MQ=255)
                                                                                                       aaaaCATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAg                    >  4:318072/1‑136 (MQ=255)
                                                                                                            ataaCCTTTATCTTAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTg               <  2:19818/133‑1 (MQ=255)
                                                                                                            aTTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTg               >  1:88209/1‑136 (MQ=255)
                                                                                                               aCCTTTATATCAAAAAACGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCg            <  3:132636/136‑1 (MQ=255)
                                                                                                               aCCTTTAGCTGAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCg            <  5:122028/136‑1 (MQ=255)
                                                                                                               aCCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGc             >  8:341605/1‑135 (MQ=255)
                                                                                                                ccGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCg            <  8:222040/135‑1 (MQ=255)
                                                                                                                 gtttttctcAAATTCCGCGACCGCGCCTCCTCCGCGTTTGCGCTGGTATCGGGCGCCGCAGCCTTGAGCTTCGATGGCCATGGGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag          <  5:60312/131‑1 (MQ=255)
                                                                                                                 cttaaTCTCAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag          <  4:318410/132‑1 (MQ=255)
                                                                                                                 cGTTTTCTTAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag          <  3:97638/136‑1 (MQ=255)
                                                                                                                 cGTTATCTCAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag          <  7:328721/136‑1 (MQ=255)
                                                                                                                      tCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCg     >  1:161943/1‑136 (MQ=255)
                                                                                                                        tGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGt   >  3:357709/1‑136 (MQ=255)
                                                                                                                        tGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGt   >  2:346369/1‑136 (MQ=255)
                                                                                                                          aaaaTCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGtt  <  4:203686/135‑1 (MQ=255)
                                                                                                                               |                                                                                                                                 
TTGCGGCATTCCATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTT  >  ACB143/3737316‑3737572

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: