Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	51.1%	53.9 / 71.3	45	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (20/2);  new base C (0/23);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 7.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-08

ACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443131‑3443427                                                                                                                                                         |                                                                                                                                                       aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTAcccc                                                                                                                                                    >  2:203805/1‑151 (MQ=255) aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTAcccc                                                                                                                                                    >  1:409204/1‑151 (MQ=255)  cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTAccccc                                                                                                                                                   <  1:70894/151‑1 (MQ=255)   gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCg                                                                                                                                                  >  2:634086/1‑151 (MQ=255)       agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCat                                                                                                                                              >  1:182787/1‑151 (MQ=255)         agTAACACATAAAGCGGTTGCCCGGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑ACATGGCGCGATGGTCGCCGGGAACAGAATAAGACTGGAAGCGAAAATTGCCGGG‑A‑TTACCCCCGCCatat                                                                                                                                            <  2:243650/151‑1 (MQ=255)           tAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTc                                                                                                                                          >  2:304177/1‑151 (MQ=255)                   acggCGGTTGCCCGGGCTCAAAATAAAGAAAATTTG‑TTGC‑CGGCATGTCCACCCC‑AGTCCCCCCCCCAAA‑CCTGGCGGCGAGGCCCGCCGGCAAAAAAAAAATACTGGAAGCGAAGTTTGCCGGT‑T‑TTCCCCCCGCCATATTCACTTTCAg                                                                                                                                  <  2:642846/148‑1 (MQ=255)                        gtttGCCGAGGCTGCAAATACAGCAAAATTG‑TTGT‑CAGCCAGTCCAAACC‑AGTCCCCCCCCGGAA‑CAAGGACGCGATGGTCGCGGGAAACAAAATAATATGGGAAGCGAAGTTTGCCGGG‑T‑TTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                             <  1:611388/149‑1 (MQ=255)                        ggTTGCCCGGGCTGAAAATACAGAGAAATTT‑TTGT‑CAGCCAGTACAAACC‑GGTCCCGCCCCCGAA‑CCATGACGGGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAGGATTGCGGGA‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                             <  1:172622/151‑1 (MQ=255)                             cccAGGCTGCAAATCCAGCAAATTTT‑TTGC‑CACCCATTTCCAACC‑AGTCCCCCCCCCGAA‑CCATGACGCGATGGTCGCGGGGAAAAAAATAATCCTGAAAGCGAAGATTGCCGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                         <  1:102541/150‑1 (MQ=255)                             cccAGGCTGCAAATACAG‑AAAATTT‑TTGTCCAGCCAGTTCCACCC‑AGTCCCGCCCCCGAC‑CCATGCCGCGAGGCCCGCCGGAAAAAAAATAATACTGGAAGCAAAGATTGCCGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGc                                                                                                                                 <  2:54589/142‑1 (MQ=255)                             cccAGGCTGAAAATACAGCAAATTTG‑TTGT‑CAGCCATTTCAAACC‑AGTCCCGCCCCCGAA‑CAATGCCGCGATGGTCGCCGGAAAAAAAATAATACTGGAAGCGAAGATTGCGGGA‑T‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                                        <  1:503040/151‑1 (MQ=255)                             cccAGGCTGAAAATACAGCAAATTTG‑TTGT‑CAGCCATTCCCACCC‑AGTCCCCCCCCCGAA‑CCATGCCGCGATGGCCGCCGGAAAAAGAATAATACTGGAAGCGAAGATTGCGGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                                        <  1:646379/151‑1 (MQ=255)                                       aaaaccagataaaTTG‑TGTC‑CAGCCTGTCCCACCG‑ATTACCGCCCCCGAA‑CAGGGAGGCGTGGGGCGCGGGAAACAGAAAAATTCGGGGAGCGAAGTTTGGCGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                 <  2:582538/138‑1 (MQ=255)                                       aaaaacaagaaaaTTT‑TTGC‑CCGCCCTTTCCACCC‑GGTCCCCCCCCCGAC‑CCATGCCGGGTTGGTCGCGGAAAAAAAAATAATGGTGGAGGGGAA‑ATTGGCGGGTT‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                 <  2:504929/138‑1 (MQ=255)                                        aaTACAGCAAATTTG‑TTGC‑CAGCCACTCCCAACC‑AGTACCCCCCCGGAA‑CAAGGACGGGATGGCCGCGGGAAAAAGAATAATACTGGAAGCGAAGATGGCGGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                             <  2:654175/151‑1 (MQ=255)                                                aaaaTTGTTTGT‑AACCCAG‑TCCAACCGAGTCCCCCCCCCGAA‑CGATGACGGGATGGTCGCGGGGAAAAAAAAAATAGGGGAGGGAAAGAGCGGGGGA‑A‑TCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgca                                                                                                            <  2:690922/144‑1 (MQ=255)                                                              aGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                         >  1:82708/1‑151 (MQ=255)                                                                     tCCAACC‑AGTACCGTCGACGAA‑CCATGACGCGGTGGCCGCTGGGCTCCGAATAATAATGGAAGAGCAGATTGCCGGG‑A‑TTACCCCCGCCATGTTAACTTTCCGCGGCAGCTGTGGTCTCTGTGCAGAATAGACCCGACGACCATGCTGAc                                                                                   >  2:468653/1‑150 (MQ=255)                                                                       cAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                                >  1:131898/1‑151 (MQ=255)                                                                       cAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                                >  1:472757/1‑151 (MQ=255)                                                                         aCC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                              >  1:608335/1‑151 (MQ=255)                                                                          cc‑AGTCCCCCCCCCAAA‑CCAGGACGGGATGGTCGCGGGAAACAGAATAATATTGGAAGCGAAGATTGCGGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                             <  1:359734/151‑1 (MQ=255)                                                                             ccgaCCCCCCCCAAACCCAGGACGGGTTGGGCGCGGGAAACAAATAAAGACGGGAAGAAAAG‑TTGGCGGT‑ATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                               <  2:680028/145‑1 (MQ=255)                                                                             aGTCCCCCCCCCAAA‑CCGGGCCGCGTGGGCCGCGGGAAACAAAAAAATACTGGAAGCAAAGATTCCGGGT‑T‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                                    <  2:38343/142‑1 (MQ=255)                                                                              gTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                          >  2:238505/1‑151 (MQ=255)                                                                                 ccccccccGAA‑CCATGACGGGGTGGTCGCGGGGAACAAAAAAATACTGAAAGGGAAGTTTGGGTGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                                           <  2:200655/144‑1 (MQ=255)                                                                                                    gcgATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                     >  2:455810/1‑151 (MQ=255)                                                                                                            cGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                             >  2:657716/1‑151 (MQ=255)                                                                                                               cGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                           >  1:552632/1‑150 (MQ=255)                                                                                                                 ggAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTACAGCATAGACACGACGACCTTGCTGACGTTTCGTGTAGTTCACCACAATGCGGCGTTGACCACGCTcaa                                         >  2:662481/1‑150 (MQ=255)                                                                                                                      aGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                    <  1:189205/150‑1 (MQ=255)                                                                                                                        aataatACTGGAAGCGAAGATTGCGGGG‑T‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                                 <  2:329099/151‑1 (MQ=255)                                                                                                                        aaaaataCGGGAAGCGAGAATGGCGGGT‑T‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                                 <  2:110524/148‑1 (MQ=255)                                                                                                                          taataCTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                               >  2:310572/1‑151 (MQ=255)                                                                                                                              acgggAAGCGAAGAGGGCCGGA‑T‑CTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACaaa                            <  2:564177/147‑1 (MQ=255)                                                                                                                                 ggAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                        >  1:198504/1‑151 (MQ=255)                                                                                                                                    aGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                     >  2:196835/1‑151 (MQ=255)                                                                                                                                      cGAAGATTGCCGGG‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCa                    <  1:618538/150‑1 (MQ=255)                                                                                                                                       gAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  >  1:377205/1‑151 (MQ=255)                                                                                                                                               cggtg‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata           <  2:217832/146‑1 (MQ=255)                                                                                                                                               ccGGA‑A‑TTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaaa             <  1:377067/148‑1 (MQ=255)                                                                                                                                                 ggg‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact         >  2:160890/1‑150 (MQ=255)                                                                                                                                                 ggg‑A‑TCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact         <  1:634314/150‑1 (MQ=255)                                                                                                                                                         aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTg  >  1:591644/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       ACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTG‑TTGT‑CAGCCAGTTCCAACC‑AGTACCGCCCCCGAA‑CCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGG‑A‑TTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443131‑3443427

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