Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,035,757	0	A	C	59.2%	30.5 / 54.4	49	G150G (GGT→GGG)	prfB	peptide chain release factor RF‑2
Reads supporting (aligned to +/- strand):  ref base A (14/6);  new base C (0/29);  total (14/35)
Fisher's exact test for biased strand distribution p-value = 5.74e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-06

TCCGGATAAACAAACGCGGAGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCA‑GGCGGTGAACGCCGGTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCA  >  NC_000913/3035608‑3035900                                                                                                                                                         |                                                                                                                                                tCCGGATAAACAAACGCGGAGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCA‑GGCGGTGAACGCCGGTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATAc                                                                                                                                                <  1:117685/151‑1 (MQ=255)                     gCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCA‑GGCGGTGAACGCCGGTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGa                                                                                                                             >  1:219674/1‑150 (MQ=255)                      cTGAACGACGTGGGGCGACGCCCGCCGGAGTAAAACGGGTTTTTACGCCCCG‑GGCGGTGACCGCCGTTTTC‑TGTCCGCACCCAGCCGTAAGCGTAACCGCCGGAGATTTTGATCGTCACGGA‑TTTAATCCCCGCCACTTCACCTTCCGAc                                                                                                                            <  2:691018/150‑1 (MQ=255)                        gAACGGGGTGTGGCGACGCCCGCCGGAGTCAACCGGGTTTTTCCGCACCG‑GGGGGTGAACGCCGTTTTC‑TGTCCGCAGCCAGCCGTAAGCGTAACCGCCGGAGTTTTTGATCGTCACGAA‑TTTAATCCCCGCCACTTCACCTTCCGActc                                                                                                                          <  2:726688/150‑1 (MQ=255)                                gggggCGACGCCCGCGGGA‑TTAAACGGGTTTTTTCCCCCCAGGGGGGTGAACGCGGTTTTC‑TGTCCGCACCCACCCGTAAGCGTAACCCCCGGAGATTTTGTTCGTCAGGGT‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCgatga                                                                                                                  <  1:647331/147‑1 (MQ=255)                                 gtgGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCA‑GGCGGTGAACGCCGGTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGAtt                                                                                                                >  2:370441/1‑151 (MQ=255)                                            cgccgGAGTCAAACGGGCTTTTACGCACCA‑GGCGGTGAACGCCGGTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGaa                                                                                                     >  2:264910/1‑151 (MQ=255)                                              cgggAGTCAACCGGGTTTTTCCGCACCG‑GGCGGTGACCGCCGTTTTC‑TGTACGCAGCCAGCCGTAAGCGTAACCGCCGGAGATTTTGATCGTCACGGA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAAc                                                                                                   <  2:75806/149‑1 (MQ=255)                                              cgggAGTCAAACGGGTTTTTCCGCCCCG‑GGCGGTGAACGCCGTTTTC‑TGTCCGCAGCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAAc                                                                                                   <  2:322444/149‑1 (MQ=255)                                                           ggCTTTTACGCACCA‑GGCGGTGAACGCCGGTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTg                                                                                       >  1:207074/1‑150 (MQ=255)                                                             tttttCCGCCCCG‑GGCGGTGAACGCCGTTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGccc                                                                                    <  2:312163/150‑1 (MQ=255)                                                             tttttCCGCACAG‑GGGGGTGAACCCCGTTTTC‑TGTCCCCCCCCAGCCGTAGGCGTAACCGCCGGATATTTTGATCGTCGCGTT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGccc                                                                                    <  2:296193/150‑1 (MQ=255)                                                              tttttccgcccggggCGGTGAACGCCGTTTTC‑TGAACGCACCCAGCCGTAAGCTTAATCGCCGGAGATTTTGACCGTCACGGA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGccc                                                                                    <  1:14499/138‑1 (MQ=255)                                                                  ccGCACCG‑GGCGGTGACCGCCGTTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGTA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc                                                                               <  2:618004/150‑1 (MQ=255)                                                                      ccccgggCGGTGAACGCCGTTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGa                                                                            <  2:591375/146‑1 (MQ=255)                                                                       ccA‑GGCGGGAACCGCCGTTTTC‑TGTACCCACCCACCCGAAAGCTTACTCGCCGGAGTTTTTGACCGTC‑CGGATTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGata                                                                          <  1:222024/151‑1 (MQ=255)                                                                             gggTGACCGCGGTTTTC‑TGTCCGCAGCCAGCCGTAAGCGTAACCGCCGGAGTTTTTGATCGTCACGAT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGatacata                                                                      <  2:476659/149‑1 (MQ=255)                                                                              ggTGACCGCCGTTTTT‑TGTCCGCACCCAGCCGTAAGCGAAACCGCCGGAGTTTTTGTTCGTCACGGT‑TTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGatacatac                                                                     <  2:403901/150‑1 (MQ=255)                                                                                taaaCGCGGTTTTC‑TGTCCGCCCCCACCCGTAGGCGAAACCGCCGGAGTTTTTGATCGTCACGGT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGc                                                                   <  1:67182/148‑1 (MQ=255)                                                                                tGAACGCGGTTTTC‑TGTACGCAGCCAGCCGAAAGCGTAACCGCCGGAGATTTTGATCGTCACGGA‑TTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCt                                                                  <  2:312311/151‑1 (MQ=255)                                                                                      cctttttC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCa                                                            <  2:556777/147‑1 (MQ=255)                                                                                      ccGTTTTC‑TGAACCCACCCAGCCGAACGCTTACCCCCCGGATTTTTTGACCGTCCAGGA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCa                                                            <  1:523008/151‑1 (MQ=255)                                                                                       ccgTTTCTTGTCCGCACCCGGCCGAAAGCTAAACCGCGGGATTTTTTGTTCGTCACGGT‑TTTAATCCCCGCCACTTCACATTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGc                                                             <  1:326685/148‑1 (MQ=255)                                                                                          tttC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATg                                                          >  2:278506/1‑149 (MQ=255)                                                                                           ttC‑TGTCCGCAGCCAGCCGAAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTc                                                       <  2:447112/151‑1 (MQ=255)                                                                                               tGTCCGCAGCCAGCCGTAGGCGTAACCGCCGGAGATTTTGATCGTCACGGA‑TTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGcg                                                    <  2:90713/151‑2 (MQ=255)                                                                                                 tACGCCGCCCGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGccc                                                   <  1:236158/150‑1 (MQ=255)                                                                                                  aCGCAGCCACCCGAAGGCGAAATCGCGGGAGTTTTTGTTCGTCCCGGT‑TTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAg                                                 <  2:724682/151‑1 (MQ=255)                                                                                                    gCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGt                                                >  1:212793/1‑150 (MQ=255)                                                                                                    gCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGt                                                >  1:601175/1‑150 (MQ=255)                                                                                                     ccgccagcgGAAGGCGAAATCGCCGGAGTTTTTGTTCGTCACGTA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTcc                                              <  1:366992/149‑1 (MQ=255)                                                                                                     cacccggccGTAAGCGTATCCGCCGGAATTTTTGATCGTCCCGTT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTcc                                              <  1:467399/145‑1 (MQ=255)                                                                                                       gccagccGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGTCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTc                                               >  2:681977/1‑148 (MQ=255)                                                                                                       gccagccGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTc                                               >  2:100769/1‑148 (MQ=255)                                                                                                        ccagccGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGT‑TTAAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCGGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTcc                                              <  2:253016/148‑1 (MQ=255)                                                                                                        ccagccGAAAGCGTAATCGCCGGAAATTTTGATCGCCACGGT‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTcc                                              <  1:494463/148‑1 (MQ=255)                                                                                                            ccGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCtt                                       >  1:226946/1‑151 (MQ=255)                                                                                                            ccGTAAGCGTAACCGCCGGAGATTTTGATCGTCACGGT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGc                                         <  2:653609/149‑1 (MQ=255)                                                                                                               tAAGCGTAATCGCGGGAGTTTTTGATCGTCACGAT‑TTAAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCg                                    <  2:632196/151‑1 (MQ=255)                                                                                                               aaaGCGTAACCGCCGGAGATTTTGTTCGTCACGGA‑TTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCg                                    <  2:207076/150‑1 (MQ=255)                                                                                                                   cGTAACCGCCGGAGTTTTTGATCGTCACGTA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTAc                                 <  2:218894/150‑1 (MQ=255)                                                                                                                   cGTAACCGCCGGAGATTTTGTCCGTCACTGT‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTAc                                 <  1:536899/150‑1 (MQ=255)                                                                                                                             ggAGTTTTTGATCGTCAGGTT‑TTAAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGAccc                       <  2:186334/150‑1 (MQ=255)                                                                                                                                       accGTCACGTA‑TTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAatat            <  1:249237/149‑1 (MQ=255)                                                                                                                                            cACGAA‑TTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAgg       <  2:435132/151‑1 (MQ=255)                                                                                                                                                gA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAg    >  2:576347/1‑150 (MQ=255)                                                                                                                                                gA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAg    >  1:476754/1‑150 (MQ=255)                                                                                                                                                 a‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCa  >  2:519376/1‑151 (MQ=255)                                                                                                                                                 a‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCa  >  1:227860/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                TCCGGATAAACAAACGCGGAGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCA‑GGCGGTGAACGCCGGTTTC‑TGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGA‑TTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCA  >  NC_000913/3035608‑3035900

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