Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,7180GA40.0% 44.5 / 10.4 30intergenic (‑99/+42)uspF/ompNstress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (17/1);  new base A (0/12);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 1.50e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGT  >  NC_000913/1435574‑1435866
                                                                                                                                                |                                                                                                                                                    
aTTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATa                                                                                                                                                  >  1:540783/1‑149 (MQ=255)
       aTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaaaa                                                                                                                                         <  2:1906/151‑1 (MQ=255)
               tatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGGGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                 >  2:21956/1‑151 (MQ=255)
                   gATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATa                                                                                                                              >  1:28127/1‑150 (MQ=255)
                    aTAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATaa                                                                                                                             >  1:69859/1‑150 (MQ=255)
                          aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                                                      >  1:279681/1‑151 (MQ=255)
                                   gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                             >  2:552017/1‑151 (MQ=255)
                                   gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                             >  2:294234/1‑151 (MQ=255)
                                         ttCAAAAACCCCCTTTTTTTAGTTTCCGCATAAATTGTACCGCCAGATGCCTTTTTTTGTGTAATGACGGATTTCAATTTTTTATTTTAAATCAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAg                                                                                                         <  2:554954/149‑1 (MQ=255)
                                           cATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACt                                                                                                     >  1:573181/1‑151 (MQ=255)
                                              aaacCCTCTTTTTTTAATACCCGCAAAAAGTGTAACGCGAGATGCCATTTTTTTTTTAAGAACGGATTTCCTATTTTAATTTTAGAAAAGGGGAGAAAAAAAAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                  <  2:149486/151‑1 (MQ=255)
                                              aaacCCTCCTGTTTTAGTACCCGCATAAAGTGAAACGCCAGATGACTTTTTTTGTGTAATGCCGGAGTTCAATTTTTTAATTAAGATCAAGGGAGAAAAAATAAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                   <  2:500138/150‑1 (MQ=255)
                                                  cctcttTTTTTAGTACCGGTAAAATGTAAAACCCTAGGTTACTCTTTTTGTGTAAGGGCGGATTACATATTTTAAATGAAGAAAAAAGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                        <  1:133192/141‑1 (MQ=255)
                                                   ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                                                             >  2:171195/1‑151 (MQ=255)
                                                     cttGTTTTAGTATCCGAAAAAAGTGAACCGCCAGATGCCATTTTTTGTGAAATGACGGATTTCACTTTTTTAATTTAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCa                                                                                            <  2:166276/148‑1 (MQ=255)
                                                       ttttttAGTACCCGCATAAAGTGTACCGCCAGATGACTTTTTTTGTGTAATGACGGAGTTCATTTTTTTAATTTAAATCAAAGGAGAAAAAAAAAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAg                                                                                           <  2:120238/147‑1 (MQ=255)
                                                          tttATTATCCAAAAAAATTGAACCGCAAGATGACTTTTTTTTTTTAATGGGGGATTTAATTTTTTTAATTAAGACAAAAGGAGAAAAAAAAAAAAAAAAACCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                       <  1:396332/150‑1 (MQ=255)
                                                                         aaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcat                                                                       >  1:344450/1‑151 (MQ=255)
                                                                                 cGCAAGAGGACTTTTTTTGGGAAATGAGGGATTTCACTTTTTTAATTAAGAAAAAGGGAGAAAAAAAAAGAAAAAAACCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                                     <  1:341845/145‑1 (MQ=255)
                                                                                   ccAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                 >  1:485077/1‑147 (MQ=255)
                                                                                       aTGACCCTTTTTGTGTAATGACGGCGTTCACATTTTTCATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCCCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTCAGGCGACAACATCATCGGTAGAGATg                                                            >  1:23555/1‑148 (MQ=255)
                                                                                            acTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                    >  2:552945/1‑151 (MQ=255)
                                                                                                    tgtAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                               >  2:223326/1‑148 (MQ=255)
                                                                                                                 ttCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc                               >  1:412822/1‑151 (MQ=255)
                                                                                                                       tttttAATTTAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa                          <  2:116388/150‑1 (MQ=255)
                                                                                                                            atttaaGAACAAAGGAAAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAg                      <  1:436788/144‑1 (MQ=255)
                                                                                                                            aaTTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                    >  1:267397/1‑151 (MQ=255)
                                                                                                                               ttAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                  >  2:118810/1‑150 (MQ=255)
                                                                                                                                 aGATAAAGGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGa                 <  1:456000/149‑1 (MQ=255)
                                                                                                                                      aaaGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtata          <  2:567867/151‑1 (MQ=255)
                                                                                                                                               agaaTCAGCACAAAAAGACCGCTATACCAGCGGGTAGGCAGATTTTGAACTGATAAACTCGACCTAACGCGACAAGCTCATCGGTAGGGATGCCATTGGCAGCGGAGAAGCTGTCCACTTCATCCAACAGGTTGATTTtat           >  1:133349/1‑141 (MQ=255)
                                                                                                                                               agaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGt  >  2:378162/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGT  >  NC_000913/1435574‑1435866

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