Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,912	0	G	C	47.8%	63.4 / 24.8	47	noncoding (60/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base G (4/20);  new base C (22/0);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.94e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTA  >  NC_000913/925765‑926054                                                                                                                                                      |                                                                                                                                                 aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATa                                                                                                                                               <  1:155879/151‑1 (MQ=255)    ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGt                                                                                                                                            <  2:558621/151‑1 (MQ=255)        gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGcc                                                                                                                                        <  2:460911/151‑1 (MQ=255)            aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGtat                                                                                                                                    <  1:585375/151‑1 (MQ=255)             gAAATTCACAGGTAACATACTCCCCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGtata                                                                                                                                   <  2:629465/151‑1 (MQ=255)             gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTTGAACC‑CCCCGTACggt                                                                                                                                           >  1:227616/1‑141 (MQ=255)                  tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑GGGGGGGGAGTTGGACC‑CCCCAGACGTTGCCGgata                                                                                                                                   >  2:84793/1‑146 (MQ=25)                          aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCTCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAgg                                                                                                                      <  1:72679/151‑1 (MQ=255)                          aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAgg                                                                                                                      <  2:391849/151‑1 (MQ=255)                               aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGGTTGAACC‑CCCCCCACGGTGCCCGATACAAACCTTCCCGGCg                                                                                                                    >  1:208102/1‑148 (MQ=21)                                        gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑AGGGGGGGGGTTGAAAC‑CCCCACACGTTGCCGTATACACACTTTTCAGGGGTGCC‑CCTTCAg                                                                                                        >  1:239690/1‑150 (MQ=18)                                            accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAAAC‑CCCCCTACGGTGCCGTATACACACCTTTCAGGCGTGCT‑CCTTCCGcc                                                                                                      >  1:545603/1‑148 (MQ=18)                                             ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTc                                                                                                  <  1:445567/151‑1 (MQ=255)                                                      cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCt                                                                                         <  2:610429/150‑1 (MQ=255)                                                               agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTCGAACC‑CCCCCTACGTTGCCGTATACACACCTTTCAGGCGGGGT‑CCTTCAGCCACCCGGAC‑ACCCCCCCaaa                                                                                  >  1:780823/1‑148 (MQ=14)                                                                 acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGtt                                                                             <  2:615680/151‑1 (MQ=255)                                                                   agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTTGCAGC‑CCCGCCACGTTTTCTCAAAAACACCTTCCCTGCGGGCT‑CCTTCTCCCCCCCCGAG‑ACCCCCCCAAAATAtttt                                                                           >  1:17214/1‑151 (MQ=9)                                                                     aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTTGAAAC‑CCCCATACGTTTCCGTATACACACTTTCCCGGCGTGCT‑CCTTCAGCCACTCCGAC‑ACCCCACCAAAATGtttt                                                                           >  1:3460/1‑149 (MQ=14)                                                                     aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AAGGGGGGGGTCGAAAC‑CCCCATAAGTTTGCGCATACACACTTTTCCGGCGGGCT‑CCTTCTGCCCCCCCGGC‑ACCCCCCCCAAATGtttt                                                                           >  2:217563/1‑149 (MQ=9)                                                                      tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTCGAACC‑CCCCCTACGTTGCCGTATAAACACTTTCCAGGCGTGCT‑CCTTCAGCCACCCCGGC‑ACCCCACCAAAATGTTTTTCt                                                                        >  2:800990/1‑151 (MQ=14)                                                                       aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑AGGGGGGGGTTCGCAAC‑CCCCCTAACTTTCCCTCTACACAATTTCCCCGCGCGGT‑CCTTCACCCACCCGGACAACCCCCCCAAATT‑TTTTGCt                                                                        >  1:154507/1‑150 (MQ=9)                                                                        aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTCGAAAC‑CCCCCTAAGTTGCCGTATAAAAACCTTCCAAGCGTGGT‑CCTTCAGCCACTCCGAC‑CCCCCACCCAAATGtttttgt                                                                        >  2:528506/1‑146 (MQ=14)                                                                          aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCAcc                                                                                     <  1:657226/134‑1 (MQ=37)                                                                             cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAAc                                                                 <  2:3460/151‑1 (MQ=255)                                                                                aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑AGGGGGGGAGTTGAACC‑CCCCACACGTTTCCGTTTAAACACTTTCCCGGGGGGGT‑CCTTCAGCCACTCGGGC‑ACCCCCCCAAAATATTTTTCTACCcaa                                                                  >  1:300581/1‑147 (MQ=2)                                                                                 aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTCGAAAC‑CCCCATACGTTGCCGGATACACACTTTCCCGGCGGGGT‑CCTTCAGCCCCACCGGC‑ACCCCACCCAAATTTTTTTCTACCaaaac                                                                >  2:448306/1‑146 (MQ=14)                                                                                   aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCCATACGTTGCCGGATACACACTTTTCCGGCGTGCT‑CCTTCCGCCACTCGGGC‑ACCTCACCCAATTGTTTTGCTACCAAAcc                                                                >  2:59470/1‑146 (MQ=14)                                                                                      aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTCGGACC‑CCCCATACGTTGCCGTATACACACTTTTCCGGCGTGGG‑CCCTCCGCCACCCCGAC‑ACCTCCCCCAAATGTTTTTCTACCAAAACTCATgg                                                          >  2:572550/1‑149 (MQ=14)                                                                                            cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGCGGGGGGGGGTTTAAACCAACCCAT‑CGTT‑CCGGAAAAACACTTTCCCGGCGGGGTCCCTCCACCCACCCACA‑‑ACCTCCCCCAATTTTTTTTTTACCAAACCTCATCGgggg                                                      >  1:442404/1‑147 (MQ=2)                                                                                              cTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTTGGACC‑CCCCACATGTTTCCGGTTATACACCTTTCAGGCGTGGT‑CCTTTAGCCCCCCGGAC‑ACCCCCCCAAAATTTTTTTCTTCCAAACCTCCTGGGGGGGaa                                                   >  1:73590/1‑148 (MQ=2)                                                                                                ttCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACggg                                               <  1:310073/150‑1 (MQ=255)                                                                                                 tCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTTGGACC‑CCCCATACGGTGCCGTATACACACTTTCCAGGCGGGCT‑CCTTCAGCCACTCGGAC‑ACCCCACCCAAATGTTTTTCTTCCAAACCTCCTGGGTGGCAAAgggg                                              >  2:175074/1‑150 (MQ=14)                                                                                                   gTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGTTCGAACC‑CCCCACTCGTTTCCGTATAGAAACCTTCCCGGCGTGCG‑CCTTTCGCCACTCGGGC‑ACCCCTCCAAAATGTTTTTCTTCCCAACCTCCTGGGGGGCGACggg                                               >  1:787984/1‑147 (MQ=9)                                                                                                    tACGAGCTCTTCTTTAAATATGGCGGGG‑AGGGGGGGGGTTTAAAC‑CCCCACACGTTGCCGTATACACAATTTCCAGGCGGGCT‑CCTTCAGCCACCCGGGC‑ACCCCACCCAAATGTTTTTCTATCCAACCTCATGGGGGGGAACgggggg                                            >  1:229195/1‑147 (MQ=9)                                                                                                      cGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGGTCGAACC‑CCCCATACGTTTCCCTATACACAATTTTCAGGCGTGGT‑CCTTCAGCCCCTCGGAC‑ACCCCACCAAAATTTTTTGCTACCAAACCTCATGGGTGGCAAAGGGGCGctac                                        >  2:333483/1‑151 (MQ=17)                                                                                                       gAGCTCTTCTTTAAATATGGCGGTG‑GGGGGGGGGGTCGAACC‑CCCGATACGTTTCCGTATACACACATTCCAGGCGGGGT‑CCTTCAACCAAGCGGGG‑ACCTCACCCAAATGTTTTGCTACCAAAACTCACCGGGGGCGAAgggggg                                            >  1:359454/1‑144 (MQ=9)                                                                                                               ctgtAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCTATACGTAGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAACCCTCATTGGTGGCAACGGGGCGCTACTATAgg                                  <  1:751536/145‑1 (MQ=255)                                                                                                                  tAAATATGGCGGTG‑AGGGGGGGGTTTGAACC‑CCCCCCAAGGTTCCCTATTAACACTTTTCAGGCGGGGG‑CCTTCCACCCCTCCGGC‑CCCCCACCAAAATGTTTTTTTTCCAAACCTCACGGGGGGGAACGGGGGGCTCCTATAggg                                 >  1:277188/1‑146 (MQ=255)                                                                                                                  tAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                             <  2:63959/150‑1 (MQ=255)                                                                                                                   aaaTATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTgg                           <  2:594827/151‑1 (MQ=255)                                                                                                                         ggCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa                     <  1:54095/151‑1 (MQ=255)                                                                                                                           cGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc                   <  2:429875/151‑1 (MQ=255)                                                                                                                            ggTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg                  >  2:65891/1‑151 (MQ=255)                                                                                                                               g‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc               <  2:41205/151‑1 (MQ=255)                                                                                                                                 aGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc               >  1:55863/1‑150 (MQ=255)                                                                                                                                  gggggggATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa             <  1:254806/151‑1 (MQ=255)                                                                                                                                             gAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTa  <  2:536552/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 ACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGATTCGAACC‑CCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGAC‑ACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTA  >  NC_000913/925765‑926054

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