Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,041,3860GC45.5% 14.9 / 12.7 22G257A (GGA→GCA) appAphosphoanhydride phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base C (0/10);  total (7/15)
Fisher's exact test for biased strand distribution p-value = 5.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATTTTATTTGCTA  >  NC_000913/1041330‑1041453
                                                        |                                                                   
gCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCg                                                          >  1:488565/1‑68 (MQ=255)
 cTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGa                                                         <  1:348254/68‑1 (MQ=255)
  tGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGGGGGGAAGGATCACCGAt                                                        <  1:369391/68‑1 (MQ=255)
  cgACGGAGATATTTCTCCTGCAACATGCTCAGGGAATGCCGGAGCCGTGGTGGGCAATGATCACCGAt                                                        <  1:306155/67‑1 (MQ=255)
      ggAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTcac                                                    <  1:270796/68‑1 (MQ=255)
       gagaTATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTcaca                                                   <  1:645529/68‑1 (MQ=255)
        agaTATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGAt                                                        >  2:567304/1‑62 (MQ=255)
        agaTATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGAt                                                        <  1:567303/62‑1 (MQ=255)
          atatTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCa                                                <  2:613179/68‑1 (MQ=255)
            atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGTCAAGGATCACCGATTCACACCAGt                                              <  2:74058/68‑1 (MQ=255)
            atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGt                                              >  2:562269/1‑68 (MQ=255)
                ctcCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGag                                          >  2:698578/1‑67 (MQ=255)
                   cTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGTGGAacac                                       <  1:403798/68‑1 (MQ=255)
                     gCAACAAGCACTGGGATTGCCGGTGCCGGGGTGGGCAAGGATCAACGATTCACACCAGTGGAACACCt                                     <  1:230086/68‑1 (MQ=255)
                     gCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGTGGAACACCt                                     <  1:511488/68‑1 (MQ=255)
                            gCACTGGTAATGCCGTAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAg                              <  2:488565/68‑1 (MQ=255)
                            gCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAg                              <  1:405994/68‑1 (MQ=255)
                             cacaGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGt                             >  2:542855/1‑68 (MQ=255)
                                   gAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCAt                       >  2:508417/1‑68 (MQ=255)
                                     aTGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATaa                     >  2:65668/1‑68 (MQ=255)
                                        ccGTAGCCGTGTTGTGGAAGTAGCACCGATTCTCATCTGTTGAACACCTTGCTAAGTTTGCATAAcgc                  <  2:487159/68‑1 (MQ=255)
                                             gCCGGGGTGGGCTAGGAGCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAAt             <  1:144635/68‑1 (MQ=255)
                                                        gAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATTTTATTTGCTa  >  1:211034/1‑68 (MQ=255)
                                                        |                                                                   
GCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATTTTATTTGCTA  >  NC_000913/1041330‑1041453

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: