Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,273,8950CA24.4% 66.0 / 13.6 45D40E (GAC→GAA) fabG_53‑oxoacyl‑[acyl‑carrier‑protein] reductase FabG
Reads supporting (aligned to +/- strand):  ref base C (19/15);  new base A (0/11);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGCATGTCAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGTGATATAGGTCCGAT  >  ACB122/3273771‑3274023
                                                                                                                            |                                                                                                                                
caggcaTGTCAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACTCCACTGGCGGCGAGATCGCCAAGCGTTTTGCTCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCTACAAACTGCAGCCg                                                                                                                      >  2:190938/1‑137 (MQ=255)
  ggcaTGTCAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCgct                                                                                                                    <  5:305437/137‑1 (MQ=255)
  ggcaTGTCAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCgct                                                                                                                    >  6:342561/1‑137 (MQ=255)
   gcaTGTCAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCgctg                                                                                                                   >  6:360650/1‑137 (MQ=255)
    caTGTCAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCgctgc                                                                                                                  >  8:120026/1‑137 (MQ=255)
         cAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTggag                                                                                                             <  2:6213/137‑1 (MQ=255)
            aacaaaAAAATGTCGGGCGGGTGATCGGGGCGGGCGACGCCAGTGGCGGCGAGATCGCAATGCGTTTTGCGGGTGTGGGATATATTCCCTGCTTAACCCGTCGCCAGGCCTAAAAACTGCCGCCGCTGCTGTAGaaa                                                                                                          <  4:10921/137‑1 (MQ=255)
                acaaAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATac                                                                                                      >  6:126525/1‑137 (MQ=255)
                    aaaGTCGTGCTGGTGATCGGGGCGGGCGCCGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCTTAACCCGTCGCCAGGCCGACAAAATGCAGCCGCTGCTGGAGGAAATACAcgc                                                                                                  >  7:291757/1‑137 (MQ=255)
                    aaaGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACt                                                                                                                             >  8:352519/1‑110 (MQ=255)
                    aaaGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACt                                                                                                                             <  7:352519/110‑1 (MQ=255)
                    aaaGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACAcgc                                                                                                  >  4:242562/1‑137 (MQ=255)
                                         gcgggctACGCAACGGGCGGCGAGAGCGCCATGCGTTTGGCGCGTGTGGGCTATATTGCGTGCGTAACCCGTCCCCAGGCCGAAAAACTGCAGCCGCTGCTGGAGGAAAAACACGCCGCTGGCGGCCAGGCCTACgg                                                                             <  4:322157/137‑1 (MQ=255)
                                                aCGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGc                                                                      >  2:330487/1‑137 (MQ=255)
                                                aCGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGc                                                                      >  1:194024/1‑137 (MQ=255)
                                                         ggggcgAGATCGCCAAGGTTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGTAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATgcg                                                             <  2:24050/135‑1 (MQ=255)
                                                           ggcgAGATCTCAATGGGTTTTGGGGGTTAGGGATATATTGCCTGCGAACCCCGGCGCCCGGCCTACAAACTGCCGCCGCTGGTGGAGCAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATgcgcg                                                           <  2:167507/137‑1 (MQ=255)
                                                           ggcgAGATCGCCATGCGTTTTGCGCGTTAGGTCTATATTGCCTGCGTAACCCGTCGCCTGGCCTACAAACTGCCGCCGCTGCTGGTGGAAAAACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATgcgcg                                                           <  7:154191/137‑1 (MQ=255)
                                                              gagaTCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCaa                                                        >  8:39636/1‑137 (MQ=255)
                                                               tgaTCGCAAAGTGTTTTGCGGGTGAGGTGTATTTTGCCTGCTTAACCCGTCGCCAGGCCTACAAACTGAAGCCGCTGCTGGAGGAAATACACGCGGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGTTGCGCGCAAg                                                       <  8:280185/136‑1 (MQ=255)
                                                                     ccatgggTTTTGCGGGTGTGTGTTTTATTGCCTGGGAAACCCGGCCCTAGGCCTAAAACGGGCTGCGGGTTGTGTTTGAAAAACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGAAAGGAagag                                                 <  6:7588/131‑1 (MQ=255)
                                                                     ccATGCGTTTGGCGCTGTTGGGCTATATGGCCTGCGTACCCCGCCGCAAGGCCTAAAAACTGCAGCCGCTGCTGGAGAAAATCCACGCCGCGGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAagag                                                 <  8:224779/137‑1 (MQ=255)
                                                                          ggTTTTGGGGGTGAGGGCTATATTGCCTGCTTAACCCTGCGCCAGGCCTAAAAACTGCAGCCGCTGCTGGAGTAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGt                                            <  8:291757/136‑1 (MQ=255)
                                                                          ggTTTTGCGCGTGAGGGCAATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGAAGCCGCTGCTGGAGTAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGt                                            <  6:316135/136‑1 (MQ=255)
                                                                          cGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCTTCGCCAGGCCTAAAAACTGCAGCCGCTGCTGGAGTAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGt                                            <  5:126525/137‑1 (MQ=255)
                                                                            ttttGCGTGTGTGGTATATATTGCGTGCGTACCCCGTGGCCAGGCCTAAAAACTGCAGCCGCTGCTGGTGTAAATCCACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTgg                                          <  8:328430/137‑1 (MQ=255)
                                                                              ttGCGGTTGAGGGATATATGGCCTGCGTACCCCGCCCCCCGCCCCAAAAACTGCCGCCGCTGCTGGATCAAAACCACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGcc                                        <  1:104472/137‑1 (MQ=255)
                                                                                   ggtgtGGGCTTTTTTGCCTGCGTAACCCGTCGCCTGGCCTAAAAACTGCAGCCGCTGCTGGAGTAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACt                                   <  5:250616/136‑1 (MQ=255)
                                                                                      gAGGGCTATTTTGCCTGCGAAACCCGGCCCCAGGCCCACAAACTGCCGCCGCTGCTGTAGTAAACCCACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGt                                <  3:161390/137‑1 (MQ=255)
                                                                                        gggCTATATTGCGTGCGTAACCGGTCGCCTGGCCGACAAAGTGCAGCCGCTGGTGTTGGAAATACCCGCCGCTGGCGGCCAGGCCTACGGTTTTGGCGCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTag                              <  7:56710/137‑1 (MQ=255)
                                                                                             ttatTGCCTGCTAACCCCGGCGCCAGCCCTAAAAACTCCAGCGGCTGCTGGAGAAAAAACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACt                         <  3:242562/136‑1 (MQ=255)
                                                                                                ttGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATc                      <  2:194024/137‑1 (MQ=255)
                                                                                                  gCCTGCGTACCCCGTCGCCTGGCCGACAAACTGCAGCCGCTGCTGGAGTAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGa                    <  6:280521/137‑1 (MQ=255)
                                                                                                  gCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATCCACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGa                    <  1:465/137‑1 (MQ=255)
                                                                                                  gCCGGCGTACCCCGGCGCCTGGCCGACAAACTGCTGCCGCGGGTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGa                    <  2:198297/137‑1 (MQ=255)
                                                                                                      gCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGt                >  1:62611/1‑137 (MQ=255)
                                                                                                      gCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGt                >  6:374829/1‑137 (MQ=255)
                                                                                                      gCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGt                >  8:111000/1‑137 (MQ=255)
                                                                                                      gCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGt                >  1:174505/1‑137 (MQ=255)
                                                                                                        gTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGTGa              >  1:229846/1‑137 (MQ=255)
                                                                                                         tAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCg                     <  7:24047/129‑1 (MQ=255)
                                                                                                         tAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCg                     >  8:24047/1‑129 (MQ=255)
                                                                                                         tAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGTGat             >  5:243985/1‑137 (MQ=255)
                                                                                                             ccGTGGCCAGCCCGAAAAACTGCAGCCGGTGCTGTATGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGTGATATAg         <  6:23413/137‑1 (MQ=255)
                                                                                                                    cAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGTGATATAGGTCCGAt  >  3:326234/1‑137 (MQ=255)
                                                                                                                            |                                                                                                                                
CAGGCATGTCAGAACAACAAAAAGTCGTGCTGGTGATCGGGGCGGGCGACGCCACTGGCGGCGAGATCGCCAAGCGTTTTGCGCGTGAGGGCTATATTGCCTGCGTAACCCGTCGCCAGGCCGACAAACTGCAGCCGCTGCTGGAGGAAATACACGCCGCTGGCGGCCAGGCCTACGGTTTTGGCTCTGATGCGCGCAAGGAAGAGGAGGTGGCCGAACTGGTAGAGACTATCGAGCGTGATATAGGTCCGAT  >  ACB122/3273771‑3274023

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