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

CGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT  >  NC_000913/3443132‑3443423                                                                                                                                                  |                                                                                                                                                   cGCATAGAGTAACACATAAAGCGGTTGCCCATGCTGCAAATGCAGCGAAATTGGTGTCAGCCAGTTCCAACCAGTACCGCCCCCGACTCATGAGGCGATGGTCGCCGGGAACAGAATAATTCTGGAAGCGAAGATTGCCGGGATTAccccc                                                                                                                                               <  2:753071/151‑1 (MQ=255)   cATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGc                                                                                                                                             <  2:768359/151‑1 (MQ=255)    aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                            >  1:612131/1‑151 (MQ=255)    aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                            >  1:388404/1‑151 (MQ=255)                  aagggggTTCCCCAGGCTGAAAATACAGGGAATTTTTTGCCACCTATTCCCAACCAGTCCCCCCCCCAACCCAGGACGGGTTGGCCGCCGGAAAAAAAATAATACTGGAGGCGAAATTGGCCGGAATTCCCCCCGCCATATTCACTTTCa                                                                                                                               <  2:343774/145‑1 (MQ=255)                      cGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGt                                                                                                                          >  1:411557/1‑151 (MQ=255)                       gtttGCCCAGGCTAAAAATACAGCAAATTTGTTGTCAGCCATTCCCAACCAGTCCCCCCCCCAAACCATGCCGCGATGGTCGCCGGAAACAAAATAATCCTGGAAGCAAAGATTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGt                                                                                                                          <  1:539646/148‑1 (MQ=255)                           gCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                     >  2:352792/1‑151 (MQ=255)                               aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                 >  1:675906/1‑151 (MQ=255)                                 gCTGAAAATACAGCAAAATTGTTGTCAGCCATTTCCACCCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTATGct                                                                                                               <  1:761530/151‑1 (MQ=255)                                   tgaaaaTACAGCAAATTTGTTGTCAGCCGGTCCCAACCAGTCCCGCCCCCGAACCATGCCGCGATGGCCGCCGGGAAAAAAATAATACTGGAAGCGAAGATTGCGGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                               <  1:314169/146‑1 (MQ=255)                                   tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                             >  1:180294/1‑151 (MQ=255)                                   tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                             >  1:47577/1‑151 (MQ=255)                                      aaaaACGGCGAAGTTTTGGTCAGCAAGTACAACCAAGCCCCCCCCCCAACCCATGACGCGTTGGTGGCGGGGAACAGAATAATACGGGAAGGGAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                           <  2:808769/150‑1 (MQ=255)                                       aatacgggaaattttttgtCGCCCATTTCCACCCATTCCCGCCCCCAAACCATGACGGGATGGTCCCCGGAACAAAAATATTACGGGAGCCGAAGTTTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                         <  2:148767/136‑1 (MQ=255)                                             gCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                    >  1:918337/1‑150 (MQ=255)                                              cGAAATGTTTGACACCCTGTTCCAACCGGACCCGCCCCCCGTCCACGACGCAAGGGCCGCCGGGAAAATAATAATACTGAAGGCAAAGATTGCCGGGAAACCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATa                                                                                                   <  2:545764/150‑1 (MQ=255)                                               gAAATGGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGAAAATGTGCGCTCTGTgcagca                                                                                                     >  2:915125/1‑147 (MQ=255)                                                      ttgtCAGCCAGTTCCACCCAGTCCCCCCCCCGAACCAGGCCGCGATGGTCGCCGGAAAAAAAATAATACTGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacga                                                                                            <  1:312677/149‑1 (MQ=255)                                                      ctgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGCGCTCTGTGCAGCATAGACacgacg                                                                                          >  2:649066/2‑151 (MQ=255)                                                          cAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                      >  1:664273/1‑151 (MQ=255)                                                                     acccAGTCCCCCCCCCAAACCATGAGGGGATGGTCGCGGGAAAAAAAAAAATCCTGGAAGCAAAGTTTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                             <  2:899171/147‑1 (MQ=255)                                                                         tgccccGCCCCCGACCCGTGACGCGAGGGCCGAGGGAAAAAAAATAATACTGGAAGGAAAAATGGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                         <  1:410774/145‑1 (MQ=255)                                                                           tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                      >  1:72549/1‑150 (MQ=255)                                                                            aCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                      >  1:870473/1‑149 (MQ=255)                                                                               gCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                                 >  2:667854/1‑151 (MQ=255)                                                                                     gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                            >  2:871241/1‑150 (MQ=255)                                                                                      aaCCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                           >  1:89899/1‑150 (MQ=255)                                                                                        ccATGACGGGATGGTCGCCGGAAAAAAAATAATCCTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                        <  1:273375/151‑1 (MQ=255)                                                                                        ccAGGACGCGATGGTCGCCGGGAAAAGAATAATACTGAAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                        <  1:257744/151‑1 (MQ=255)                                                                                             aCGCGATGGTCGCCGGAAACAAAAAAATACTGAAAGCGAATATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                      <  1:957790/148‑1 (MQ=255)                                                                                              cgcgtgggTCGCCGGGAACAAAATAATACTGGAAGCGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                    <  1:848175/143‑1 (MQ=255)                                                                                                 gATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                >  2:526510/1‑150 (MQ=255)                                                                                                  aTGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                              >  1:23621/1‑151 (MQ=255)                                                                                                       cGCCGGAAACAAAATAATAGTGGAAGCAAAGTTGGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                          <  1:139066/150‑1 (MQ=255)                                                                                                         ccgaaaaCAAAATAATCCTGAAGGCGAATATTGCCGTGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCCCGc                                         <  2:592075/144‑1 (MQ=255)                                                                                                                 aGAATAATCCGGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                               <  2:634024/151‑1 (MQ=255)                                                                                                                     taataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                           >  1:733370/1‑151 (MQ=255)                                                                                                                        tacggaaaGCGAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                         <  1:369134/144‑1 (MQ=255)                                                                                                                         aCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                         >  1:594943/1‑149 (MQ=255)                                                                                                                            ggaaggaaAGAGTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa                     <  2:706197/150‑1 (MQ=255)                                                                                                                            ggAAGCGAAGTTTGCGGGCATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                      <  2:951040/149‑1 (MQ=255)                                                                                                                            ggAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg                       <  1:150498/148‑1 (MQ=255)                                                                                                                            gaaaGCAAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACaaa                        <  1:403224/145‑1 (MQ=255)                                                                                                                              aggCGAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                   <  2:55019/148‑1 (MQ=255)                                                                                                                              aaGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                  >  1:629514/1‑151 (MQ=255)                                                                                                                                   aaaaTTGCCGGGTTGCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg             <  2:388224/148‑1 (MQ=255)                                                                                                                                   aaGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg             >  2:611936/1‑151 (MQ=255)                                                                                                                                      aTTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa           >  2:839090/1‑150 (MQ=255)                                                                                                                                       ttcccGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa           <  2:503959/146‑1 (MQ=255)                                                                                                                                         gCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa           <  2:787494/147‑1 (MQ=255)                                                                                                                                            ggtaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac      <  1:267976/146‑1 (MQ=255)                                                                                                                                            ggcaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat        <  2:702589/144‑1 (MQ=255)                                                                                                                                               aTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGACGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat  >  1:568503/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT  >  NC_000913/3443132‑3443423

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