Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,035,757	0	A	C	43.4%	87.7 / 11.6	53	G150G (GGT→GGG)	prfB	peptide chain release factor RF‑2
Reads supporting (aligned to +/- strand):  ref base A (25/5);  new base C (0/23);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 2.61e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGATAAACAAACGCGGAGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCAGTC  >  NC_000913/3035608‑3035903                                                                                                                                                      |                                                                                                                                                   tCCGGATAAACAAACGCGGAGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATAc                                                                                                                                                   >  2:741136/1‑151 (MQ=255)             aacGCGGAGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCAc                                                                                                                                       >  1:159439/1‑151 (MQ=255)                  ggAGCGGAACGCGGGGGGGCGCCGCCCGCGGGAGTCAAGGGGTTTTTTACGCCCCAGGCGGTGACCGCCGTTTTCTGTACGCAGCCAGCCGTAGGCGTAATCGCCGGAGTTTTTGATCGCCAGGGTTTAAATACCCGCCACTTCACCtt                                                                                                                                    <  2:422580/149‑1 (MQ=255)                   gggcgaaaCGGCGGGTGGCGCCGCCCGCGGGAGTCAACCGGTTTTTCCCGCACGGGGGGGTGAGCGCGGGTTTCTGTCCGCAGCCAGCCGAAAGCTTAACCGCGGGAGATTTTGAGCGTCGCGGTTTAAATACCCGCCACTTCACCTTCCg                                                                                                                                 <  2:730042/145‑1 (MQ=255)                    aGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCg                                                                                                                                 >  2:431221/1‑150 (MQ=255)                      cTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGAct                                                                                                                              >  1:293783/1‑151 (MQ=255)                         accggcggggggCGACGCCCGCGGGATTCAACCGGTTTTTTCCCCACCAGGGGGTGACCGCCGTTTTCTGTACGCACCCAGCCGAAAGCGTAATCGCGGGAGATTTTGATCGTCACGGATTAAATCCCCGCCACTTCACCTTCCGACTCtt                                                                                                                           <  1:450501/141‑1 (MQ=255)                         aaCGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCtt                                                                                                                           >  2:180109/1‑151 (MQ=255)                            gacgTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCg                                                                                                                         >  1:227452/1‑150 (MQ=255)                               gtggGGCGACGCCCGCGGGAGTAAAGGGGGTTTTTACGCCCCGGGGGGTGACCGCCGTTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGTTTTAATACCCGCCACTTCACCTTCCGACTCTTCgatga                                                                                                                     <  1:683537/151‑1 (MQ=255)                               gggtgGCGCCGCCCCCCGGAGCCAAGCGGGTTTTACCGCCCGGGGCGGTGAACCCCGTTTTTTGTCGGCAGCCAGCCGTAAGGGTAACCGCCGGAGATTTTGATCGTCCCGAATTTAATCCCCGCCACTTCACCTTCCGACTCTTc                                                                                                                          <  2:106888/144‑1 (MQ=255)                                 gtgCCGACGCCCGCCGGAGTCAAGCGGGTTTTTACGCCCCGGGGGGTGACCGCCGTTTTTTGTACGCGCCCAGCCGTAGGCGTATTCGCCGGAGTTTTTGATCGTCACGGTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCgatga                                                                                                                     <  2:764844/149‑1 (MQ=255)                                  tgGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGAttt                                                                                                                  >  2:586373/1‑151 (MQ=255)                                   ggCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTc                                                                                                                 >  2:235366/1‑151 (MQ=255)                                         gacCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGtttt                                                                                                           >  2:52164/1‑151 (MQ=255)                                           ccgccgGAGTCAACCGGTTTTTTACGCCCCGGGGGGTGACCGCGGTTTTCGGCACGCACCCAGCCGAAAGCTTAACCCCCGGAGTTTTTGATCGTCACGTATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGtt                                                                                                             <  2:850059/147‑1 (MQ=255)                                                ggagcaaaaGGGTCTTTTCCGCACCGGGGGGTAACCGCCGTTTTCTGTCCGCACCCAGCGGAAGGCTTACCCGCCGGAGTTTTTGATCGTCACGATTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCa                                                                                                    <  1:513197/145‑1 (MQ=255)                                                ggATTCAAACGGTTTTTTACGCCCCAGGCGGTGAACGCCGTTTTCTGTACGCACCCAGCCGTAGGCGTAATCGCGGGAGTTTTTGTTCGTCACGGATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGaaa                                                                                                       <  1:218127/148‑1 (MQ=255)                                                  aGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcg                                                                                                  >  1:635209/1‑151 (MQ=255)                                                   gTCAAACGGTTTTTTCCGCCCCGGGGGGGGAACCCTGTTTTCTGTCCGCAGCCAGCCGAAAGCGTAACCGCCGGAGTTTTTGACCGTCAGGGTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAAcc                                                                                                     <  2:96231/147‑1 (MQ=255)                                                      aaaCGGGTTTTTACGCCCGGGGGGGTGACCCCCGTTTTCTGTCCGCAGCCAGCCGAAAGCGTAATCCCCGGAGTTTTTGCTCGTCACGGTTTTAAACCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAc                                                                                                   <  1:405973/146‑1 (MQ=255)                                                      aaaCGGGTTTTCACGCCCCGGGCGGTGACCGCCGTTTTTTTTCCGCAGCCAGCCGTAAGCTTAACCGCCGGAGTTTTTGATCGTCACGGTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAt                                                                                              <  1:200758/151‑1 (MQ=255)                                                       agcGGGTTTTTCCGCCCCGGGGGGTAACCGCCGGTTTCTGTCCGCAGCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGTATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAtt                                                                                             <  1:797574/149‑1 (MQ=255)                                                       aaCGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAtt                                                                                             >  1:669740/1‑151 (MQ=255)                                                           ggCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGc                                                                                         >  1:738533/1‑151 (MQ=255)                                                                     cACCGGGGGGTGACCGCGGTTTTCTGTACGCAGCCAGCCGTAAGCGAAATCGCCGGAGTTTTTGATCGTCAGGGATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGa                                                                               <  1:83199/151‑1 (MQ=255)                                                                           ggggTGACCGCCGTTTTCTGTCCGCACCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGatacata                                                                         <  1:355694/149‑1 (MQ=255)                                                                                gAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTc                                                                    >  1:57447/1‑151 (MQ=255)                                                                                    gCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGc                                                                >  1:680605/1‑151 (MQ=255)                                                                                     ccGTTTTCTGTCCGCACCCCGCCTTAAGCGTAATCGCCGGATTTTTTGATCGTCACGGATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCaa                                                                  <  1:658643/148‑1 (MQ=255)                                                                                              gTACGCACCCGGCCGTAAGCGTAATCGCGGGAGTTTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCg                                                         <  2:364364/148‑1 (MQ=255)                                                                                                ccGCACCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCCCGATTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGccc                                                      <  1:60819/148‑1 (MQ=255)                                                                                                  gCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTc                                                  >  1:592352/1‑151 (MQ=255)                                                                                                   cagccagcCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTcc                                                 >  1:118016/1‑151 (MQ=255)                                                                                                       cagccGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCtgtg                                             >  1:169601/1‑151 (MQ=255)                                                                                                          ccGAAAGCTTAACCGCCGGAGTTTTTGATCGTCACGTATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCtt                                          <  2:554244/151‑1 (MQ=255)                                                                                                              aaggttAATCGCCGGATATTTTGACCGTCACGGATTAAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAACCATGCTCGCCCAGTCCTGTGCTTCCGt                                      <  1:586456/146‑1 (MQ=255)                                                                                                                gCGAAATCGCCGGAGTTTTTGATCGTCACGGTTTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTa                                     <  1:139692/150‑1 (MQ=255)                                                                                                                    aaTCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGcc                                >  1:331172/1‑151 (MQ=255)                                                                                                                        gCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGAc                            >  2:805613/1‑151 (MQ=255)                                                                                                                         ccGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGAcc                           >  1:327293/1‑151 (MQ=255)                                                                                                                           ggAGTTTTTGATCGTCACGGTTTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGAcccc                         <  1:247912/151‑1 (MQ=255)                                                                                                                           ggAGTTTTTGACCGTCCCGTATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGAcccc                         <  2:696460/151‑1 (MQ=255)                                                                                                                           ggAGATTTTGATCGTCACGGATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGAccc                          <  2:505636/150‑1 (MQ=255)                                                                                                                           ggAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGAcccc                         >  2:447302/1‑151 (MQ=255)                                                                                                                                tttaGATCGTCCCGTATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCt                     <  2:822425/150‑1 (MQ=255)                                                                                                                                  ttGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGa                   >  2:803109/1‑150 (MQ=255)                                                                                                                                   tGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAat                 >  1:649542/1‑151 (MQ=255)                                                                                                                                        gTCACGGTTTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCg             <  1:715500/150‑1 (MQ=255)                                                                                                                                            cGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGt         >  2:843896/1‑150 (MQ=255)                                                                                                                                             gtttttACTCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTa        <  2:56847/147‑1 (MQ=255)                                                                                                                                                 ttAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCAGt   >  1:391800/1‑151 (MQ=255)                                                                                                                                                   aaTCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCAGTc  <  2:718325/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TCCGGATAAACAAACGCGGAGCTGAACGACGTGTGGCGACGACCGCCGGAGTCAAACGGGCTTTTACGCACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCTCGCCCAGTCCTGTGCTTCCGTACCGCCAGACCCCGCCTGAATATCGAGGTAGCAGTC  >  NC_000913/3035608‑3035903

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