Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG53.6% 26.0 / 37.0 29T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (4/9);  new base G (15/0);  total (20/9)
Fisher's exact test for biased strand distribution p-value = 1.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-04

GTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGC‑ATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356058‑1356350
                                                                                                                                               |                                                                                                                                                      
gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTaaa                                                                                                                                                  >  2:478621/1‑150 (MQ=255)
gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTaaa                                                                                                                                                  <  2:273305/150‑1 (MQ=255)
   aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGTCAGCGGGATTATTTAGCGGCGGTATGGCGGTGTTTAATg                                                                                                                                                >  2:643926/1‑149 (MQ=255)
       gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCCGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGATTGGCGGCGTTTAATGc                                                                                                                                               >  1:289564/1‑146 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                     <  2:441613/151‑1 (MQ=255)
                                       cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAAAAAAGGGACGTTTATCCCCGGGCACAGCGttta                                                                                                           >  1:180398/1‑148 (MQ=255)
                                          ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAAATCCAGATAGGGGGCCTTTCTCCCACGACGCAGCgctta                                                                                                           >  1:592736/1‑143 (MQ=255)
                                            cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCCGAAAAGCGACGTTCATCCCCGGACGCAGCGGTAAAAGCa                                                                                                      >  1:130006/1‑150 (MQ=255)
                                              gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAAAGCCAGATAAGAGACGTTCATCCCAAGACGCCGCGGTAAACGCaa                                                                                                     >  2:569480/1‑149 (MQ=255)
                                                    ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                             <  2:619196/151‑1 (MQ=255)
                                                      ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTAAATGCCAGAAAAGCGAAGTTCCTCCCAAGGCGCCGCGGTAAAACCAAGCGCcggg                                                                                             >  1:247796/1‑148 (MQ=255)
                                                      ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtcg                                                                                           <  1:629352/151‑1 (MQ=255)
                                                       taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGTTAAATGCCACAAAAGGGACGTTCATCCCCGGACGCAGCGTTAAAAGCAAACGCggg                                                                                              >  1:644305/1‑147 (MQ=255)
                                                       taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTATATGCCAAAAAAGCGGCGTTCATTCCCGGGAGCAGCGGTAAACGCAAGCGCGGggcg                                                                                           >  1:701551/1‑150 (MQ=255)
                                                         ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCa                                                                                        <  1:207148/151‑1 (MQ=255)
                                                                gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAAAGCCAAATAAGCGAAGGTCATCCCCGGAAGCCGCGTTAAAAGCAAACGCGGGGCGTCACGcca                                                                                   >  2:200344/1‑147 (MQ=255)
                                                                gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAGATAGGGGACGGTTATCCCACGGCGCAGCGGTAAAAGCAAGCGCGGGGGGGCACGCaaa                                                                                  >  1:229383/1‑150 (MQ=255)
                                                                        cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGCCGTGGTAAATGTCAGCGAGGGGACGCTCATTCCACGGCGCACCGTTAAAGACTAACGCGGGGCGGAACACAAAATGGaa                                                                            >  1:650768/1‑147 (MQ=255)
                                                                           gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctgg                                                                      <  1:483639/151‑1 (MQ=255)
                                                                                     ttaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGATAAGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACCCAAAATGGATAGCTGGCTGgcag                                                               >  2:85620/1‑148 (MQ=255)
                                                                                           gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggccagg                                                      <  2:783205/151‑1 (MQ=255)
                                                                                                    aGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGCCAGAAAGGCGGCGTTCATCCCAGGGGGCAGCGGTAAACGCAAGGGCGGGTGGGCACCCAAAAAGGATAGGTGGCTGGGCGCagggcaggcccaaa                                                >  2:772012/1‑144 (MQ=255)
                                                                                                           ggCGGACAGCGGGATTATTTAGCGGCGGCTTGGCGGGGGTTAAAGCCCGACAGGGGACGTTTATCCCCGGACGCAGCGTTAAACAGAAGCGCGGGGCGCCACGCAAAATGGATAGGGGGCTGGCAGCagggccggccccaaccgcgca                                         >  2:466977/1‑137 (MQ=255)
                                                                                                                  aGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAACCAGATAAGCGACGTTTATCCCCGGACGCAGCGGTTAAAGCAAGCGGGGGGCGGCACGCAAAAAGGGTAGCTTGGCGGGAGCAGGGCAAGCCAGAACGTCGG‑ATTcccc                                  >  1:413124/1‑147 (MQ=255)
                                                                                                                  aGCGGGATTATTTAGCGGCGGGTTGGCGGGGTTAAATGGCAGATAGGGGGCGTTTATTCCAGGAAGGAGCGGTAAAAGCAAGAGGGGGTCGTCACGGAAAAAGGATAGGTGGGTGGCAGCAGGGCCGGGAAGAAAGTCGC‑ATTCCCCgg                                >  2:258754/1‑149 (MQ=255)
                                                                                                                        attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGC‑ATTCCCCGGTCAGGAGt                        <  1:359405/151‑1 (MQ=255)
                                                                                                                          tattTAGCGGCGGTTTTGCGGGGTTAAATTCCAGATAGGCGACGTTCATCCCCGGGCGCAACGTTAAACGCAAACGCCGGGCGTCACGCAAAATGGATTACCTGCCGGCAGCACGCCAAGCCCGAACGCCGCAATTCCCCGTTccgagg                         >  1:758240/1‑144 (MQ=255)
                                                                                                                                 cggcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGC‑ATTCCCCGGTCAGGAGTTTCGACAGa               <  2:99125/151‑1 (MQ=255)
                                                                                                                                             ggTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGC‑ATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgcc    >  2:231079/1‑150 (MQ=255)
                                                                                                                                              gtgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGC‑ATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc  >  1:636555/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
GTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGC‑ATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356058‑1356350

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