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

ACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATAC  >  NC_000913/3443131‑3443425                                                                                                                                                    |                                                                                                                                                     aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTTGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTAcccc                                                                                                                                                  >  2:217145/1‑151 (MQ=255)       agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCat                                                                                                                                            >  2:771536/1‑151 (MQ=255)             acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCAc                                                                                                                                      >  2:384867/1‑151 (MQ=255)                   aaaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTc                                                                                                                                  >  2:23715/1‑149 (MQ=255)                    aaGGGGTTGCCCGGGCTGCAAATACAGAAAAATTGTTGTCAGCCGGTTCCAACCAGTCCCGCCCCCGACCCATGACGCGATGGTCGCCGGGAACAAAATAATACTGAAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGc                                                                                                                               <  1:153494/151‑1 (MQ=255)                      gCGTTTGCCCGGGCGAAAAATACGGAAAATTTTTTGTCAGCCATTCCCAACCGGCCCCGCCCCCGACCCATGCCGCGAGGGTCCCCGGAAACAAAATAATCTTGGAAGC‑GAAGTTTGCGGGATTTCCCCCCGCCATATTCACTTTCAGCg                                                                                                                              <  2:636594/150‑1 (MQ=255)                           tGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                        >  1:911668/1‑151 (MQ=255)                             cccAGGCTAAAAATACAGCAAATTTGTTGTCCCCCAGTTCCACCCAGCCCCGCCCCCGACCCATGCCGCGATGGTCGCCGGAAACAAAATAATACTGGAAGC‑GAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                       <  1:160153/150‑1 (MQ=255)                                aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                   >  2:667734/1‑151 (MQ=255)                                       aaaaACAGCAAAGTTGTTGTCACCCGTTCCCAACCAGTCCCCCCCCAAAACCATGACGGGAGGGCCGCCGGAAACAAAATAATACTGGAGGC‑GAATATTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                              <  2:473330/149‑1 (MQ=255)                                       aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                            >  1:260022/1‑151 (MQ=255)                                                       ttgcCAGCCAGTCCCACCCAGTCCCCCCCCCGAACCATGACGCGATGGCCGCCGGAAACAAAATAATACTGAAAGC‑GAAGATTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                             <  1:869614/150‑1 (MQ=255)                                                        tgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                            >  2:320618/1‑150 (MQ=255)                                                          tCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc                                                                                         >  1:737026/1‑151 (MQ=255)                                                             gCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                       >  2:792384/1‑150 (MQ=255)                                                              ccAGTTCCACCCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGAAACAGAAAAATACTGGAAGC‑GAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc                                                                                         <  1:384875/147‑1 (MQ=255)                                                                  ttCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                                  >  1:207183/1‑150 (MQ=255)                                                                  tcccAACCAGCCCCGCCCCCAACCCATGCGGGGAGGGCCCCGGGAAAAAAAATAATACTGGAAGC‑GAGGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                        <  1:519576/142‑1 (MQ=255)                                                                    caacccAGTCCCCCCCCCGAACCAGGACGGGATGGTCGCGGGAAAAAAAATAATACGGGAAGC‑GAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                                    <  1:25661/142‑1 (MQ=255)                                                                        cctgtcccccccccAAACCATGACGCGATGGTCGCCGGAAAAAAAATAATACTGGAAGC‑GAAGATTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                                 <  2:635418/136‑1 (MQ=255)                                                                        ccggtcccccccccAAACCAGGACGGGAGGGTCGCCGGAAAAAAAAAAATCCGGGAAGC‑GAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                           <  2:630731/142‑1 (MQ=255)                                                                        ccagccccccccccAACCCATGACGCGAGGGCCGCCGGAAAAAAAAAAATACTGGAAGC‑GAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                            <  2:600643/141‑1 (MQ=255)                                                                              ccccccccAACCCATGCCGGGTGGGTCGCGGGGAAAAGAATAATACTGGAGGC‑GAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAg                                                                      <  2:370824/147‑1 (MQ=255)                                                                               cGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                                    >  2:784587/1‑151 (MQ=255)                                                                                gCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                                   >  2:650087/1‑151 (MQ=255)                                                                                   cccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc                                                                >  1:745453/1‑151 (MQ=255)                                                                                       aaCCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                             >  2:468257/1‑150 (MQ=255)                                                                                         caaTGAGGGGATGGTCGCCGGAAAAAAAATAATCCTGGAGGC‑GAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                          <  1:330034/149‑1 (MQ=255)                                                                                            ggccgggtgggccgcgGGAAAAAGAAAAATACGGGAAGC‑AAAGTTTGCCGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                           <  1:585987/135‑1 (MQ=255)                                                                                                gcgATGGTCGCGGGAAACAAAATAATACTGAAAGC‑GAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                   <  1:61766/151‑1 (MQ=255)                                                                                                gcgATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                   >  2:185788/1‑151 (MQ=255)                                                                                                gcgATGGCCGCCGGAAAAAAAAAAATACTGGAAGC‑AAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                   <  1:697833/151‑1 (MQ=255)                                                                                                  gATGGCCGCCGAGAACAAAAAAAT‑CTGGGAGCGGAAGATGGCGGTTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                   <  1:314983/149‑1 (MQ=255)                                                                                                     ggTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc                                               >  2:349943/1‑150 (MQ=255)                                                                                                          cgggaaaaaaaaaaataCGGGAAGC‑AAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                            <  1:289730/135‑1 (MQ=255)                                                                                                           cGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                         >  2:176763/1‑150 (MQ=255)                                                                                                            gggAAAAAAAAAATACGGGAAGC‑GAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaa                                       <  1:418721/151‑1 (MQ=255)                                                                                                                         taCTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                           >  1:149268/1‑150 (MQ=255)                                                                                                                          acggggaGC‑AAAGATGGGGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg                         <  2:494203/145‑1 (MQ=255)                                                                                                                            tGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                        >  1:87576/1‑150 (MQ=255)                                                                                                                             ggAAGC‑AAATATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                      <  2:622753/151‑1 (MQ=255)                                                                                                                             gaaaGC‑AAAGATTGCGGGAATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                        <  2:93733/147‑1 (MQ=255)                                                                                                                               aaGA‑GAAGATGGCGGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                    <  2:357121/151‑1 (MQ=255)                                                                                                                                     aaGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg               >  2:41983/1‑151 (MQ=255)                                                                                                                                       gATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc              >  2:39004/1‑150 (MQ=255)                                                                                                                                           gCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat          >  2:201079/1‑150 (MQ=255)                                                                                                                                             cGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact       >  2:238776/1‑151 (MQ=255)                                                                                                                                              gggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact       >  2:597019/1‑150 (MQ=255)                                                                                                                                                  ttACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatac  >  2:92421/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGC‑GAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATAC  >  NC_000913/3443131‑3443425

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