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

GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCG  >  NC_000913/2151538‑2151826
                                                                                                                                                    |                                                                                                                                              
gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttgtt                                                                                                                                              <  2:471233/151‑1 (MQ=255)
gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttgtt                                                                                                                                              <  2:33665/151‑1 (MQ=255)
                                   atccatttGTTCTTTTATTGGTCCCCTACCAAAAAAAAGACGCCGACGAGTTAAGAAAAAAATTAGTTCCCAAATCAGCCGTTATATTTTTTGTTAAACCGGGCGTTTAAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCca                                                                                                          <  2:286629/144‑1 (MQ=255)
                                    taCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTGATTGTTGTCTTCAACCGCCATGAGAAAATCACACTTGCcac                                                                                                         >  1:598170/1‑150 (MQ=255)
                                       cATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGCGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTGTCTTTAAACGCCCACAGAAACTCACTCTGc                                                                                                             >  1:945042/1‑143 (MQ=255)
                                         ttGATCCTTTATTGGTCCCCTACCCAAAAGAAGAGCCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCACCCGTTATATTTATTTTTTAAACGGGGGTTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgcg                                                                                                     <  2:633029/149‑1 (MQ=255)
                                                 ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa                                                                                            >  2:451358/1‑150 (MQ=255)
                                                    ttGTTCCCCTATCAAAAAGAAGAAGCCGACGAGCAAAGAAAAAAATTATTTCTCTAATCAGCCGTTATATTTATTTTTAAACCGGGCGGTTAATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACg                                                                                         <  1:560588/150‑1 (MQ=255)
                                                         ccccTATCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTGTTAAAACGGGCGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTa                                                                                    <  2:665254/148‑1 (MQ=255)
                                                           ccTACCCAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCTCTAATCAGCCTTTATATTTATTGTTTAACCTGGCGTTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatca                                                                                 <  2:974386/151‑1 (MQ=255)
                                                             tATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGGGCCAACGtttt                                                                                     >  1:305485/1‑145 (MQ=255)
                                                              gccaaaaaGAAAAACCCATCGAGATAAGAAAAAAATATTTTCGCTAATCACCTTTTATATTTATGTATTAAACTGGCGGTTAAATTTTTTTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatca                                                                                 <  1:236685/144‑1 (MQ=255)
                                                              aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGGGCCCACGTTTTatcat                                                                                >  1:853031/1‑149 (MQ=255)
                                                                  aaaaaaAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCGTTATATTTATTTTTAAAACGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                <  2:985247/145‑1 (MQ=255)
                                                                    aagaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGt                                                                         >  2:1089512/1‑150 (MQ=255)
                                                                                  gaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccg                                                           >  2:874509/1‑150 (MQ=255)
                                                                                   aGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccgcc                                                         >  1:355513/1‑151 (MQ=255)
                                                                                    gcaaaGAGAAAATTTATTTGGCAAACCCGCCTTTATTTTTTTTTTTAAACCGGGCGTTTAATTTTTTGTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCACCCGTACGCTCGTTTAccgccg                                                        <  2:1055621/148‑1 (MQ=255)
                                                                                       aaGAGAAAAATTGTTTCCTTAATCACCCTTTATTTTTATTTTTAAACCTGGGGTTTAATTTTTTTCCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGc                                                     <  1:915903/151‑1 (MQ=255)
                                                                                        agaaaaaaaTTAGTTCTTAAATCCGCTTTTATTTTTTTTGTTTAAACGGGCGGTGTATTTTTGGTTTCCAGCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCt                                                    <  1:294200/146‑1 (MQ=255)
                                                                                          agCAAAATTTGTTCGCTAAGCCGCCTTTATTTGTTTTTTTAAAACTGGGGGTTAATTTTTTGTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCgg                                                       <  1:438292/146‑1 (MQ=255)
                                                                                                atttttgcGCTAATCCGCCTTTATATTTATTGTAAAAACGGGGGTTTGAATTGTTTTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaa                                            <  1:379342/144‑1 (MQ=255)
                                                                                                     ttCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaaa                                           >  2:553323/1‑147 (MQ=255)
                                                                                                            aTCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCtttta                                  >  2:665224/1‑147 (MQ=255)
                                                                                                                 ccGTTATATTTATTGTTAAATCGGGCGGTTTATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCa                           <  1:465733/151‑1 (MQ=255)
                                                                                                                   gTTATTTTTTTTGTTTAAACTGGGGTTTAAATTGTTGTTTAAAACCGCCATGCGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                            <  2:260865/148‑1 (MQ=255)
                                                                                                                    ttATTTTTTTTTTTAAACCGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCtt                                     <  2:125911/138‑1 (MQ=255)
                                                                                                                    ttATATTTATTGTTAAACCTGGCGTTTGAATGTTGTTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAgg                         <  1:485790/150‑1 (MQ=255)
                                                                                                                     tataTTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCTACCTCCATTAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAACGCCGGGCTGGACAAAAACTGCTTTaaaa                                >  1:5670/1‑140 (MQ=255)
                                                                                                                      atatTTATTTTTTAACCTGGGGGTTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                       <  2:684541/150‑1 (MQ=255)
                                                                                                                       tttttattGTATAACCGGGCGGTTAATTTTTTTTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAgg                         <  1:554429/145‑1 (MQ=255)
                                                                                                                       tatttattTTTTAACCGGGCGTTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                       <  2:465672/149‑1 (MQ=255)
                                                                                                                         ttttttttttAAACTGGGGGTTAATTTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                       <  2:775312/147‑1 (MQ=255)
                                                                                                                         ttttttttaaaaccgggCGTTTAATTTGTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaata                    <  1:916243/135‑1 (MQ=255)
                                                                                                                         tttattTTTTAATCTGGCGTTTAATTTTTTCTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaata                    <  2:139171/150‑1 (MQ=255)
                                                                                                                          tttttttttAAAGTGGCGGTTAATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaataa                   <  1:769823/144‑1 (MQ=255)
                                                                                                                             ttttttAACCGGGCGGTTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                 <  1:589918/146‑1 (MQ=255)
                                                                                                                             ttGTTTAACCTGGCGTTTAATTTTTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATc                <  1:921017/150‑1 (MQ=255)
                                                                                                                             ttGTTAAAACGGGGGGTTGATTTGTTGTTTTCGACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATc                <  2:267204/150‑1 (MQ=255)
                                                                                                                                         gCGGTTAAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGc   <  2:359070/151‑1 (MQ=255)
                                                                                                                                           ggTTGAATTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCg  <  2:700197/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCG  >  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: