Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,253,4660AC40.5% 23.0 / 12.0 44L178R (CTG→CGG) PP_3726enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (8/17);  new base C (15/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGAA  >  NC_002947/4253336‑4253599
                                                                                                                                  |                                                                                                                                     
gcCGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGgcg                                                                                                                                  >  2:8821/1‑136 (MQ=255)
 ccGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGgcgc                                                                                                                                 >  1:6474/1‑136 (MQ=255)
 ccGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGgcgc                                                                                                                                 >  1:222213/1‑136 (MQ=255)
 ccGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGgcgc                                                                                                                                 >  1:498929/1‑136 (MQ=255)
 ccGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGgcgc                                                                                                                                 <  2:230468/136‑1 (MQ=255)
 ccGTCGCATATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGGTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCggggcc                                                                                                                                 >  1:83057/1‑133 (MQ=255)
              gCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCca                                                                                                                    <  1:254861/136‑1 (MQ=255)
                 gCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGGGCCCCACaa                                                                                                                 >  1:213466/1‑136 (MQ=255)
                  gttGGTCTGCGCATAGGCGAAGGGGGGGCCATTGGCCAGGCGCCCGGCCAGTGCGGCGGCACCGGCAGGCAACTGCGGGGCGGCGCAGACTTCGCTCACCAACCGGGCGGCCAGGGCGCGATCGGCGCCCCACAAc                                                                                                                <  2:441292/135‑1 (MQ=255)
                     ggTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCCATCGGCGCCCCACAACtcc                                                                                                             >  1:497830/1‑135 (MQ=255)
                              aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACCACTCGTCGAGGAAc                                                                                                    >  1:568231/1‑136 (MQ=255)
                              aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAc                                                                                                    >  2:487115/1‑136 (MQ=255)
                              aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAc                                                                                                    >  1:396489/1‑136 (MQ=255)
                               ccggCGAAGCTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAca                                                                                                   <  1:132309/134‑1 (MQ=255)
                                   cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGAACGGCGCCCCCCAACTCGGCGCGGCAccccag                                                                                               >  1:429354/1‑131 (MQ=255)
                                   cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTCGGCGCCCCACAACTCGTCGGGGGAcagccg                                                                                               >  2:10934/1‑134 (MQ=255)
                                    gAAGGTCGGGCCATTGGCAAGGCGCCCGGCCTGGTCGGCGGCAACGGCAGGCAACTGCTGGGCGGCGCAGACTTCGCGGACCAACCCGGCGGCCAGGGCGCTATGGGCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                              <  2:382320/136‑1 (MQ=255)
                                               cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGCCCCACCCGGCGGCCCGGGCGCGCTCGGGGCCCCCCCACTCGGCGGGGGACCGCCGGCGCGTGGccgg                                                                                   >  2:91764/1‑134 (MQ=255)
                                               cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACTCGTCGAGGAACCGCAGGCGCTTGGccgg                                                                                   >  1:362375/1‑134 (MQ=255)
                                               cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAAGGCGCGAACGGCGCCCCCCAACTCGACGACGAACAGCACGCGCTTGGccag                                                                                   >  1:235605/1‑134 (MQ=255)
                                                     ccAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACAACGCGTCGAGGGACACCCCGCCCTTGGGCTTCTCGcg                                                                             >  1:517230/1‑135 (MQ=255)
                                                     ccAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCt                                                                             >  1:326248/1‑136 (MQ=255)
                                                        ggCGCCGGGCATGTGCGGCGGCAACGGCTGGCAAGTGCGCGGGGGCGCATACTTCGCTCACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                          <  1:351595/136‑1 (MQ=255)
                                                                      tCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCa                                                            >  1:193805/1‑136 (MQ=255)
                                                                        ggcggcAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAgg                                                          >  1:384153/1‑136 (MQ=255)
                                                                                        aaCTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                                          <  2:193805/136‑1 (MQ=255)
                                                                                        aaCGGCGCGGCGGCGCAGACTGCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                                          <  1:311903/136‑1 (MQ=255)
                                                                                          ctgcgggtcgGCGAATACTTCGCTGACCACCCCGGCGGCCCGGGCGCGAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                        <  1:161491/130‑1 (MQ=255)
                                                                                          cTGCTCGTCGGCGCAGACTTCGCTGACCACCCGGGCGCCCAGGCCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                        <  2:346872/136‑1 (MQ=255)
                                                                                               ggtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                                   <  2:481138/135‑1 (MQ=255)
                                                                                                gtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCCGGCGCTTGGCCTGCTCGCTGCCGCTCAGCCGGGGGCGGTGGCCGCGGGCGCCGGGCTCgggg                                  >  1:460610/1‑136 (MQ=255)
                                                                                                         gACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGccc                         <  1:556808/136‑1 (MQ=255)
                                                                                                          aCTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCa                        <  2:22227/136‑1 (MQ=255)
                                                                                                          aCTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCa                        >  2:34130/1‑136 (MQ=255)
                                                                                                                  gACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtg                >  2:494862/1‑136 (MQ=255)
                                                                                                                   cccACCCCGGGGGCCAGGGCGGGTTCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  2:95410/135‑1 (MQ=255)
                                                                                                                   cccACCCCGGCGGCCTGGGCGCGTTGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  1:573671/135‑1 (MQ=255)
                                                                                                                   cccAACCCGGCGCCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  2:196926/135‑1 (MQ=255)
                                                                                                                   accaccCCGGCGGCCTGGGCGCTATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  1:234033/136‑1 (MQ=255)
                                                                                                                   aCCAACCGGGCGGCCAGGGCGGGTGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  1:331313/136‑1 (MQ=255)
                                                                                                                   aCCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  2:365747/136‑1 (MQ=255)
                                                                                                                   aCCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  2:497830/136‑1 (MQ=255)
                                                                                                                     cAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAg             <  2:549347/136‑1 (MQ=255)
                                                                                                                                ccAGGGCGCCATAGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATAAGCCCTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGaa  >  1:494351/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GCCGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGAA  >  NC_002947/4253336‑4253599

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: ‑