Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,020,3160AC35.0% 24.4 / 11.5 20T31P (ACC→CCC) ubiBubiquinone biosynthesis protein UbiB
Reads supporting (aligned to +/- strand):  ref base A (13/0);  new base C (0/7);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 1.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTT  >  NC_000913/4020249‑4020383
                                                                   |                                                                   
ccTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCa                                                                     <  1:1108603/68‑1 (MQ=255)
     atTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTg                                                                >  1:1899275/1‑68 (MQ=255)
       ttCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGcc                                                              >  1:1500393/1‑68 (MQ=255)
        tcatcaCTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCg                                                             >  1:42823/1‑68 (MQ=255)
               tCGCATTTTTTTAAGCTACGGACTTGATGAACTGATCCCAAAAATGCGTATCCCCCTGCCGCTACGGc                                                      <  1:916301/68‑1 (MQ=255)
               tCGCACTTTTTTAAGCTACGGACTTGATGAACGGATCCCCAAAATGCGTATCCCCCTGCCGCTACGGc                                                      <  1:1175502/68‑1 (MQ=255)
                   ccTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATg                                                  <  1:1042784/67‑1 (MQ=255)
                          taggCTACGGACTTGATGAACTGATCCCAAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATAc                                           <  1:802314/65‑1 (MQ=255)
                             gCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCa                                        >  1:855944/1‑68 (MQ=255)
                             gCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTAAGGCTATACTCa                                        >  1:273725/1‑68 (MQ=255)
                             gCTACGGACTTGATGAACTGATCCCCAAAATGCGTACCCCCCTGCCGCTACGGCTATGGCGATACTCa                                        <  1:3444943/68‑1 (MQ=255)
                             gCTACGGACTTGATGAACGGATCCCCAAAAGGCGTACCACCCTGCCGCTACGGCTATGGCGATACTCa                                        <  1:2917915/68‑1 (MQ=255)
                             gCTACGGACTTGATGAACGGACCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCa                                        <  1:2565233/68‑1 (MQ=255)
                                 cGGACTTGATGAACTTACCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCattat                                    <  1:94998/68‑1 (MQ=255)
                                 cGGACTTGATGAACTGATCCCAAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCattat                                    <  1:1498595/68‑1 (MQ=255)
                                    aCTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCt                                 >  1:1473070/1‑68 (MQ=255)
                                     cTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTg                                >  1:2745105/1‑68 (MQ=255)
                                       tgatgaACTAACCCCCAAAATGCTTTCCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGa                              <  1:2835724/68‑1 (MQ=255)
                                          tgaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGc                           >  1:800179/1‑68 (MQ=255)
                                            aaCTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCa                         >  1:3291415/1‑68 (MQ=255)
                                            aaCTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCa                         >  1:1447813/1‑68 (MQ=255)
                                                     ccAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCAt                >  1:1642493/1‑68 (MQ=255)
                                                       aaaaTGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATaa              >  1:1892577/1‑68 (MQ=255)
                                                               tATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACGAAc      <  1:582200/68‑1 (MQ=255)
                                                               tATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAAc      <  1:1657081/68‑1 (MQ=255)
                                                               tATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAAc      <  1:155428/68‑1 (MQ=255)
                                                               tATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAAc      <  1:1487317/68‑1 (MQ=255)
                                                               tATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAAc      <  1:695136/68‑1 (MQ=255)
                                                                 tCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGCTGCCAAATCGGCATAAAGACAACCtt    >  1:2976107/1‑68 (MQ=255)
                                                                   aCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGGCGAAAGACAAACtttt  >  1:1310988/1‑68 (MQ=255)
                                                                   aCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACtttt  >  1:1738087/1‑68 (MQ=255)
                                                                   |                                                                   
CCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTT  >  NC_000913/4020249‑4020383

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