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

GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCG  >  NC_000913/1356055‑1356343
                                                                                                                                                   |                                                                                                                                               
gTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGATTTGCGGGGtt                                                                                                                                             >  1:622883/1‑150 (MQ=255)
          gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAATGccc                                                                                                                                     >  1:428019/1‑147 (MQ=255)
           cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATa                                                                                                                                 <  1:106562/151‑1 (MQ=255)
               cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCg                                                                                                                             <  1:472460/151‑1 (MQ=255)
                          gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAATGCCCGAAAGGCGGCGTTctttcc                                                                                                                  >  1:272594/1‑146 (MQ=255)
                           tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGTTTAATACCAGATAGGGGACGTTCATccc                                                                                                                  >  1:360684/1‑150 (MQ=255)
                                tAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGCCCGAAAGGGGGCGTTCATtcc                                                                                                                  >  2:584527/1‑145 (MQ=255)
                                       cgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGtt                                                                                                     <  2:5918/151‑1 (MQ=255)
                                        gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAAATGCCGATAGGCGGCGGTTATTCCAgg                                                                                                               >  1:330948/1‑140 (MQ=255)
                                                       ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCgg                                                                                       <  2:110455/149‑1 (MQ=255)
                                                       ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                     <  2:510702/151‑1 (MQ=255)
                                                           aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAATGCCAGATAAGCGACGTTTATTCCCGGACGCAGCGGTTAAAGCAAGCGCggg                                                                                      >  1:282402/1‑146 (MQ=255)
                                                            ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTTAAAGCCAGATAAGCGGCGGTTCTCCCCGGGCGCAGCGTTTAAAGCCAACGCggg                                                                                      >  1:644599/1‑145 (MQ=255)
                                                             taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAATCCAAAAAGGGGGGGGTCATCCCAAGGAGCAGCGGTAAAAGACAACGCgggg                                                                                     >  2:652783/1‑144 (MQ=255)
                                                                 tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAGGCGACGGTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGc                                                                             >  1:288487/1‑149 (MQ=255)
                                                                 gggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaa                                                                           <  2:424634/150‑1 (MQ=255)
                                                                       gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGGTAAACGCAAGCGCGGGTCGTCACGCAAAAAgg                                                                      >  1:425183/1‑150 (MQ=255)
                                                                                  tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTAAATACCCAATAGGCGGCACTCTTCCCACGGAGCAGCGGTAAACGGCAGGGCGGGGGGGCCCGCAAAAAgggaagacg                                                               >  2:664005/1‑139 (MQ=255)
                                                                                   ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGctggctgg                                                          <  2:196333/150‑1 (MQ=255)
                                                                                     ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCAAAATGGATAGctggctgggcg                                                       >  1:625890/1‑148 (MQ=255)
                                                                                     ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGCCCGATAGGCGACGTTCATCCCAGGGCGCAGCGGTAAAAGCAAGAGCGGGGCGTCACGCCAAAAGGATAGatgggtgg                                                          >  1:57275/1‑148 (MQ=255)
                                                                                      tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATTCCAGATAAGCCACGGTCATCCCCGGCCCCAGCGTTAAACGCAAGCGCGGGGGGTGACGCAAAAAGGATAACtggctggccg                                                       >  1:650982/1‑148 (MQ=255)
                                                                                        ttaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAGAAAAGCGACGGTCATCCCCGGACGCCGCGTTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGATAGCtggctggccg                                                       >  2:144759/1‑146 (MQ=255)
                                                                                             cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCCGGAAGGGCGCGTTCCTCCCACGGCGCCGGGGTTAAAGCAAGCGCGGGGCGTCACCCAAAAAGGGCAAGCGGCTGGCCGcagg                                                   >  2:602454/1‑147 (MQ=255)
                                                                                                    gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGaa                                        <  2:470326/151‑1 (MQ=255)
                                                                                                            aTGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAATGCCCGATAGGGGACGGTTATTCCAGGGCGCAGCGGTTAACGCAAACGCGGGGGGGTCCGCAAAAAGGAAAGGCGGGTGGGAGCCGGGCCAGCCCGAAAGTCGCa                                 >  2:243177/1‑150 (MQ=255)
                                                                                                            aTGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTAAATGCCCGATAAGCGACGTTTATCCCCGGACGCAGGGGTTAACGCAAACCCGGGGGGTCACGCCAAAAGGGGAGGCGGCTGGGAGGAGGGCCGGGCAGAAAGgcgc                                  >  2:692411/1‑149 (MQ=255)
                                                                                                                cGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAAAGCCAGATAGGCGGCGTTCATCCCCGGGCGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCCAAAAGGATAGgtgggtgggcg                                                       >  1:517529/1‑121 (MQ=255)
                                                                                                                cGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagcag                                                    <  2:517464/127‑1 (MQ=255)
                                                                                                                   aCAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCgg                         <  2:317390/151‑1 (MQ=255)
                                                                                                                       cGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAg                     <  2:644639/151‑1 (MQ=255)
                                                                                                                         ggATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGa                   <  1:758466/151‑1 (MQ=255)
                                                                                                                          gATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAg                  <  2:733631/151‑1 (MQ=255)
                                                                                                                                            ttGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCg  <  1:201024/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCG  >  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: