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.2%	80.9 / 20.0	51	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (29/2);  new base C (0/20);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 2.98e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCATAGAGTAACACATAAAGCGGTTGCCCAGG‑CTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443132‑3443427                                                                                                                                                       |                                                                                                                                                       cGCATAGAGTAACACATAAAGCGGTTGCCCAGG‑CTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTAccccc                                                                                                                                                   >  1:634525/1‑151 (MQ=255) cGCATAGAGTAACACATAAAGCGGTTGCCCAGG‑CTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTAccccc                                                                                                                                                   >  2:828854/1‑151 (MQ=255)          tAACACATAAAGCGGTTGCCCAGG‑CTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTc                                                                                                                                          >  2:169866/1‑151 (MQ=255)                  aaaGCGGTTGCCCAGG‑CTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTc                                                                                                                                    >  1:1001332/1‑149 (MQ=255)                          tttcccgggagtaaaaTACAGCAAAATTGTTGCCACCGAGTCCC‑AC‑CCG‑TTCCCGCCCCCGACCCAGGGCGGGATGGTCGCCGAAAACAGAAAAATATGGGAACGAAAGAT‑TCCCGGAATTCCCCCCGCCATATTCACTTTCAGCgg                                                                                                                               <  1:1206322/134‑1 (MQ=255)                             ccAGGCCTGAAAATCCAGCGAATTTTTTGTCGCCCTGTCCC‑AA‑CCA‑GCCCCCCCCCCGAACCATGACGGGATGGCCGCCGGAAAAAAAAAAATCCGGGAAGGGAAGAT‑GGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTaaa                                                                                                                           <  2:637451/148‑1 (MQ=255)                                   cTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                                  >  1:384506/1‑151 (MQ=255)                                     gCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                                >  2:1126813/1‑151 (MQ=255)                                      cAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                               >  1:386288/1‑151 (MQ=255)                                       aaaaACAGCAAATTTGTTGTCACCCA‑TTCC‑AACCCA‑GTCCCCCCCCGGAACCATGCCGGGATGGTCGCGGGAAAAAGAAAAATCCTGGAAGCAAAGTT‑TGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                                <  1:641189/149‑1 (MQ=255)                                       aaaTACAGCGAATTTTTTGTCAGCCAGTCCC‑AA‑CCA‑GTCCCGCCCCCAAACCATGACGCGATGGTCGCCGGAAACAAAATAATACAGGAAGCAAAGTT‑GGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                               <  1:1196115/150‑1 (MQ=255)                                        aaacgaGCGAATTGTTTGCCAGCAGGTCCC‑AA‑CCA‑GGGCCCCCCCCGACCGACGCGGTGAGGGCCGCCGGGAAAAGAATAATGGGGGAGGGGAAGTT‑TGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                             <  2:900335/146‑1 (MQ=255)                                                   aTTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                                  >  2:643333/1‑151 (MQ=255)                                                     ttgtgctggCCTGTTCC‑AA‑CCAGGCCCCCCCCCCGAATAAGGGCGGGGTGGCCGCCGGGAAGAAAAGAATAGAGGAGGGAAAGAT‑GGCGGGGTTGCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAg                                                                                                      <  2:710865/138‑1 (MQ=255)                                                       ttgcCACCCGGTCCC‑AA‑CCA‑TCCCCCCCCCCGACCCAGGCCGCGTTGGTCGCCGGGAAAAAAAAAATACTGGAAGCGAAGTT‑TGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacga                                                                                                <  1:1140186/149‑1 (MQ=255)                                                        tgggcgCCAGCTCCAAA‑CCA‑GACCCCCCCCCAAACCATGGCGGGGTGCTCGCGGGAAAAAGAAAAAAGCGGAAAGAGAAGAT‑GGCGGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                               <  2:1007976/145‑1 (MQ=255)                                                            aGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                         >  1:1307080/1‑151 (MQ=255)                                                              ccattccc‑AC‑CCA‑GTCCCCCCCCCGAACCATGACGCGAGGGTCGCCGGAAAAAAAATAATACGGGAAGCGAAGAT‑TGCCGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                                       <  1:934920/145‑1 (MQ=255)                                                              ccAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                                       >  1:925550/1‑151 (MQ=255)                                                              ccAGTACA‑AA‑CCA‑GTCCCCCCCCCGAACCATGACGCGAGGGTCGCCGGAAACAAAATAATACGGGAAGCGAATAT‑TGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                         <  1:359497/149‑1 (MQ=255)                                                                        acccA‑GTCCCCCCCCCGACCCATGACGCGAGGGTCGCGGGAAAAAAAAAAATACTGGAAGCAAAGAT‑TGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                                   <  1:209627/145‑1 (MQ=255)                                                                        a‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                              >  1:585601/1‑151 (MQ=255)                                                                            a‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                                           >  2:27822/1‑151 (MQ=255)                                                                                  cGCCCCCGAACCATGACGCGATGGTCGCGGGAAACAAAATAATACTGAAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATACTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtggcgaggaggt                                                                      <  1:1333174/151‑12 (MQ=255)                                                                                     ccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAc                                                                   >  1:483053/1‑151 (MQ=255)                                                                                        cgagccagcaCGCGAGGGCCGCGGGAAAAAAAATAATACTGAAAG‑GTAGATGTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACca                                                                 <  1:169872/141‑1 (MQ=255)                                                                                              agagagggATGGGCGCCGGAAAAAAAATAAGCTTGGAGCGGATGTT‑TGCCGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                             <  2:970784/141‑1 (MQ=255)                                                                                               tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg                                                         >  2:714887/1‑151 (MQ=255)                                                                                                           cGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                             >  2:619988/1‑151 (MQ=255)                                                                                                              cGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                           >  2:873661/1‑150 (MQ=255)                                                                                                                ggAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                        >  2:284980/1‑151 (MQ=255)                                                                                                                gaaaaaaaaaaataCTGAAAGCGAAGAT‑GGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                           <  2:1121508/138‑1 (MQ=255)                                                                                                                  aaCAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                                       >  1:1209021/1‑150 (MQ=255)                                                                                                                  aaCAAAAAAATACTGGAAGCAAAGTT‑GGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                                      <  1:469511/151‑1 (MQ=255)                                                                                                                      gAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                  >  1:133206/1‑151 (MQ=255)                                                                                                                      aaaaacaattGGAAGCAAGGAT‑TGCCGCATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACaaa                                     <  1:480815/139‑1 (MQ=255)                                                                                                                      aaaaaatcCTGGAAGCGAAGTT‑TGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                  <  2:716358/147‑1 (MQ=255)                                                                                                                        ataatcCTGGAAGCGAATAT‑TGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaa                                <  1:390588/151‑1 (MQ=255)                                                                                                                           ataCTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                              >  1:514268/1‑150 (MQ=255)                                                                                                                              cTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             >  2:20444/1‑148 (MQ=255)                                                                                                                                      gAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  >  1:461769/1‑151 (MQ=255)                                                                                                                                      gAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  >  2:160765/1‑151 (MQ=255)                                                                                                                                      gAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  >  2:1312481/1‑151 (MQ=255)                                                                                                                                      gAAGAT‑TGCCGGGATTACCACCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCCTAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  >  2:912836/1‑151 (MQ=255)                                                                                                                                       aaGAT‑TGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGTGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                     <  1:424780/147‑1 (MQ=255)                                                                                                                                        aGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc                >  2:120871/1‑151 (MQ=255)                                                                                                                                         gtt‑tGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg                 <  1:1073531/147‑1 (MQ=255)                                                                                                                                           t‑tGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa              >  1:1154755/1‑150 (MQ=255)                                                                                                                                           t‑tGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa              >  2:1129129/1‑150 (MQ=255)                                                                                                                                                 gcaaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  <  2:534454/138‑1 (MQ=255)                                                                                                                                                      tACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatact   >  2:586190/1‑151 (MQ=255)                                                                                                                                                       aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTg  >  2:843035/1‑151 (MQ=255)                                                                                                                                                       aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTg  >  1:827343/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CGCATAGAGTAACACATAAAGCGGTTGCCCAGG‑CTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCC‑AA‑CCA‑GTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGAT‑TGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443132‑3443427

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