Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,610,466	0	T	G	50.9%	66.8 / 89.4	53	G351G (GGT→GGG)	prfC	peptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base T (6/20);  new base G (27/0);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 1.14e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-09

AAATTCAGGCCAACATGGACCCGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGCACCAGAACTGTTCCGTCGTATC  >  NC_000913/4610318‑4610616                                                                                                                                                     |                                                                                                                                                         aaaTTCAGGCCAACATGGACCCGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGa                                                                                                                                                        <  2:851025/151‑1 (MQ=255)          ccAACATGGACCCGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTgggcggtcg                                                                                                                                                 >  1:172569/1‑141 (MQ=255)          ccAACATGGACCCGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCggt                                                                                                                                                   >  2:656901/1‑145 (MQ=255)                     ccGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGaagaag                                                                                                                                    <  2:270569/151‑1 (MQ=255)                      cGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGTTTCCAAGGTgaag                                                                                                                                       >  1:362991/1‑147 (MQ=255)                      cGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGACGGTCGCGCGTTgaaga                                                                                                                                      >  1:537793/1‑148 (MQ=255)                       gAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCCCGTTGAAGAAGCg                                                                                                                                  >  2:464468/1‑151 (MQ=255)                        aaaCACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGt                                                                                                                                 <  2:45129/151‑1 (MQ=255)                        aaaCACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGTTTGCCCGTTggagaa                                                                                                                                     >  1:672324/1‑147 (MQ=255)                          acacCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGTCCGTTCTCAcacttacgcaccgggt                                                                                                                               >  1:758443/1‑136 (MQ=255)                           cacCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATc                                                                                                                              <  2:99334/151‑1 (MQ=255)                             ccgcgaccgcgTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCg                                                                                                                            <  2:536984/151‑1 (MQ=255)                              cgcgaccgcgTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGAGCGGTCGCCCGTTGAAGAAACgggaaccg                                                                                                                           >  1:540532/1‑144 (MQ=255)                                  accgcgTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGTTCGCACCGTGGAGAAGGGGGTCCggg                                                                                                                          >  1:325113/1‑148 (MQ=255)                                         ggCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTTGCACGTTGAAGAAGCGGATCCGGGCggggtcc                                                                                                                  >  1:563137/1‑143 (MQ=255)                                                aTGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTTGCACGTTGAAGAAGGGTATCCGGGGGGTATCCTccg                                                                                                              >  2:41178/1‑144 (MQ=255)                                                 tGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCACGGTGAAGAAGAGGATCCGGGCGATATCCTCGGcc                                                                                                            >  2:279591/1‑147 (MQ=255)                                                   cGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGcaca                                                                                                      <  2:157426/151‑1 (MQ=255)                                                       tggtgTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCCCGTTGAAGAAGCGGATCCGGGCGGTATCCTCGGCCTGCACaa                                                                                                     >  2:114612/1‑148 (MQ=255)                                                           gtCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGGTCGCACGCTGGAGAAGAGAAACCGGGCGATATACTTGGCCCGCCCcaccac                                                                                                 >  1:517407/1‑148 (MQ=255)                                                            tCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGGCCGGTTTCACCGTGGAGGAAAGTATCCGGGCGGTATCCTCCGCCCGcaccaac                                                                                                   >  1:254318/1‑141 (MQ=255)                                                            tCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTCGCACGTTGAAGAAACGTATCCGGGCGGTATCCTCCGCCCTCACaaac                                                                                                   >  2:674959/1‑143 (MQ=255)                                                             ccGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGACGCTTCCACGTTTGAGAAAAGGATCCGGGCGGGGTCCTCGGCCTGCACAACCACgg                                                                                               >  1:799240/1‑148 (MQ=255)                                                               ggTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTTGCGCGTGGGAGAAGGGGAGCCCGGGGGGATCCTCGGCCTGCCCAAACCCGGCAcc                                                                                           >  1:411242/1‑150 (MQ=255)                                                                gTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTCGCCCGTTGAAGAAGCGGATCCGGGCGGTATCCTCGGCCTGCCCAACCACGGGAcc                                                                                           >  2:797521/1‑149 (MQ=255)                                                                    atatGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAg                                                                                     <  2:573563/151‑1 (MQ=255)                                                                    atatGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTTGCCCGGTGTTGGAGCGTATCCGGGCGATATTCTCCGCCCGCCCAACCACCGCACCAtt                                                                                        >  2:526719/1‑148 (MQ=255)                                                                      atGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGAt                                                                                   <  1:913602/151‑1 (MQ=255)                                                                        gAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCg                                                                                 <  1:454853/151‑1 (MQ=255)                                                                        gAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCg                                                                                 <  1:850796/151‑1 (MQ=255)                                                                               gCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTCGCCCGGTGGAGAAGGGGGTCCGGGCGGTATCCTCTGCCTGCACAACCACCGCCCCCCTCAGaaccgggaccc                                                                          >  1:507370/1‑141 (MQ=255)                                                                                   gAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTTGCACGTTGAAGAAAGGTATCCCGGGGGTATTCTCGGCCTGCACAACCACCGGCCCAATCAGATCGGCgcccccct                                                                       >  1:115126/1‑143 (MQ=255)                                                                                    aaaCTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTCGGACGCTTAAGAAAAGGAAACGGGGGGTATTCTCGGCCTGCACAACCACGGCACCATTCAAATCGGCGACA‑CCttt                                                                     >  2:473158/1‑150 (MQ=255)                                                                                         gcgcCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCa                                                                <  1:922026/150‑1 (MQ=255)                                                                                         gcgcCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCACGTTGGAGAAGCGTATCCGGGCGGTATCCTCGGCCCGCACAAACACGGCACCAATCAGATCGGGGACA‑CCCTTAccc                                                                 >  2:188223/1‑149 (MQ=255)                                                                                           gcCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAggg                                                             <  2:656335/151‑1 (MQ=255)                                                                                              aGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGTTCGCACGTTGTAGAAGCGGATACCGGGGGTATTCTCGGCCTGCACAACCACCGCCCCATTCCGCTCGGCGGCA‑CCTTTAcccccgggg                                                           >  1:910120/1‑144 (MQ=255)                                                                                              aGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGGTCGCCCGGTGAAGAAGCGGATCCGGGGGGTATCCCCGGCCCGCACCACCACCGCACCCCTCCGATCGGCGACACCCTTTACCCggggg                                                            >  1:873480/1‑149 (MQ=255)                                                                                                gTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTgaga                                                        <  2:914877/151‑1 (MQ=255)                                                                                                    gcACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGAACGTTCGCACGTTGAAGAAGCGGATCCCGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGGGACC‑CCTTTTccccagggga                                                          >  2:61772/1‑138 (MQ=255)                                                                                                        tGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTCGCACGGTGTAGAAGAGGATTCGGGCGCTAATCTCCGCCCGCACAACCACCGCACCAATCAGATCGGCGGCA‑CCCTTACCCCGGGGGGGGAGATGaa                                                 >  2:654081/1‑150 (MQ=255)                                                                                                                  gtggtgATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGt                                       <  2:230750/150‑1 (MQ=255)                                                                                                                    ggtgATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑tt                                   <  1:725000/151‑1 (MQ=255)                                                                                                                      tgATCTCCGACGCGCTGACCTTTATGGCGGGGGAACGTTCGCACGGTGAAGAAGCGGATCCGGGCGGTATCCCCGGGCTGCACAAACACGGCACCATTCAGAACGGCGACA‑CCTTTACCCAAGGGGAGATGATGAAAATCACCGGGA‑tttc                                 >  2:648551/1‑149 (MQ=255)                                                                                                                       gATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGCTTGCACGTTTAAGAAAAGGATCCCGGCGGGATCCTCGGCCTGCACAACCACCGCACCATTCAGATCGGCGGCA‑CCCTTATCCAGGGGGAGATGATGAAGGTCACCGGGA‑tt                                   >  2:163975/1‑148 (MQ=255)                                                                                                                           tcCGACGCGCTGACCTTTATGGCGGGGGGGCGTTCGGACGGTGGAGGAACGGATCCTGGGGGGGTTCTCGGGCCGCACCACAACCGGACCCTTTAGATCGGCGGCA‑CCTTTCACCCGGGGGAGAAGATGAAGATCAcccggat                                     >  1:97344/1‑138 (MQ=255)                                                                                                                            ccGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGaa                              <  1:561146/148‑1 (MQ=255)                                                                                                                             cGACGCGCTGACCTTTATGGCGGGGGGCCGTTTGCCCCTTGAAGAAGCGGATCCGGGCGATATCCTCGGCCTGCACAACCACCGGACCATTCAGATCGGCGACC‑CCTTTTCCCCGGGGGGGGAGATGAAAGTCACCGGGATTTCCGaatt                            >  1:951173/1‑148 (MQ=255)                                                                                                                               aCGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGc                        <  1:284964/151‑1 (MQ=255)                                                                                                                                                gCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGCACCAGAACTGTTCcgtc       <  2:367788/151‑1 (MQ=255)                                                                                                                                                 cGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGCACCAGAACTGTTCcgtcg      >  2:204939/1‑151 (MQ=255)                                                                                                                                                  gggTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGCACCAGAACTGTTCcgtcgt     <  2:838377/151‑1 (MQ=255)                                                                                                                                                   ggTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGCACCAGAACTGTTCCGTCGTa    <  1:186086/151‑1 (MQ=255)                                                                                                                                                     tGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGCACCAGAACTGTTCCGTCGTATc  >  2:586360/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                         AAATTCAGGCCAACATGGACCCGAAACACCGCGACCGCGTGGCGTTTATGCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACA‑CCTTTACCCAGGGTGAGATGATGAAGTTCACCGGTA‑TTCCGAACTTCGCACCAGAACTGTTCCGTCGTATC  >  NC_000913/4610318‑4610616

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