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

CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACT  >  NC_000913/3443128‑3443426                                                                                                                                                        |                                                                                                                                                      cAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACAATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTAc                                                                                                                                                      >  2:622643/1‑151 (MQ=255)  aGACGCATAGAGTAACACATAAAGCGGTTGCCAAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTAcc                                                                                                                                                     <  2:763987/151‑1 (MQ=255)        aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGcc                                                                                                                                               >  1:636489/1‑151 (MQ=255)        aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGcc                                                                                                                                               <  2:799331/151‑1 (MQ=255)        aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGcc                                                                                                                                               >  2:436128/1‑151 (MQ=255)            agTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCatat                                                                                                                                           >  2:746243/1‑151 (MQ=255)                           gtttGCCCAGCCGGAAAACACGAGGAAATTTTGTTCACCCAGTCCCAACCAGTCCCCCCCCGA‑AACCAGGCCGGGATGGCCCCCGAAACAAAAATAATACT‑GAGGGGGAGATTGGGGGAGATTCCCCCCGCCATATTCACTTTCAGCGGt                                                                                                                             <  2:129952/148‑1 (MQ=255)                                    ggctgaaaaacaaGCGAATTTGTTGTCCGCAATTCCCAACCAGCCCCCCCCCCA‑ACCCATGACGCGATGGTCGCGGGGAACAAAATAATTCTGGAAGAGAAATTTGCGGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc                                                                                                                   <  2:840748/139‑1 (MQ=255)                                       tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                >  2:403947/1‑151 (MQ=255)                                       tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                >  2:44568/1‑151 (MQ=255)                                       tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                >  1:573710/1‑151 (MQ=255)                                         cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                              >  2:393282/1‑151 (MQ=255)                                          aaaGACAGCGAAATTTTTCTCAGCCAGTTCCAACCAGTCCCGCCCCCA‑AACCATGACGCGGTGGCCGCCGGGACCAGAAAAATCCTGGAGGCAAAGATTGCCGG‑ATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                <  1:422473/148‑1 (MQ=255)                                               cAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                         >  1:542810/1‑150 (MQ=255)                                                   gAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                       >  2:335304/1‑148 (MQ=255)                                                       tttttgtCAGCCATTCCCACCCAGTCCCCCCCCCA‑AACCATGGCGCGATGGCCGCGGGGAACAAAATAATACGGGAAGCGAAGTTGGCGGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                                 <  1:75265/147‑1 (MQ=255)                                                           tgtCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                             >  2:169107/1‑150 (MQ=255)                                                              cAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                         >  2:92991/1‑151 (MQ=255)                                                               aGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                        >  1:510306/1‑151 (MQ=255)                                                                       ccAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                                >  1:422366/1‑151 (MQ=255)                                                                       caacccagccccccccccA‑ACCCAGGACGGGTTGGTCGCGGGAAAAAAAAAAATACTGGAAGAAAAGTTTGCGGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACg                                                                                 <  1:114682/137‑1 (MQ=255)                                                                           ccagccccGCCCCCA‑ACCCATGACGCGTTGGTCGCCGGGAAAAAAATAATGCTGGAGGCGAAGATGGCCGG‑GTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                                  <  1:20041/139‑1 (MQ=255)                                                                            cAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgc                                                                           >  1:13337/1‑151 (MQ=255)                                                                                 ccccccccGAAACCATGACGGGTGGGTCGGGGGAAAAAAAAAAAAATTGGAAGCGAAATTTGCCGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                               <  2:221332/140‑1 (MQ=255)                                                                                     ccccG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAc                                                                  >  1:329875/1‑151 (MQ=255)                                                                                      cccG‑ACCCAGGACGGTGTTGTCGCCGGCAAAAAAATAATCCGGGAAGCGAAGATTGCGTG‑TATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc                                                                 <  2:339504/151‑1 (MQ=255)                                                                                                tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg                                                        >  1:309313/1‑151 (MQ=255)                                                                                                 gACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                       >  1:86998/1‑151 (MQ=255)                                                                                                    gggttGGTCGCGGGAAAAAAAATAATACTGAAAGCGAAGATTGCGGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgc                                                          <  1:771654/141‑1 (MQ=255)                                                                                                    gcgATGGTCGCGGGGAACAAAATAATACGGGAAGCGAAGTTTGCCGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg                                                        <  1:746433/147‑1 (MQ=255)                                                                                                         ggTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                                 >  2:47683/1‑149 (MQ=255)                                                                                                          gTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                              >  1:606370/1‑151 (MQ=255)                                                                                                           tCGCCGGGAAAAAAAAAATACGGGAAGAAAAGTTTGCCGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                              <  1:62327/150‑1 (MQ=255)                                                                                                              ccGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                          >  2:349169/1‑151 (MQ=255)                                                                                                                   aaCAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                                      >  1:595183/1‑150 (MQ=255)                                                                                                                    aCAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTGTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                                     >  1:783579/1‑150 (MQ=255)                                                                                                                      aaaataatCCTGGAAGCGAAGTTTGCCGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACaaa                                    <  1:385034/147‑1 (MQ=255)                                                                                                                        aataatCCGGAAAGCAAAGTTTGCCGG‑GATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                 <  1:816000/150‑1 (MQ=255)                                                                                                                           aataGTGGAAGCGAAGTTTGCCGG‑GATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                              <  1:648851/150‑1 (MQ=255)                                                                                                                             taCTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                            >  1:356087/1‑150 (MQ=255)                                                                                                                                gggaagtgagtttggCCGG‑GACCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa                        <  2:274395/137‑1 (MQ=255)                                                                                                                                 ggAAGAAAATATTGCCGG‑GATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                       <  1:222313/151‑1 (MQ=255)                                                                                                                                     gagagggtcgCCGG‑GATTCCCCCCGCCACATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                      <  1:833494/139‑1 (MQ=255)                                                                                                                                            tggCCGG‑ATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc               <  2:541837/146‑1 (MQ=255)                                                                                                                                                cGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact        >  1:721673/1‑151 (MQ=255)                                                                                                                                                 gg‑gATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatacta       <  1:154947/151‑1 (MQ=255)                                                                                                                                                    gATTACTCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGACCGTAACTGCAaatacta       >  2:272440/1‑149 (MQ=255)                                                                                                                                                       tACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatact  >  1:758714/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGG‑GATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACT  >  NC_000913/3443128‑3443426

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