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

GTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGT  >  NC_000913/1435570‑1435866
                                                                                                                                                     |                                                                                                                                                     
gTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGgaagaa                                                                                                                                                     <  1:399428/151‑1 (MQ=255)
              ggAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGc                                                                                                                                      >  1:574436/1‑151 (MQ=255)
                  atatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCg                                                                                                                                   >  2:22848/1‑150 (MQ=255)
                        aTAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAAc                                                                                                                            >  2:73055/1‑151 (MQ=255)
                           gggACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGc                                                                                                                         >  1:719341/1‑151 (MQ=255)
                              aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                                                      >  2:565999/1‑151 (MQ=255)
                                aaaaaTCTTTCTGTTAATAAACCCTCCTTTTTTAGTACCCGCAAAAAGTATAACCCCAGATGCCATTTTTTGTAAAAGGAGGAAGTTCCCATTTTTATTTA‑AAAACAAAGGGGAAAAAAA‑AACAAAAAAAAGCCCGCCATAACAGCgg                                                                                                                       <  2:429730/145‑1 (MQ=255)
                                        tttttttCAAAAACCCTCTTTTTTTAGTACCCGCAAAAAGTGAAACGCAAGATGACCTTTTTTGTGTAATGAGGGAGTTCAATTTTTTAATTA‑AGAAAAAAGGAGAAAAAAT‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCaggagg                                                                                                               <  1:615673/143‑1 (MQ=255)
                                        ttCTTTTAAAAAACCCTCTTTTTTTAGTACCCGAAAAAAGTGTAACGCAAGATGACATTTTTGGTAAAAGGAGGTATTACACATTTTAAATTA‑AGAACAAGGGAGAAAGAAT‑AACAAGAAAAAGCCCGCCATAACAGCGGGCaggag                                                                                                                <  2:456096/147‑1 (MQ=255)
                                             ttCAAAAACCCTCCTTTTTTAGTATCCGAAAAAAGTGTAACGCAAGATGACTTTTTTTGTGTAATGACGGATTTCACTTTTTTAATTT‑AGAAAAAAGGAGGAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                        <  1:617180/150‑1 (MQ=255)
                                              taaaaaacCCTCTTTTTTTATTACCCGAAAAAAGTGTAACGCAAAATGACTTTTTTTGTGAAAGGACGGATTCCACATTTTAAATTA‑AGAGAAAGGGAGAAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                        <  2:239504/145‑1 (MQ=255)
                                              tCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAAc                                                                                                      >  1:497294/1‑151 (MQ=255)
                                                  aaacCCTTCTGTTTTAGTATCCGAATAATGTGTAACGCAAGATGACTTTTTTGGTGTAATGAGGGATTTCACTTTTTTAATTT‑AAATCAAGGGAGAAAAAAT‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACt                                                                                                     <  2:784747/148‑1 (MQ=255)
                                                  aaacCCTCTTGTTTGAGTAACCACAGAAAGTTAAACCACAGATGACTTTTTTTGTGAAATGAGGGAGGCCACATTTTAAATTA‑AGATCAAAGG‑GAAAAAATGAAACAAAAAAAGTCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                   <  2:577358/150‑1 (MQ=255)
                                                  aaacCCTCCTTTTTTATTACCCTAAAAAAGTAAAACGCAAAATGACTTTTTTGGTGAAATGACGGAGTTCAC‑TTTTTAATTTAAGAAAAAAGGAGAAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTa                                                                                                          <  2:401517/143‑1 (MQ=255)
                                                   aaaCCTCGTGTTTGATTACCCGAAAAAATTGAAACGCCAGAGGACATTTTTTGGGAAAGGACGGAGTTCGCATTTTAAATTA‑AGAAAAAGGGAGAAAAAAT‑AAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa                                                                                                 <  2:643589/151‑1 (MQ=255)
                                                    aaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa                                                                                                 >  1:642576/1‑150 (MQ=255)
                                                      cctcttGTTTGAGTACCGGAATAAAGTGTAACGCCAGATGACATTTTTTGTGAAAGGCGGGAGTTCACATTTTTAATTT‑AGATAAAGGGAGAAAAAAA‑AAGAAGAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa                                                                                               <  1:356172/150‑1 (MQ=255)
                                                       ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGGGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                                                             >  2:467563/1‑151 (MQ=255)
                                                       ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                                                             >  2:107763/1‑151 (MQ=255)
                                                        tcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAc                                                                                              >  1:402686/1‑149 (MQ=255)
                                                            gTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                                        >  1:168376/1‑151 (MQ=255)
                                                                 aGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa                                                                                     >  1:398442/1‑149 (MQ=255)
                                                                   tATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                                 >  1:855855/1‑151 (MQ=255)
                                                                      ccGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                                                               >  2:96323/1‑150 (MQ=255)
                                                                          aaaaagggaaccGCCAAATCACTTTTTTTGTAAAAAGAGGGATTTCATATTTTAAATAA‑AAATGAAGGAAAAAAAAAA‑AAGAAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACaa                                                                             <  2:645183/137‑1 (MQ=255)
                                                                          aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                          >  1:130944/1‑151 (MQ=255)
                                                                          aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                          >  1:110600/1‑151 (MQ=255)
                                                                                   accGCCAGATGACAATTTTGTTATAATGAGGGAGTTCACTTTTTTATTTT‑AGAACAAGGGAGAAAAATT‑AAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCgg                                                                    <  2:771604/146‑1 (MQ=255)
                                                                                         agagGACATTTTTTGTGTAAGGACGGAGTTCATTTTTTTAATTA‑AGAAAAAAGGAGAAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                            <  2:551526/150‑1 (MQ=255)
                                                                                           aTGACTCTTTTTGTGAAAGGAGGGGTGTCAAATTTTAAATTA‑AAAAAAAAGGAGAAAAAAA‑AAGAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTagag                                                               <  2:123515/145‑1 (MQ=255)
                                                                                            tGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAt                                                        >  1:484924/1‑151 (MQ=255)
                                                                                              actttttttGTGAAATGAGGGATTTCAATTTTTTAATTA‑AGAAAAAGGGAGAAAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGc                                                           <  1:784722/142‑1 (MQ=255)
                                                                                                tttttttGTATAGGGGCGGAGTACACATTTTAAATTT‑AAAAAAAAGAAGAAAAAAA‑AAAAAAAAAAGTCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                     <  1:421148/148‑1 (MQ=255)
                                                                                                    ttggtatAATGAGGGATTCCAATTTTTTAATAT‑AAATAAAAGGAGAAAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                                 <  2:651500/144‑1 (MQ=255)
                                                                                                       gtgtAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                               >  2:512508/1‑149 (MQ=255)
                                                                                                                   atttCACATTTTAAATTG‑AGAAAAAAGGAAGAAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATc                                  <  2:336481/148‑1 (MQ=255)
                                                                                                                    gTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCtt                                >  1:165839/1‑151 (MQ=255)
                                                                                                                     ttCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc                               >  1:648091/1‑151 (MQ=255)
                                                                                                                         caTTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATcc                           >  2:390252/1‑151 (MQ=255)
                                                                                                                                aaTTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                    >  2:167803/1‑151 (MQ=255)
                                                                                                                                  ttt‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                  >  1:49292/1‑151 (MQ=255)
                                                                                                                                         tCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt             >  2:640786/1‑150 (MQ=255)
                                                                                                                                              ggagaaAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt        <  2:227005/150‑1 (MQ=255)
                                                                                                                                                  gaagaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACg   >  2:378158/1‑151 (MQ=255)
                                                                                                                                                   aagaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtata          >  2:94060/1‑143 (MQ=255)
                                                                                                                                                   aagaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGt  >  1:817862/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGT  >  NC_000913/1435570‑1435866

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