Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6860GT41.0% 8.2 / 16.0 39intergenic (+40/+25)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (9/14);  new base T (2/14);  total (11/28)
Fisher's exact test for biased strand distribution p-value = 8.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCG  >  NC_000913/2151538‑2151826
                                                                                                                                                    |                                                                                                                                              
gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttgtt                                                                                                                                              <  2:471233/151‑1 (MQ=255)
gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttgtt                                                                                                                                              <  2:33665/151‑1 (MQ=255)
                                   atccatttGTTCTTTTATTGGTCCCCTACCAAAAAAAAGACGCCGACGAGTTAAGAAAAAAATTAGTTCCCAAATCAGCCGTTATATTTTTTGTTAAACCGGGCGTTTAAATTTTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCca                                                                                                          <  2:286629/144‑1 (MQ=255)
                                    taCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTGATTGTTG‑T‑CTTCAACCGCCATGAGAAAATCACACTTGCcac                                                                                                         >  1:598170/1‑150 (MQ=255)
                                       cATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGCGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTG‑T‑CTTTAAACGCCCACAGAAACTCACTCTGc                                                                                                             >  1:945042/1‑143 (MQ=255)
                                         ttGATCCTTTATTGGTCCCCTACCCAAAAGAAGAGCCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCACCCGTTATATTTATTTTTTAAACGGGGGTTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAgcg                                                                                                     <  2:633029/149‑1 (MQ=255)
                                                 ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa                                                                                            >  2:451358/1‑150 (MQ=255)
                                                    ttGTTCCCCTATCAAAAAGAAGAAGCCGACGAGCAAAGAAAAAAATTATTTCTCTAATCAGCCGTTATATTTATTTTTAAACCGGGCGGTTAATTTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACg                                                                                         <  1:560588/150‑1 (MQ=255)
                                                         ccccTATCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTGTTAAAACGGGCGGTTGATTTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTa                                                                                    <  2:665254/148‑1 (MQ=255)
                                                           ccTACCCAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCTCTAATCAGCCTTTATATTTATTGTTTAACCTGGCGTTTGATTTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatca                                                                                 <  2:974386/151‑1 (MQ=255)
                                                             tATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTG‑T‑CTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGGGCCAACGtttt                                                                                     >  1:305485/1‑145 (MQ=255)
                                                              gccaaaaaGAAAAACCCATCGAGATAAGAAAAAAATATTTTCGCTAATCACCTTTTATATTTATGTATTAAACTGGCGGTTAAATTTTTT‑T‑TTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatca                                                                                 <  1:236685/144‑1 (MQ=255)
                                                              aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGGGCCCACGTTTTatcat                                                                                >  1:853031/1‑149 (MQ=255)
                                                                  aaaaaaAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCGTTATATTTATTTTTAAAACGGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                <  2:985247/145‑1 (MQ=255)
                                                                    aagaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGt                                                                         >  2:1089512/1‑150 (MQ=255)
                                                                                  gaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccg                                                           >  2:874509/1‑150 (MQ=255)
                                                                                   aGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccgcc                                                         >  1:355513/1‑151 (MQ=255)
                                                                                    gcaaaGAGAAAATTTATTTGGCAAACCCGCCTTTATTTTTTTTTTTAAACCGGGCGTTTAATTTTTTG‑T‑CTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCACCCGTACGCTCGTTTAccgccg                                                        <  2:1055621/148‑1 (MQ=255)
                                                                                       aaGAGAAAAATTGTTTCCTTAATCACCCTTTATTTTTATT‑TTTAAACCTGGGGTTTAATTTTTT‑TCCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGc                                                     <  1:915903/151‑1 (MQ=255)
                                                                                        agaaaaaaaTTAGTTCTTAAATCCGCTTTTATTTTTTTTGTTTAAACGGGCGGTGTATTTTTGG‑T‑TTCCAGCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCt                                                    <  1:294200/146‑1 (MQ=255)
                                                                                          agCAAAATTTGTTCGCTAAGCCGCCTTTATTTGTTTTTTTAAAACTGGGGGTTAATTTTTTG‑T‑TTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCgg                                                       <  1:438292/146‑1 (MQ=255)
                                                                                                atttttgcGCTAATCCGCCTTTATATTTATTGTAAAAACGGGGGTTTGAATTGTTT‑T‑TTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaa                                            <  1:379342/144‑1 (MQ=255)
                                                                                                     ttCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaaa                                           >  2:553323/1‑147 (MQ=255)
                                                                                                            aTCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCtttta                                  >  2:665224/1‑147 (MQ=255)
                                                                                                                 ccGTTATATTTATTGTTAAATCGGGCGGTTTATTTTTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCa                           <  1:465733/151‑1 (MQ=255)
                                                                                                                   gTTATTTTTTTTGTTTAAACTGGGGTTTAAATTGTTG‑T‑TTAAAACCGCCATGCGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                            <  2:260865/148‑1 (MQ=255)
                                                                                                                    ttATTTTTTTTTTTAAACCGGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCtt                                     <  2:125911/138‑1 (MQ=255)
                                                                                                                    ttATATTTATTGTTAAACCTGGCGTTTGAA‑TGTTGTT‑CTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAgg                         <  1:485790/150‑1 (MQ=255)
                                                                                                                     tataTTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCTACCTCCATTAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAACGCCGGGCTGGACAAAAACTGCTTTaaaa                                >  1:5670/1‑140 (MQ=255)
                                                                                                                      atatTTATTTTTTAACCTGGGGGTTGATTTTTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                       <  2:684541/150‑1 (MQ=255)
                                                                                                                       tttttattGTATAACCGGGCGGTTAATTTTTTT‑T‑CTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAgg                         <  1:554429/145‑1 (MQ=255)
                                                                                                                       tatttattTTTTAACCGGGCGTTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                       <  2:465672/149‑1 (MQ=255)
                                                                                                                         ttttttttttAAACTGGGGGTTAATTTGTTT‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                       <  2:775312/147‑1 (MQ=255)
                                                                                                                         ttttttttaaaaccgggCGTTTAATTTGTTG‑T‑CTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaata                    <  1:916243/135‑1 (MQ=255)
                                                                                                                         tttattTTTTAATCTGGCGTTTAATTTTTTC‑T‑TTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaata                    <  2:139171/150‑1 (MQ=255)
                                                                                                                          tttttttttAAAGTGGCGGTTAATTTTTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaataa                   <  1:769823/144‑1 (MQ=255)
                                                                                                                             ttttttAACCGGGCGGTTAATTTTTTT‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                 <  1:589918/146‑1 (MQ=255)
                                                                                                                             ttGTTTAACCTGGCGTTTAATTTTTTG‑T‑TTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATc                <  1:921017/150‑1 (MQ=255)
                                                                                                                             ttGTTAAAACGGGGGGTTGATTTGTTG‑T‑TTTCGACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATc                <  2:267204/150‑1 (MQ=255)
                                                                                                                                         gCGGTTAAATTGTTT‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGc   <  2:359070/151‑1 (MQ=255)
                                                                                                                                           ggTTGAATTTTTT‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCg  <  2:700197/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTG‑T‑CTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCG  >  NC_000913/2151538‑2151826

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