Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,793,5710TG48.6% ‑3.3 / 17.6 38D571A (GAC→GCC) PP_3354acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (9/10);  new base G (13/5);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCG  >  NC_002947/3793438‑3793705
                                                                                                                                     |                                                                                                                                       
cGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGGGGGCgcggcg                                                                                                                                       >  1:618379/1‑136 (MQ=255)
     cgttgtTCGGCGGCCCGCTGCATGCGATACCCCCCCTGGGCACGCCTGCCGCTGGGGGTGTACATCCCGACCTCGCCGGTGGGTTTCAGACCTTTTCGGAATTTAATGGCAACGGCGTGGCGGCCGCGTCTTTGg‑c                                                                                                                                 <  2:15967/133‑1 (MQ=255)
           tCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGGCGCGGCGCTGg‑ccgcccc                                                                                                                           >  2:861884/1‑134 (MQ=255)
             ggcggACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGATTG‑ccccccccc                                                                                                                         >  2:457375/1‑135 (MQ=255)
                ggACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGGGGCCGCCGCGCTGG‑CCGCccccacgc                                                                                                                      >  1:1313023/1‑132 (MQ=255)
                ggACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGGGGGGGCCGCGGCGATGG‑CCGCCACCTGcc                                                                                                                      <  2:787826/136‑1 (MQ=255)
                  aCCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGGGGGCCCGGCGGTAG‑CCGGCACCAGCCgg                                                                                                                    >  1:676871/1‑136 (MQ=255)
                  aCCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGGGCTGG‑GCGCCCCCAGCCgg                                                                                                                    >  1:851055/1‑136 (MQ=255)
                     gCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGACGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCg                                                                                                                 >  2:15594/1‑136 (MQ=255)
                          aTGCGCTACCGCCCCTCGCCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCAGGTGGGTTTCATAACTATTCGGCAATGAATGGCAACGGCGGGGCGGCCGCGTCTATGG‑CCGCCACCTGCCGGTCGAGGTc                                                                                                            <  2:96030/136‑1 (MQ=255)
                                  ccGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGCCCACCTCCCACGGGGGTTTCACACCGATTCGGAAATTAAGGGCACCGGCGTGGGGGCCGCGTCGAGGC‑CCGCCACCTGCCGGTCGAGGTCGGgcagca                                                                                                    <  2:1152283/136‑1 (MQ=255)
                                   cGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCGCCCCCAGCCGGTCGCGGTCGGgccgcac                                                                                                   >  2:300942/1‑135 (MQ=255)
                                    gCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCGCCCCCTGCCGGTCGAGGTCGGGCAGCAgg                                                                                                  >  2:331796/1‑136 (MQ=255)
                                    gCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCGCGGTCGGGCAGCAgg                                                                                                  >  2:333314/1‑136 (MQ=255)
                                    gCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCCGGCGCGGGCCTGG‑CCGCCCCCGCCCGGGGGGGGGGGGGCCGCAgg                                                                                                  >  1:947828/1‑136 (MQ=255)
                                      cccTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCTATGAATGGCAACGGCGTGGCGGCCGCGTGGATTG‑CCGCCCCCAGCCGGTCGAGGGCGGGCAGCAggga                                                                                                >  2:75845/1‑134 (MQ=255)
                                       ccTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACGTCGCACGTGGGTTTCTCAACGATTCGGCAATTAATGGCAAGGGGGTGGGGGCCGCGTCGAGGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAg                                                                                               <  1:1063245/136‑1 (MQ=255)
                                        cTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGATGG‑CCGCCCCCAGCCGGTCGAGGTCGGGCAGCAGGTAGc                                                                                              >  1:864947/1‑136 (MQ=255)
                                           ggCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGGGATGG‑GCGCCCCCAGCCGGGGGAGGGCGGGCCGCCGGGTGCCgc                                                                                           >  1:578694/1‑135 (MQ=255)
                                            ccACGCCAGCCGCAGCCGGTGTTCATCCCGCCCGCGCACGGGGGTTCCAGACCGATCCGCCAATCAAGGCCAACGGCGGGGCGGCCGCGGCGAGGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                                                          <  1:790537/135‑1 (MQ=255)
                                             cACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGGCGCGGGGGTTG‑GCGCCCCCGGCCGGTGGGGGGCGGGGCGCCGGTTGCCGAgg                                                                                         >  2:532291/1‑136 (MQ=255)
                                                   agccgcagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGGTGG‑CCGCCCCCAGCCGGGCGCGGGCGGGCCGCAGGGAGCCGCGGGAGCAc                                                                                   >  2:744491/1‑136 (MQ=255)
                                                    gccgcagcggGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATTAATGGCACCGGCGTGGCGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAAcc                                                                                  <  1:1231756/136‑1 (MQ=255)
                                                     ccgcagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGGCGCGGCGGGGG‑CCGCCACCCGCCGGGCGCGGTGGGGCCGCAGGGTGCCGCGGGAGGAccc                                                                                 >  1:440415/1‑136 (MQ=255)
                                                      cgccgccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGCCAACGGGGTGGCGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCa                                                                                <  1:20881/136‑1 (MQ=255)
                                                      cgcagccgGTGTTCATGACCACCTCGCAGGTGGGTTTCAGAACGATTCGGCAATGAATGGAAACGGCGTGGCGGCCGCGGCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCa                                                                                <  2:314562/136‑1 (MQ=255)
                                                        cagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGATGG‑GCGCCCCCCGCCCGGCGCGGGCGGGGCGCCGGGTGGCGAGGCAGAACCCAAc                                                                              >  1:1341747/1‑136 (MQ=255)
                                                              gtgtTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTc                                                                        <  2:540915/136‑1 (MQ=255)
                                                                     tGACGACCTCGCACGGGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGGGGCGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCttt                                                                 <  2:937420/136‑1 (MQ=255)
                                                                              cGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCCGGGATGG‑CCCCCCCCAGCCGGTCGACGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCATTACAGCTTTTACCGgccc                                                        >  2:450086/1‑134 (MQ=255)
                                                                                 aCGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCCCGTAGATGG‑CCGCCACCAGCCCGGGGAGGTAGGGCAGCACGTATCCCAGGAAGACCCCAACGCGTTACTGCTTTTACCCGCGCGgc                                                     >  2:1290597/1‑136 (MQ=255)
                                                                                   gTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCCCCACCAGCCGGTCGAGGTCGGGCAGCACGTTGCCGAGGAAGAACCCAACGCTTTACAGCTTTTACCGGCGCAgcgc                                                   >  1:894913/1‑136 (MQ=255)
                                                                                    tGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGGGTGGCGGCCGCGGCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGcgcc                                                  <  1:1298617/136‑1 (MQ=255)
                                                                                             gAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCCCCTCGATGG‑CCGCCACCAGCCGGGCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCGTTCCTGCTTTTACCCGCGCATCGCCCCCCCCAca                                         >  2:629722/1‑135 (MQ=255)
                                                                                               agtaTTCGGAAATGAATGCCATCGGGGGGGCGGCCGCTACTTGGC‑CCCCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTgg                                       <  2:1116407/133‑1 (MQ=255)
                                                                                               aCGATTGGGCAATTAAGGGCAAGGGGGGGGCGGCCGCGTCTAGGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTgg                                       <  1:300942/136‑1 (MQ=255)
                                                                                                cGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCCCCCCCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTACCGGCGCTGCGCCGCCCCCCCAGGc                                      >  1:74039/1‑136 (MQ=255)
                                                                                                                 ccAACGGGGGGGGGGGCGGGGCGGTGGCCCGCCCCCTGCCGGTCGAGGGCGGGCAGCAGGTAGCCGCGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAgc                      <  2:1363527/135‑1 (MQ=255)
                                                                                                                                    gTCGTTGG‑CCGCCACCGGCCGGTCGAGGGCGGGCAGCTGGTCCCCGAGTAAGAACCCAACGCTTTCCTGCTTTCCCCGGCGCTGCGCCGCCCCCACTGCCAACGCCTGGATGCCGCGCGCGCGGCGGGCCCAGGCg  <  2:864947/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGG‑CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCG  >  NC_002947/3793438‑3793705

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