Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG54.5% 20.4 / 44.6 33T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (4/11);  new base G (18/0);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 7.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-05

AGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356063‑1356349
                                                                                                                                           |                                                                                                                                                    
aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAAAGcc                                                                                                                                            >  2:647553/1‑149 (MQ=255)
    ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATg                                                                                                                                              >  2:128421/1‑143 (MQ=255)
        tattctctctgTAAAATATAGTCCTCCTATTCAGCCGTACCATAACAGAGGGATAGAATAAGGCTCTTGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTTGGGGGTTTAATGCCAGATAAGCg                                                                                                                                   >  1:366692/11‑150 (MQ=255)
            cgcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTTTATGGCCGATAAgcggcg                                                                                                                                >  2:646399/1‑149 (MQ=255)
            cgcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGGGTAAAtaaccgaaaag                                                                                                                                     >  2:417827/1‑134 (MQ=255)
             gcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTc                                                                                                                             <  2:349022/151‑1 (MQ=255)
                cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAAGGCCAAAAAggg                                                                                                                                    >  2:492251/1‑141 (MQ=255)
                  gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATccc                                                                                                                        <  2:692921/151‑1 (MQ=255)
                                gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGAAAGGCGACGTTCCTCCCAGGGCGCAGCGtt                                                                                                           >  1:19296/1‑150 (MQ=255)
                                   gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAAc                                                                                                       <  1:225257/151‑1 (MQ=255)
                                     ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAATGGCAAATAAGGGGCGTTCATCCCAGGACGCAGCGTTtaa                                                                                                        >  1:633698/1‑148 (MQ=255)
                                               ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                           <  1:71460/151‑1 (MQ=255)
                                               ccACAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                           <  1:628924/151‑1 (MQ=255)
                                                   aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGGTTAATGCCAGATAAGCGACGTTCCTCCCAAGGCGCCGCGTTAAAAACAAGCGCGGggcg                                                                                         >  1:1287/1‑149 (MQ=255)
                                                         tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAATGGCAGAAAGGCGAGGTTCATCCCCGGGCGCAGCGGTAAAAGCAAGCGCGGGGCGTCAc                                                                                     >  1:523606/1‑147 (MQ=255)
                                                          ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaa                                                                                 <  2:12918/150‑1 (MQ=255)
                                                             ttGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGGTTAAATGCCAGAAAGGCGACGTTCATCCCAGGACGCAGCGTTAAAAGCAAACGCGGGGCGGCACGCcaaa                                                                               >  1:584206/1‑149 (MQ=255)
                                                                   cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGGGGGGTTTAAAGCCAAATAAGCGACGGTCCACCCACGACGCCGCGGTAAAAGCAAAAGCCGGGCGGCACCCaaaa                                                                               >  2:470773/1‑143 (MQ=255)
                                                                      gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTGGCGGGGGTAAATGCCAGATAGGCGACGTTCCTCCCCCGACCCAGCGTTAAAAGCCAACGCGGGGTGTCCCGCAAAAAGGATAgctgg                                                                    >  2:692549/1‑151 (MQ=255)
                                                                       cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggc                                                                   <  2:233474/151‑1 (MQ=255)
                                                                           ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggctg                                                                 <  2:795330/149‑1 (MQ=255)
                                                                            gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgca                                                              <  1:273851/151‑1 (MQ=255)
                                                                                         aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTGGCGGGGGTTAAAGCCAGAAAGGCGACGTTCATCCCAGGGCGCCGCGCTAAAAGCAAACGCGGGGGGTCACGCAAAAAAGATAGGGGGCTGggag                                                             >  2:660444/1‑139 (MQ=255)
                                                                                            gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTGAATGACAGGTAGGCGGCGGTCATCCCCGGACGCCGCGGTAAACGCCAGAGCGGGTCGTCGCGCCAAAAGGATAACCGGCGGGCAGgaggcccg                                                     >  2:23582/1‑142 (MQ=255)
                                                                                                 cgcATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTGGCGGGGGTAAAAGCCAGATAGGGGACGGTCATCCCGGGACGCCGCGGTAAACCAAACCGCGGGGCGACCCCCCAAATGGGAAAGTTGCCGGCCGCAAGGCCGGCCAGaa                                              >  2:749760/1‑146 (MQ=255)
                                                                                                         ggACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAATGGCCGCTAGGCGACGTTTATCCCAGGACGCAGCGGTTAACGCAAACGCGGGGCGGCACGCAAAATAGATAACTGGCTGGCAGCAGGCCAGGCCAGAACGCCGCATTcccc                                 >  2:749044/1‑151 (MQ=255)
                                                                                                          gACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAAAGCCAGATAAGCGACGTTCATCCCCGGACGCAGCGGTAAACGCCAGCGCGGGGCGGCACGCCAAATGGATAGGTGGCTGGCAGCAGGGCAGGGCAGAACGTCGCATTtccc                                 >  1:540041/1‑150 (MQ=255)
                                                                                                                gggAGTATTTAGCGGAGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAgg                          <  2:424902/151‑1 (MQ=255)
                                                                                                                     tattTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGATAAGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAACGCGGGGCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCCGAACGGCGCATTTCCCCGGCAGGAGgtt                     >  1:253186/1‑151 (MQ=255)
                                                                                                                              gcgg‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg            >  1:587282/1‑151 (MQ=255)
                                                                                                                                 g‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc         >  1:343570/1‑151 (MQ=255)
                                                                                                                                     tGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg       <  2:670114/150‑1 (MQ=255)
                                                                                                                                         ggTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg  >  1:235482/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
AGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356063‑1356349

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