Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,041,3860GC26.1% 29.3 / 12.0 23G257A (GGA→GCA) appAphosphoanhydride phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base C (0/6);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATT  >  NC_000913/1041326‑1041443
                                                            |                                                         
cAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATc                                                    <  2:228441/68‑1 (MQ=255)
           gagaTATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTcaca                                         <  1:326826/68‑1 (MQ=255)
           gagaTATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTcaca                                         <  1:866921/68‑1 (MQ=255)
              atatTTCTCCTGCAACAGGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCa                                      <  1:510465/68‑1 (MQ=255)
                atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGTGTGGGCAAGGATCACCGATTCACACCAGt                                    <  2:969063/68‑1 (MQ=255)
                atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGt                                    <  1:558954/68‑1 (MQ=255)
                atTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGCAAGGATCACCGATTCACACCAGt                                    <  2:882305/68‑1 (MQ=255)
                   tctcCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGa                                 <  2:698726/68‑1 (MQ=255)
                    ctcCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGaa                                >  2:593302/1‑68 (MQ=255)
                    ctcCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATAACCGATTCACACCAGTGGaa                                >  2:591957/1‑68 (MQ=255)
                     tcCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCa                                      <  1:361330/61‑1 (MQ=255)
                     tcCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCa                                      >  2:361330/1‑61 (MQ=255)
                      ccTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAaca                              >  2:192730/1‑68 (MQ=255)
                      ccTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAaca                              >  1:918414/1‑68 (MQ=255)
                           aacaaGCACAGGGAATGCCGGAGCCGGGGCGGGCAAGGATCACCGATTCACACCAGTGGAACACCTTg                         <  2:270033/68‑1 (MQ=255)
                           aacaaGCACAGGGAATGCCGGAGACGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTg                         >  2:848961/1‑68 (MQ=255)
                                gCACAGGGAATGCCGGAGCCGGGGTGGGGACGGATCACCGATTCACACCAGTGGAACACCTTGCTAAg                    <  1:602998/68‑1 (MQ=255)
                                       gAATGCCGGAGCCGGGGTGGGGTAGGATCACCGATGCACACAAGTGGAACACCTGGCTAGGTTTGCAt             >  2:377128/1‑68 (MQ=255)
                                       gAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACCCCAGTGGAACACCTTGCTAAGTTTGCAt             >  1:419067/1‑68 (MQ=255)
                                          tGCCGGAGCCGTGGTTGCGTAGGATCACCGATTCACACCAGTGGATCACCTTGCTACGGTTGCATAAc          <  1:737874/68‑1 (MQ=255)
                                                 gCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACAGCTTGCTAAGTTTGCATAACGCGCAAt   >  1:795804/1‑68 (MQ=255)
                                                 gCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAAt   <  2:297251/68‑1 (MQ=255)
                                                  ccGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAAtt  >  1:857576/1‑68 (MQ=255)
                                                            |                                                         
CAATGCTGACGGAGATATTTCTCCTGCAACAAGCACAGGGAATGCCGGAGCCGGGGTGGGGAAGGATCACCGATTCACACCAGTGGAACACCTTGCTAAGTTTGCATAACGCGCAATT  >  NC_000913/1041326‑1041443

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: