Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,793,5710TG39.5% 16.5 / 11.3 43D571A (GAC→GCC) PP_3354acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (14/12);  new base G (10/7);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGG  >  NC_002947/3793438‑3793703
                                                                                                                                      |                                                                                                                                    
cGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGGGGCCgcggcg                                                                                                                                    >  1:117057/1‑136 (MQ=255)
cGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACAGCCCCTCCGCACGCCAGCCGCAGCCGGTGTTCATGACGACCACGCACGTGGGTTTCAGAACGATACGGCAATGAATGGCAAC‑GGCGTGGGGGCCgcggcg                                                                                                                                    >  1:702614/1‑136 (MQ=255)
 gCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGGCgcggcgg                                                                                                                                   >  1:860024/1‑135 (MQ=255)
           tCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGGCGCTGgccgcccc                                                                                                                         >  1:290650/1‑134 (MQ=255)
           tCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGGCGCTGgccccccc                                                                                                                         >  1:651383/1‑134 (MQ=255)
           tCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGGAGGAGGCCGCCAc                                                                                                                         >  1:1136940/1‑136 (MQ=255)
                      cTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGa                                                                                                              <  2:455237/136‑1 (MQ=255)
                          aggcgcTACCGCCCCGCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGGGGGTTTCAGAACGATTCGGCAATGAATGCCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTc                                                                                                          <  1:18316/134‑1 (MQ=255)
                                 aCCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCCCCTGCCGGGCGCGGTCGGgcagc                                                                                                   >  2:673103/1‑136 (MQ=255)
                                   cGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCCGGGCTGGGCGCCCCCAGCCGGTCGAGGTCGGgcagcac                                                                                                 >  2:977901/1‑135 (MQ=255)
                                    gCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCCGCGATGGCCGCCCCCGGCCGGGCGCGGTGGGGGaggagg                                                                                                >  1:597403/1‑136 (MQ=255)
                                           ggCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCAGCCGGTCGAGGTCGGGCAGCAGGTAGCCgc                                                                                         >  2:21641/1‑135 (MQ=255)
                                            gCACGCCAGCCGCAGCGGGTGTTCTTGACGACCTCGCAGGTGGGTTTCAGAACGATTCGGCATTGAATGGCAAC‑GGGGGGGGGGCCGCGGCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                                                        <  1:986925/136‑1 (MQ=255)
                                            ccACGCCTGCCGCCGCCGTTGTTAATGCCGACCGCGCACGGGGGTTTCAGACCGATTGGCAAATGATGGGCAAG‑GGCGGGGCGGCGGGGTCGTTGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                                                        <  2:97531/135‑1 (MQ=255)
                                            ccACCCCAGCCGCAGCGGGTGTTCATGACGCCCTCGCACGTGGGTTTCACACCGATGCGGAATTCAATGGCACC‑GGGGGGGGGCCGGGGGCTAGGCCCCCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                                                        <  2:200142/135‑1 (MQ=255)
                                             cACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGC‑GCCGCGTCGAAGGCCGCCCCCTGCCGGTCGAGGTGGGGCCGCAGGTATCCGAGga                                                                                      >  2:5076/1‑136 (MQ=255)
                                                 ccAGCCGCAGCCGGTGTTCATGACGACCCCCCACGGGGGTTTCAGAACGATTCGGCAATCAATGGCACC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGgaaga                                                                                   <  1:618804/136‑1 (MQ=255)
                                                   agccgcagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGGCGGTGGCCGCCCCCAGCCGGGCGCGGGCGGGCAGCAGGGATCCGCGGAAGGAc                                                                                 >  2:93514/1‑136 (MQ=255)
                                                     ccgcagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGGCGATGGCCCCCCCCAGCCCGGCGAGGGCGGGCCGCCGGTTGCCCGGGGAGAAccc                                                                               >  2:1054763/1‑136 (MQ=255)
                                                                gtTCATGACGCCCTCGCACGTGGGTTTCAGAACGATGCGGCAATGAATGGCAAC‑GGCGGGGCGGCCGCGGCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCt                                                                    <  1:88769/136‑1 (MQ=255)
                                                                gtTCATGACGACCGCGCACGGGGGTTTCACAACTATTCGGCAATGAATGCCACC‑GGGGGGGCGGCCGCGTCTAGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCt                                                                    <  1:518625/136‑1 (MQ=255)
                                                                     tcccccccTCCCCCGGGGTTTACACAACGATGCGCAAATTAATGCCACC‑GGCGGGGCGGCCGCGGCGAGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCttt                                                               <  1:1003344/130‑1 (MQ=255)
                                                                     tGACGACCTCGCAGGGGGGTTTCACAACTATTCGGCAATTATTGGCACC‑GGCGTGGCGGCCGCGTCTATGGCCCCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCttt                                                               <  2:701297/136‑1 (MQ=255)
                                                                     cgaccacCTCGCCCGGGGGTTTCAGACCGATTGGGCAATGAAGGGCACC‑GGCGTGGGGGCGGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCttt                                                               <  2:1150284/135‑1 (MQ=255)
                                                                            cTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTGGATGGCCGCCCCCTGGCGGTGGAGGTAGGGCAGCAGGTTGCCGCGGTAGAACCCAAAGCGTTACTGGTTTTACCggg                                                        >  2:559852/1‑135 (MQ=255)
                                                                            cTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCAGCCCGTCGAGGTCGGGCAGCAGGTTGCCGCGGAAGAACCCAACGCTTTACAGCTTTTACCCgc                                                        >  1:313348/1‑136 (MQ=255)
                                                                            cTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCCCGGGGGTTGCCCCCCCCAGCCGGGCGAGGGCGGGCCGCACGTTTCCGCCGAAGAACCCAACCCTTTACAGGTTTTACCCgc                                                        >  1:47092/1‑136 (MQ=255)
                                                                              cGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGGCGCTGGCCGCCCCCAGCCCGTCGAGGTCGGGCAGCAGGTAGCCGCGGCAGAACCCAACCCGTTGCTGCTTTTGCCCgcgc                                                      >  1:106683/1‑136 (MQ=255)
                                                                               ccccGTGGGTTTCACAACGATCCGCCAATGAAGGCCAAC‑GGCGGGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCt                                                     <  2:3152/133‑1 (MQ=255)
                                                                                 aCGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCCGCGCTgc                                                   >  1:602592/1‑136 (MQ=255)
                                                                                  ggggggTTTCAGACCGATGCGGCAATGAAGGGCAAC‑GGCGGGGCGGCCGCGGCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTgcg                                                  <  2:3178/133‑1 (MQ=255)
                                                                                       gTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTGGATGGCCGCCACCAGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGGAGAACCCAAAGCATTACAGCTTTTACCGGCGCTTCgccgcc                                             >  2:846431/1‑136 (MQ=255)
                                                                                       gTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCgccgcc                                             >  1:287406/1‑136 (MQ=255)
                                                                                            aGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCAc                                        <  2:1042451/136‑1 (MQ=255)
                                                                                             caaCGATTCGGAAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGTAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACt                                       <  1:478478/135‑1 (MQ=255)
                                                                                             caaCGATGCGGTAATGAAGGGCAAC‑GGCGTGGCGGCCGCGTCGAGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACt                                       <  2:875977/135‑1 (MQ=255)
                                                                                             aaaCGATTCGGCAAT‑CATGGCCACGGGGGGGGCGGCGGCGGCGAGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACt                                       <  1:673103/135‑1 (MQ=255)
                                                                                              aaCGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCCACCCTTTCCAGCTTTTACCGGCGCAGCGCCGCCCCCAcag                                      >  2:1148486/1‑134 (MQ=255)
                                                                                                    tCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCCTCGATGGCCCCCACCTGCCGGTCGACGTCGGGCAGCAGGTAGCCCAGGAAGAACCCAACCCGTTCCTGCTTTTCCCGGCGCAGCCCCGCCCCCCCAGGCAAcc                                >  2:392190/1‑135 (MQ=255)
                                                                                                      ggCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGAAGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAAACCAACGCTTTACTGCTTTTCACCGCGCAGCGCCCCCCCCACAGGCAACGcc                              >  2:420972/1‑136 (MQ=255)
                                                                                                                     c‑GGCGTGGCGGCCGCGGCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCCCAGGCAACGCCAGGATGCAgcgcgcc               >  1:460926/1‑135 (MQ=255)
                                                                                                                                gcCGCGTCGATGGCCGCCCCCAGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTGCTGCTTTTACCCGCGCTGCGCCGCCCCCCCAGGCAACCCCGGGATGCAGCGCGCGCGGCGGGccc     >  2:824457/1‑136 (MQ=255)
                                                                                                                                gcCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGccc     >  2:1035276/1‑136 (MQ=255)
                                                                                                                                  cgcgGCGTTGCCCGCCCCCGGCGGTTCGCGGTCGGGCAGCAGGTAGCCGCGGAAGAACCCACCGGTTTCCTGCTTTTCCGGGCGGGGGGCCGCCCCCCCTGCCAACGCGTGGATGCAGCGCGCGCGGCGGGCCCAg   <  2:452052/136‑1 (MQ=255)
                                                                                                                                   gcgTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTACTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCGGCGGGcccccg  >  1:22723/1‑133 (MQ=255)
                                                                                                                                   gcgTCGAGGGCCGCCACCTGCCGGTCGAGGGCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGCCAACGCCTGGATGCCGCGCGCGCGGCGGGCCCAgg  <  2:1095453/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAAC‑GGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGG  >  NC_002947/3793438‑3793703

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