Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,103,7160AC44.8% ‑0.8 / 11.5 30L186R (CTG→CGG) PP_1878hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/3);  new base C (4/9);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 9.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  NC_002947/2103606‑2103844
                                                                                                              |                                                                                                                                
gTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCAGGCCGGGCGGGccagt                                                                                                         >  1:660326/1‑133 (MQ=255)
         cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCCCCCCGGGCCGGGCGGGCCGGTTTGGAgggg                                                                                                >  2:669693/1‑135 (MQ=255)
         cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGTGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGc                                                                                                >  2:568541/1‑136 (MQ=255)
         cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGCCGTTTTGCAgggg                                                                                                >  2:613362/1‑135 (MQ=255)
         cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGGCGTTTTGGAGGGc                                                                                                >  1:538916/1‑136 (MQ=255)
              gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGCCGGTTTGGAggggcgggt                                                                                           >  1:239550/1‑130 (MQ=255)
              gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTTGAGGGCCAGtt                                                                                           >  1:263786/1‑136 (MQ=255)
              gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCGGGGCGTTTTTAAGGGccgggt                                                                                           >  1:268316/1‑132 (MQ=255)
              gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCGGGCCGGGCTGGGCGTTTTTGAGGGccgggt                                                                                           >  1:588674/1‑132 (MQ=255)
                      cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGCCTTTTTTCAGGGCCGGTTGGCCGAgg                                                                                   >  2:198207/1‑136 (MQ=255)
                             ttCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGTGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTga                                                                            <  1:568541/136‑1 (MQ=255)
                                tcAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga                                                                         >  2:401342/1‑136 (MQ=255)
                                          cccgccgctCATCATCCCCCTTCGGCGCCGGCGCCTTGCGCAAATACTCCCCGCCCCGCGCCAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                               <  2:57761/130‑1 (MQ=255)
                                           cAGCTCCTCTTCAACCCCCTGCGGCGCCGGCGGCTTGCCCAAATCCTCCCCGCCGCGCTCGACGGCCCGCCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgc                                                              <  1:613362/136‑1 (MQ=255)
                                               tgctCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGGCGGGCGGGCCGTTTTTCAGGGGCTGTTGGGCGCGGTTGTTGAAGAACTCAATATgaagca                                                          >  1:746926/1‑136 (MQ=255)
                                                 ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCTGGCCTTTTTTCAGGGCCTGTTAGCCGATGTTGATGAACAACACAATAAtcaacagc                                                        >  1:484777/1‑136 (MQ=255)
                                                 ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCTGGCCTTTTTTCAGGGCCGGTTGGCCGAGGTTGTTTAAGAACTCTATATtaagcagc                                                        >  1:377065/1‑136 (MQ=255)
                                                 ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCCCCTGGCCGGGCGGGCCTTTTTGCAGGGCCGGTTGGGCGAGGTTGTTGTAGAACTCTATATggagcagc                                                        >  1:141813/1‑136 (MQ=255)
                                                              ttCGGCTCCGGCTGCTTGCGAAAATACTGCCCGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                           <  2:405374/136‑1 (MQ=255)
                                                                 ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGGCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCggg                                        >  2:568142/1‑135 (MQ=255)
                                                                                     aTACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGGAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCATGACCCAATGCTTg                    >  1:90496/1‑136 (MQ=255)
                                                                                            ccgccgcGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa             <  2:560631/135‑1 (MQ=255)
                                                                                              gccgcGCGCGCCCGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACa           <  1:669693/136‑1 (MQ=255)
                                                                                               ccgcGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAg          <  2:431498/136‑1 (MQ=255)
                                                                                                       gacggCCCGGCGCGCGGGTCCTGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGAGGGCaag  <  2:431872/133‑1 (MQ=255)
                                                                                                       gAAGGCCAGGCGCGGCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  1:198207/136‑1 (MQ=255)
                                                                                                       cccggCCGGGCGCGCCGGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  1:275325/133‑1 (MQ=255)
                                                                                                       cccggCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  1:186296/133‑1 (MQ=255)
                                                                                                       cacggCCCGGCGGCCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  2:141813/133‑1 (MQ=255)
                                                                                                       cacggCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  1:534629/133‑1 (MQ=255)
                                                                                                              |                                                                                                                                
GTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  NC_002947/2103606‑2103844

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑