Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,7180GA43.2% 58.3 / 22.2 44intergenic (‑99/+42)uspF/ompNstress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (21/4);  new base A (0/19);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 6.43e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTAT  >  NC_000913/1435575‑1435857
                                                                                                                                                 |                                                                                                                                            
ttCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGc                                                                                                                                       >  2:639822/1‑151 (MQ=255)
    gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAg                                                                                                                                     <  1:531766/149‑1 (MQ=255)
         ggAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGc                                                                                                                             >  1:168538/1‑151 (MQ=255)
         ggAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGc                                                                                                                             >  1:505311/1‑151 (MQ=255)
             atatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGc                                                                                                                         <  2:462573/151‑1 (MQ=255)
             atatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGc                                                                                                                         <  1:746988/151‑1 (MQ=255)
                    tAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAAc                                                                                                                   >  2:444421/1‑150 (MQ=255)
                         aaaaGAATCGTTCTTTTAAAAAACCCTCCTTTTTTAGTACCCGAAAAAAGTTTAACGCCAGATGACTTTTTTTGTGTAAGGACGGATTTCACATTTTTAATTTAGATAAAGGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCggg                                                                                                             <  1:452004/149‑1 (MQ=255)
                         aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                                             >  1:782675/1‑151 (MQ=255)
                           aaGAATTTTTCTGTACAAAAACCCTCTTGTTTTATTATCCAAAAAAAGTGTAACGCAAGAAAACTTTTTTTGTGAAATGACGGATTTCACATTTTTAATTAAGAGAAAGGGAGAAAAAATAAGAAAAAAAAGCCCGCCATAACAGCg                                                                                                               <  2:1036836/147‑1 (MQ=255)
                             aaaTCTTTCTGTTAAAAAACCCCCTTGTTTAATTATCCGAAAAAGGTGTACCCCCAGATGACTTTTTTTTTGTAATGAGGGATTTCACTTTTTTAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCagg                                                                                                         <  1:581287/150‑1 (MQ=255)
                                   ttttGTTAATAAACCCTCCTTTTTTTGTACCCGAAAAAAGTGTACCGCCAGATGACTTTTTTTGTGTAATGACGGATTTCATTTTTTTAATTTAAAACAAAGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCaggagg                                                                                                      <  1:639831/145‑1 (MQ=255)
                                    tCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                  >  2:441024/1‑151 (MQ=255)
                                           aaaaacCCTCTTTTTTTAGTATCCGCAAAAAGTGAACCGCCAGATGACATTTTTTGTGTAAGGACGGATTCCATATTTTTAATTTAGAAAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa                                                                                              <  1:195787/146‑1 (MQ=255)
                                            aaaacCCTCTTGTTTTAGTACCCGAAAAAAGTGAAACCCCAGATGACTCTTTTTGTGAAAGGACGGATTTCACTTTTTTAATTTAAACAAAAGGAGGAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                          <  1:539399/150‑1 (MQ=255)
                                              aaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                         >  2:807766/1‑150 (MQ=255)
                                                    cctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                                                    >  1:749891/1‑149 (MQ=255)
                                                        ttttAGTATCCGCAAAAAGTGTAACGCAAGATAACTCTTTTTGTGTAAGGCGGGAGTTCACATTTTTAATTTAGATAAAAGGAGAAAAAAAAACAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                               <  2:642975/150‑1 (MQ=255)
                                                        ttttAGTACCCGAATAAAGTGTAACGCAAGATGACCTTTTTTGTGTAAGGACGGAGTCCACTTTTTTAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                                                                <  1:34648/149‑1 (MQ=255)
                                                               aTCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                        >  1:276166/1‑150 (MQ=255)
                                                                    cATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                   >  1:653951/1‑150 (MQ=255)
                                                                       aaaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatc                                                                >  2:790944/1‑150 (MQ=255)
                                                                          gtgtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                            >  1:57914/1‑151 (MQ=255)
                                                                            gaaaCGCCAGATGACTTTTTTTGTGTAATGACGGAGTTCACTTTTTAAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                            <  1:653644/147‑1 (MQ=255)
                                                                              aaCCCCAGATAACCTTTTTTTGTTAATGACGGATTTAACTTTTTTAATTAAGAGAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                        <  1:941315/151‑1 (MQ=255)
                                                                                cGCCAAATAACATTTTTGGTGAAATGAGGAAATACACTTTTTTATTTAAGAAAAAGGGAGGAAAAAAAAGAAAAAAACCACCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTa                                                         <  2:226851/148‑1 (MQ=255)
                                                                                   cAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                   >  1:100527/1‑151 (MQ=255)
                                                                                        gACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAtt                                              >  2:83070/1‑151 (MQ=255)
                                                                                             tttttGTGAAAGGAGGGGTTTCATTTTTTAAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAtt                                              <  2:920475/146‑1 (MQ=255)
                                                                                               tttGTGTAAGGAGGGAGTTCATTTTTTAAATTAAGAAAAAGGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc                                        <  2:803138/150‑1 (MQ=255)
                                                                                                      tAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                   >  2:754139/1‑149 (MQ=255)
                                                                                                        attgccggagtcctaTTTTTAATTTAAAACAAAAGGAGAAAAAATAGGAAAAAAAATCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGaa                                  <  2:673976/135‑1 (MQ=255)
                                                                                                           acgaggTTCACATTTTTAATTTAAACAAAGAAGGAAAAAAAAGAAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGt                             <  1:402705/145‑1 (MQ=255)
                                                                                                            cGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTc                            >  1:1008123/1‑150 (MQ=255)
                                                                                                                gTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCtt                       >  2:349648/1‑151 (MQ=255)
                                                                                                                       tttttAATTAAGATAAAGGAAAGAAAAATAACGAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa                 <  2:301045/150‑1 (MQ=255)
                                                                                                                           tAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg            >  1:572770/1‑151 (MQ=255)
                                                                                                                            aaTTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt           <  2:198247/151‑1 (MQ=255)
                                                                                                                             aTTAAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGtt          <  2:612899/151‑1 (MQ=255)
                                                                                                                              tttAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg         >  1:1030809/1‑151 (MQ=255)
                                                                                                                                 aGATAAAAGGAGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAt       <  1:865611/150‑1 (MQ=255)
                                                                                                                                   aTCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt    >  1:991453/1‑151 (MQ=255)
                                                                                                                                       aaaGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat  >  2:829473/1‑150 (MQ=255)
                                                                                                                                       aaaGGAGAAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat  <  1:145865/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTAT  >  NC_000913/1435575‑1435857

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