Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	52.2%	58.7 / 77.1	46	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (21/1);  new base C (0/24);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 3.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-09

CGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATAC  >  NC_000913/3443132‑3443419                                                                                                                                                  |                                                                                                                                               cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                            >  2:410084/1‑150 (MQ=255) cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                            >  2:242105/1‑150 (MQ=255)     tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCa                                                                                                                                       >  1:97900/1‑151 (MQ=255)       gagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCACCata                                                                                                                                     >  2:140782/1‑151 (MQ=255)         gTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACTAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATAtt                                                                                                                                   >  2:157305/1‑151 (MQ=255)          aaaCACAAAAGGCGTTTGCCCGGGCTGAAAATACAGCGAAATTGTTGTCACCCAGTCCCAACCAGTCCCCCCCCCGACCCATGACGCGTTGGTCGCCGGGAACAGAATAATACGGGAAGCGAAGATGGCGGGGATCCCCCCCGCCATATTc                                                                                                                                  <  1:207947/150‑1 (MQ=255)                          tGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                  >  1:471484/1‑151 (MQ=255)                            ccgaGGCTGCAAAAACAGCAAAATTGTTGTCAGCCAGTACAAACCAGTCCCGCCCCCGACCAAGGACGCGAGGGTCCCCGGGAACAAAATAATACTGGAAGGGAAGATCGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                 <  1:349105/147‑1 (MQ=255)                             ccgggctcaaaaTACAGCAAAATTGTGGCCGCCCAGTCCCAACCAGTCCCCCCCCCAACCCAGGCCGCGATGGCCGCGGGGAACAAAATAATTCTGGAGGGAAAGTTTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgtg                                                                                                               <  1:100712/142‑1 (MQ=255)                             ccAGGCTGAAAATACAGCAAATTTGTTGTCAGCCAGTCCCAACCAGTCCCCCCCCCGACCCAGGACGCGGTGGCCGCCGGGAAAAAAATAATACTGGAAGCGAAATTTGCGGGGTTTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg                                                                                                               <  1:160427/151‑1 (MQ=255)                                  cTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                          <  2:146869/151‑1 (MQ=255)                                  cTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                           >  1:212730/1‑150 (MQ=255)                                  cTGAAAATACAGCGAAATTTTTGTCAGCCATTTCCAACCAGTCCCCCCCCCGAACCATGCCGCGATGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                          <  2:221768/151‑1 (MQ=255)                                              caaaaTTGTTGTCAGCCAGTCCCAACCATTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGGAACAAAATAATACTGAAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                <  2:31881/147‑1 (MQ=255)                                              caaaaTTGTTGCCAGCCATTCCAAACAAGTCCCCCCCCCGACCCATGACGCGATGCCCCCCGGAAACAAAAAAATACGGAAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                <  2:404531/147‑1 (MQ=255)                                              cGAAATTGTTGTCAGCCATTTCCAACCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                <  2:386660/149‑1 (MQ=255)                                              aaaatttgtttCCAGCCATTTCAAACCATCCCCCCCCCCAAACTATGCCGCGATGCTCGCCGGGAACAAAAAAATTCTGAAAGCAAAGTGTGCGGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                <  1:294345/144‑1 (MQ=255)                                                  tttgttgTCAGCCAGTTCCAACCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                          <  2:342086/150‑1 (MQ=255)                                                         tCAGCCAGTTCCACCCAGTCCCGCCCCCAAACCATGACGCGATGGTCGCCGGAAACAAAAAAATACTGAAAGCAAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac                                                                                    <  2:20505/150‑1 (MQ=255)                                                           aGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                 >  2:112135/1‑151 (MQ=255)                                                             ccAGTCCCAACCAGTCCCGCCCCCGAACCAGGACGGGAGGGCCGCGGGAAAAAAAATAATACTGGAAGCGAAGTTTGCCGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                               <  1:453244/151‑1 (MQ=255)                                                                gTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                            >  2:405306/1‑151 (MQ=255)                                                                     aaCCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                       <  2:175599/151‑1 (MQ=255)                                                                     aaCCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                       >  2:99319/1‑151 (MQ=255)                                                                          gTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                  >  2:138642/1‑151 (MQ=255)                                                                          gTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                  >  1:190684/1‑151 (MQ=255)                                                                          gTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                  >  1:286548/1‑151 (MQ=255)                                                                                cccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa                                                            >  1:418717/1‑151 (MQ=255)                                                                                    cGACCCATGACGCGATGGTCGCCGGAAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                        <  2:466786/151‑1 (MQ=255)                                                                                    cGAACCATGACGGGTTGGCCGCCGGAAAAAAAATAATCCTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                        <  2:471924/151‑1 (MQ=255)                                                                                    cGAACCAGGAGGGGATGGCCGCCGGAAAAAAAAAAATACGGGAAGGGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                        <  2:452935/151‑1 (MQ=255)                                                                                          aTGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                  >  1:149684/1‑151 (MQ=255)                                                                                             aCGCGACGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                               >  1:214062/1‑151 (MQ=255)                                                                                              cgcgATGGCCGCCGGAAAAAGAATAATATTGAAAGCGAATATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGt                                              <  2:140496/151‑1 (MQ=255)                                                                                                        cccGGAAACAAAATAATCCTGGAAGCGAAGATTCCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                    <  2:89914/150‑1 (MQ=255)                                                                                                             gAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCGCCGCCATATTCACTTTCAGCGGTAAATGTGTGATCTGTGCAGCATAGACACGACCACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                               >  1:346862/1‑151 (MQ=255)                                                                                                             gAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                               >  1:430414/1‑151 (MQ=255)                                                                                                             gAACAAAAAAATAGGGGAAGGAAAGGTTGCCGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCACAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                <  1:356004/150‑1 (MQ=255)                                                                                                             aaaCAAAATAATCCGGGAAGCGAAGTTTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                               <  2:212734/150‑1 (MQ=255)                                                                                                              aaaaaaataatACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGGCCACGCTcaaca                               <  2:112549/145‑1 (MQ=255)                                                                                                                  gAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                           >  1:367866/1‑150 (MQ=255)                                                                                                                  gAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                           >  1:410828/1‑150 (MQ=255)                                                                                                                  gAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGGAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                           >  1:54130/1‑150 (MQ=255)                                                                                                                     taatcCTGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                       <  2:434577/151‑1 (MQ=255)                                                                                                                                      aTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa      <  2:470617/151‑1 (MQ=255)                                                                                                                                          cggggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac  <  2:108609/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               CGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATAC  >  NC_000913/3443132‑3443419

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