Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,308,7490GC38.1% 29.5 / 13.6 21intergenic (‑104/‑39)clsA/yciYcardiolipin synthase 1/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (13/0);  new base C (0/8);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 4.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGAAC  >  NC_000913/1308688‑1308804
                                                             |                                                       
cATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATaa                                                   >  1:486681/1‑68 (MQ=255)
   atGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATgg                                                >  2:637590/1‑67 (MQ=255)
       aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                            >  1:220667/1‑68 (MQ=255)
       aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                            >  2:66619/1‑68 (MQ=255)
          caccTGTCAAATAAACCCGCTGGAAAAGCCCTGGTCTGCGGCGGGGGTGGGCCCACAATGACCATTCt                                         <  2:486681/65‑1 (MQ=255)
          cATCTGCCAAATAAAACCGCTTGAAACGCACTGGTCTGCGGCGGGGGTGGGCCCACAATGACCATTCt                                         <  2:222648/68‑1 (MQ=255)
            tCTGTCAAATAAAACCGCGTTTAACTCACTGGGCCGCGTCGGGGGTGGGCCCACAATGACCATTCTGt                                       <  2:1131314/68‑1 (MQ=255)
            tCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGt                                       >  1:323889/1‑68 (MQ=255)
             cTGGCAATTATATCCGCGGTAACCGCCCTGGTCTGCGGCGGGGGTGGGCCCATAATGACCATTCTGtt                                      <  2:1140598/68‑1 (MQ=255)
                    aTAAAATCGCGTTAAATGCACTGGTCAGCGTCGGGGGTGGGCGTATAATGACCATTCTGTTATTGCAt                               <  1:340505/68‑1 (MQ=255)
                    aTAAAAACGCGTGAAAAGCACTGGTCAGCGTCGGGGGTGGGCGTATAATGACCATTCTGTTATTGCAt                               <  1:88770/68‑1 (MQ=255)
                       aaaaCGCGGGAAAAGCACTGGGCTGCGGCGGGGGGGGCCCTATAATGACCATTCTGTTATTGCATaga                            <  1:746881/68‑1 (MQ=255)
                           cgcgGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGagtag                        >  1:379100/1‑68 (MQ=255)
                           cgcgGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGagtag                        >  2:1112790/1‑68 (MQ=255)
                                    aGCGCTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGa               >  2:31348/1‑68 (MQ=255)
                                            tCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt       >  1:650245/1‑68 (MQ=255)
                                            tCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt       >  1:75008/1‑68 (MQ=255)
                                               gCGTCGGGGGTGGGCGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGa    <  1:97397/68‑1 (MQ=255)
                                                 gTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGAAc  >  2:522094/1‑68 (MQ=255)
                                                 gTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGAAc  >  2:713617/1‑68 (MQ=255)
                                                 gTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCAAAGAGTAGTTAACATGAAGCGGAGTAGAAc  >  1:167062/1‑68 (MQ=255)
                                                             ggtATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt       <  1:970529/51‑1 (MQ=255)
                                                             ggtATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt       >  2:970525/1‑51 (MQ=255)
                                                             |                                                       
CATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGAAC  >  NC_000913/1308688‑1308804

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: