Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273711,0810AC53.3% 9.3 / 13.9 15T72P (ACC→CCC) BW25113_RS03560lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (0/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAACCAGAGGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAA  >  NZ_CP009273/710992‑711165
                                                                                           |                                                                                    
aaaCCAGAGGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGataata                                                                                      >  1:469101/1‑90 (MQ=255)
      gagGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACccgcc                                                                                >  1:276584/1‑90 (MQ=255)
              aggggTAGCTTTGCGGAAAATATTGCCCTTTCCGGGTGAAAAAAAAATTGCCGTGGTGATAACTTCCCTGATAATCCCCGCCGCTGGCAg                                                                        <  1:451372/88‑1 (MQ=255)
                     gCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGGTGTCAg                                                                 >  1:437221/1‑90 (MQ=255)
                                  tatTCCCCCTTTCCGGGGAAAGCAAAATTCGCGGGGTAACAACTTCCCTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGaaaa                                                    <  2:27452/90‑1 (MQ=255)
                                    ttGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAg                                                  >  1:723107/1‑90 (MQ=255)
                                     tGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGa                                                 >  2:60153/1‑90 (MQ=255)
                                          cTTTCCGGGTGTAAGAAAAATTGCGGGGGTAAAAACTTACCTAAAAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAactt                                             <  2:280574/90‑1 (MQ=255)
                                           tttCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAacttac                                           >  2:374120/1‑90 (MQ=255)
                                              ccGGTGTAAGAAAAATTGCCGTGGTGAAAACTTCCCTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAActtactta                                        <  2:276584/90‑1 (MQ=255)
                                               cGGTGAAAGACAAATTGCGGGGGGAAAAACTAACCTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTAt                                       <  2:519294/90‑1 (MQ=255)
                                                 ggTGTAAGACAACTTGCCGTGGTAATAACTTCCCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATc                                      <  1:665738/90‑1 (MQ=255)
                                                  gtgaAAGAAAAATGGCCGTGGTGATAACTTACCTAAAAACCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCt                                     <  1:438085/90‑1 (MQ=255)
                                                      aaGACAAATTGCCGTGGTGATAACTTACCTGATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCt                                     <  1:55348/86‑1 (MQ=255)
                                                      aaGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCt                                     >  2:55348/1‑86 (MQ=255)
                                                          aaaaTTGGCCGTGGTAATAATTTCCCTGATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCataat                              <  1:558437/89‑1 (MQ=255)
                                                                                      ataataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTaa  >  1:133401/1‑90 (MQ=255)
                                                                                           |                                                                                    
AAACCAGAGGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAA  >  NZ_CP009273/710992‑711165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: