Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	47.2%	88.9 / 14.8	53	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (28/0);  new base C (0/25);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.11e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAAT  >  NC_000913/3443134‑3443417                                                                                                                                                |                                                                                                                                             cATAGAGTAACACATAAAACGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGc                                                                                                                                       >  2:767343/1‑151 (MQ=255)  aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                      >  2:468433/1‑151 (MQ=255)              aTAAAGCGGTTGCCCAGGCTGCAAATACAGCAAAATTGTTGTCAGCCAGTTCCAACCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTc                                                                                                                          <  2:732589/151‑1 (MQ=255)                aaggCGTTTGCCCGGGCTAAAAAAACAGCAAATTTTTTGTCACCCAGTCCCACCCAGTCCCCCCCCCAAACCAGGCCGCGATGGCCGCCGGAAAAAAAATAATACTGGAAGCGAAGATTGCCGGAATCCCCCCCGCCATATTCACTTTCAg                                                                                                                        <  1:780470/148‑1 (MQ=255)                aaaGCGGTTGCCCAGGCTGAAAATACAGCAAATTTGTTGTCAGCCAGTTCCAACCAGTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAg                                                                                                                        <  2:440919/151‑1 (MQ=255)                   gCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCgg                                                                                                                     >  2:48626/1‑151 (MQ=255)                     ggTTCCCCGGGCTGAAAAAACAGCAAAATTGTGGTCAGCCAGTCCCAACCAGTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACGGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGt                                                                                                                    <  1:662229/150‑1 (MQ=255)                          cccAGGCTGAAAAAACAGCAAAATTGTTGCCAGCCATTCCCAACCAGTCCCCCCCCCGACCCATGCCGGGGTGGCCGCGGGGAACAAAAAAATACTGGAGGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                <  1:912225/149‑1 (MQ=255)                          cacaggcgaaaaaTACAGCAAATTTTTTGCCAGCCAGTCCCAACCAGTCCCCCCCCCGAACAAGGCCGGGATGGCCGCGGGGAAAAAAATAATCCTGGAAGGGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                              <  1:788854/141‑1 (MQ=255)                              ggCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc                                                                                                          >  1:726237/1‑151 (MQ=255)                                cTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                        >  2:150118/1‑151 (MQ=255)                                   cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                     >  2:314324/1‑151 (MQ=255)                                    aaaTAAAGCAAAATTGTTGTCAGCAATTCCCACCCAGACCCCCCCCCAAACCATGACGCGATGGTCGCGGGAAAAAAAATAATACTGGAAGCGAAATTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                       <  2:805678/148‑1 (MQ=255)                                             gAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                              >  1:274144/1‑148 (MQ=255)                                               aaTTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                          >  2:873354/1‑150 (MQ=255)                                                 ttgttgTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacg                                                                                       >  2:815923/1‑151 (MQ=255)                                                    ttgtCAGCCAGTTCCAACCGGCCCCCCCCCCGAACCATGCCGCGATGGTCGCCGGGAACAAAAAAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                     <  2:221351/150‑1 (MQ=255)                                                     tgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                    >  1:842263/1‑150 (MQ=255)                                                       tCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc                                                                                 >  2:360846/1‑151 (MQ=255)                                                         aGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                               >  2:828373/1‑151 (MQ=255)                                                          gCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTg                                                                              >  2:260565/1‑151 (MQ=255)                                                              gtcccAACCAGTCCCCCCCCCGACCCAGGACGGGATGGTCGCGGGAAACAGAATAATGCTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg                                                                           <  2:564613/147‑1 (MQ=255)                                                               tcccACCCAGTCCCCCCCCCGAACCATGCCGGGATGGTCGCGGGAAAAAAAATAATACTGGAAGCGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                          <  2:684121/148‑1 (MQ=255)                                                                tCCAACCAGTCCCCCCCCAAAACCATGACGCTTGGTTCGCCGGAAACAGAAAAATACTGAAAGGGAAGTTTGCGTGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                          <  2:307545/149‑1 (MQ=255)                                                                 caaaCCAGTCCCCCCCCCGAACCATGACGGGAGGGTCGCGGGGAAAAAATAAATACTGGAAGGGAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                       <  1:180495/149‑1 (MQ=255)                                                                     ccAGTACCGCCCCCGACCCATGACGCGATGGTCGCGGGAAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                   <  2:632157/151‑1 (MQ=255)                                                                              cccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                           >  2:225615/1‑150 (MQ=255)                                                                              cccccGAAAAATGACGGGATGGTCGCCGGAAAAAAAAAAATATTGGAAGAAAAGATTCCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa                                                          <  2:226974/151‑1 (MQ=255)                                                                                   gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                      >  1:896071/1‑150 (MQ=255)                                                                                    acccATGACGCGATGGTCGCGGGAAACAAAAAAATACTGGAAGCAAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                      <  2:465419/147‑1 (MQ=255)                                                                                      ccATGACGGGATGGTCGCCGGGAAAAAAATAATACTGGAAGCGAAGATGGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                   <  1:233260/150‑1 (MQ=255)                                                                                      ccATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                  >  1:786795/1‑151 (MQ=255)                                                                                       cATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgc                                                 >  1:347266/1‑151 (MQ=255)                                                                                        aTGCCGCGATGGCCGCGGGGAACAAAATAATACTGGAAGCAAAGTTTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                <  1:222033/151‑1 (MQ=255)                                                                                         tGACGCGATGGTCGCCGGAAAAAAAATAATACTGAAAGCGAAATTTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAGTGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                   <  1:592651/147‑1 (MQ=255)                                                                                                aTGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                        >  1:9168/1‑151 (MQ=255)                                                                                                 tGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc                                       >  1:823198/1‑151 (MQ=255)                                                                                                  ggTCGCGGGAAAAAAAATAATACTGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc                                       <  2:432755/150‑1 (MQ=255)                                                                                                      gCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                  >  2:389692/1‑151 (MQ=255)                                                                                                      gCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                  >  1:617542/1‑151 (MQ=255)                                                                                                       ccGGAAAAAAAAAAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                  <  2:946726/150‑1 (MQ=255)                                                                                                         ggaaaaaaaaaaataGTGGAAGCGAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                     <  1:933558/134‑1 (MQ=255)                                                                                                                 aataatACTGGAAGCGAAGATTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                       <  2:907232/151‑1 (MQ=255)                                                                                                                    aataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                     >  2:346829/1‑150 (MQ=255)                                                                                                                       acggaaaGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg                 <  2:338466/146‑1 (MQ=255)                                                                                                                          ggAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc              >  1:319369/1‑151 (MQ=255)                                                                                                                            aaGAGAAGATTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt            <  2:26071/151‑1 (MQ=255)                                                                                                                             aGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc           >  2:847739/1‑151 (MQ=255)                                                                                                                                gAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt        >  1:279498/1‑151 (MQ=255)                                                                                                                                 aaGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg       >  1:542502/1‑151 (MQ=255)                                                                                                                                  aGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc      >  2:492105/1‑151 (MQ=255)                                                                                                                                     tggCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa     <  1:316345/147‑1 (MQ=255)                                                                                                                                       gCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat  >  2:631569/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                             CATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAAT  >  NC_000913/3443134‑3443417

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: ‑