Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG52.9% 34.1 / 51.8 34T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (2/14);  new base G (18/0);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 8.62e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-06

GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCG  >  NC_000913/1356055‑1356343
                                                                                                                                                   |                                                                                                                                               
gTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑ATTTGCGGGG‑tt                                                                                                                                             >  1:622883/1‑150 (MQ=255)
          gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑GTAAATGccc                                                                                                                                     >  1:428019/1‑147 (MQ=255)
           cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATa                                                                                                                                 <  1:106562/151‑1 (MQ=255)
               cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCg                                                                                                                             <  1:472460/151‑1 (MQ=255)
                          gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTTAATGCCCGAAAGGCGGCGTTctttcc                                                                                                                  >  1:272594/1‑146 (MQ=255)
                           tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGTG‑TTTAATACCAGATAGGGGACGTTCATccc                                                                                                                  >  1:360684/1‑150 (MQ=255)
                                tAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑GTTAATGCCCGAAAGGGGGCGTTCATtcc                                                                                                                  >  2:584527/1‑145 (MQ=255)
                                       cgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGtt                                                                                                     <  2:5918/151‑1 (MQ=255)
                                        gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGGTTTAAATGCC‑GATAGGCGGCGGTTATTCCAgg                                                                                                               >  1:330948/1‑140 (MQ=255)
                                                       ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCgg                                                                                       <  2:110455/149‑1 (MQ=255)
                                                       ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                     <  2:510702/151‑1 (MQ=255)
                                                           aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑TTTAATGCCAGATAAGCGACGTTTATTCCCGGACGCAGCGGTTAAAGCAAGCGCggg                                                                                      >  1:282402/1‑146 (MQ=255)
                                                            ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTTAAAGCCAGATAAGCGGCGGTTCTCCCCGGGCGCAGCGTTTAAAGCCAACGCggg                                                                                      >  1:644599/1‑145 (MQ=255)
                                                             taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGG‑TTAAAT‑CCAAAAAGGGGGGGGTCATCCCAAGGAGCAGCGGTAAAAGACAACGCgggg                                                                                     >  2:652783/1‑144 (MQ=255)
                                                                 tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑TTAAATGCCAGATAGGCGACGGTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGc                                                                             >  1:288487/1‑149 (MQ=255)
                                                                 gggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaa                                                                           <  2:424634/150‑1 (MQ=255)
                                                                       gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑TTAAAAGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGGTAAACGCAAGCGCGGGTCGTCACGCAAAAAgg                                                                      >  1:425183/1‑150 (MQ=255)
                                                                                  tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGGGGGG‑TTAAATACCCAATAGGCGGCACTCTTCCCACGGAGCAGCGGTAAACGGCAGGGCGGGGGGGCCCGCAAAAAgggaagacg                                                               >  2:664005/1‑139 (MQ=255)
                                                                                   ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGctggctgg                                                          <  2:196333/150‑1 (MQ=255)
                                                                                     ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑TTAAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCAAAATGGATAGctggctgggcg                                                       >  1:625890/1‑148 (MQ=255)
                                                                                     ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTTAAAGCCCGATAGGCGACGTTCATCCCAGGGCGCAGCGGTAAAAGCAAGAGCGGGGCGTCACGCCAAAAGGATAGatgggtgg                                                          >  1:57275/1‑148 (MQ=255)
                                                                                      tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑TTAAATTCCAGATAAGCCACGGTCATCCCCGGCCCCAGCGTTAAACGCAAGCGCGGGGGGTGACGCAAAAAGGATAACtggctggccg                                                       >  1:650982/1‑148 (MQ=255)
                                                                                        ttaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTAAATGCCAGAAAAGCGACGGTCATCCCCGGACGCCGCGTTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGATAGCtggctggccg                                                       >  2:144759/1‑146 (MQ=255)
                                                                                             cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTAAATGCCCGGAAGGGCGCGTTCCTCCCACGGCGCCGGGGTTAAAGCAAGCGCGGGGCGTCACCCAAAAAGGGCAAGCGGCTGGCCGcagg                                                   >  2:602454/1‑147 (MQ=255)
                                                                                                    gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGaa                                        <  2:470326/151‑1 (MQ=255)
                                                                                                            aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑GTAAATGCCCGATAGGGGACGGTTATTCCAGGGCGCAGCGGTTAACGCAAACGCGGGGGGGTCCGCAAAAAGGAAAGGCGGGTGGGAGCCGGGCCAGCCCGAAAGTCGCa                                 >  2:243177/1‑150 (MQ=255)
                                                                                                            aTGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTTGCGGGG‑GTAAATGCCCGATAAGCGACGTTTATCCCCGGACGCAGGGGTTAACGCAAACCCGGGGGGTCACGCCAAAAGGGGAGGCGGCTGGGAGGAGGGCCGGGCAGAAAGgcgc                                  >  2:692411/1‑149 (MQ=255)
                                                                                                                cGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGG‑TTTAAAGCCAGATAGGCGGCGTTCATCCCCGGGCGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCCAAAAGGATAGgtgggtgggcg                                                       >  1:517529/1‑121 (MQ=255)
                                                                                                                cGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagcag                                                    <  2:517464/127‑1 (MQ=255)
                                                                                                                   aCAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCgg                         <  2:317390/151‑1 (MQ=255)
                                                                                                                       cGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAg                     <  2:644639/151‑1 (MQ=255)
                                                                                                                         ggATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGa                   <  1:758466/151‑1 (MQ=255)
                                                                                                                          gATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAg                  <  2:733631/151‑1 (MQ=255)
                                                                                                                                            ttGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCg  <  1:201024/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCG  >  NC_000913/1356055‑1356343

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: ‑