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.7%	54.0 / 13.1	38	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (21/0);  new base C (1/16);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 7.64e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTA  >  NC_000913/3443144‑3443421                                                                                                                                        |                                                                                                                                                 cacaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACt                                                                                                                                 >  1:291658/1‑151 (MQ=255)   caTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACttt                                                                                                                               >  2:577448/1‑151 (MQ=255)     aaagggggTTGCCCGGGCTCAAAAAACAGAAAAATTGTTCCCGCCCATTCCCAACC‑AGTCCCCCCCCCAAACCATGACGCTGGGGTCGCGGGA‑AACAAAAAAATCCTGGAAGCAAAGTTGGCCGGAATTCCCCCCGCCATATTCACTTTCa                                                                                                                             <  2:664708/145‑1 (MQ=255)      agaggtgTTGCCCAGGCTGAAAATACAGCGAAATTTTTGGCACCCA‑TTCCAACCAAGGCCCCCCCCCGACCCAGGGCGGGGTGGTGGCCGGG‑AAAATAATAATTTTGGAGACGAGAATTGCGGGGATACCCCCCGCCATATTCACTTTCa                                                                                                                             <  2:106951/144‑1 (MQ=255)             ttGCCCAGGATGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTaa                                                                                                                      >  1:102255/1‑150 (MQ=255)              tGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGAACCGCCCCCGAACCATGACGCGATGGTAGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACACCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                    >  1:605477/1‑151 (MQ=255)               gCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                   >  1:590379/1‑151 (MQ=255)                      cTGCAAATACAGCGAATTTTTTGCCAGCCAGTCCCAACC‑AGTCCCGCCCCCGACCCAGGACGCGATGGT‑GCCGGGAAAAAAAATAATTCTGGAGGGGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                             <  2:3220/150‑1 (MQ=255)                          aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTCCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                        >  1:463494/1‑151 (MQ=255)                          aaaTACAGCGAAATTGTTGCCAGCCATTTCCACCC‑AGTCCCGCCCCCGAACCATGACGCGATGGTCGCGGGA‑AACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                        <  1:128912/151‑1 (MQ=255)                           aaTACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                         >  2:272366/1‑149 (MQ=255)                             tACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                       >  2:652340/1‑149 (MQ=255)                                    aaaTTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                              >  2:120274/1‑151 (MQ=255)                                       ttgtggTCAGCCAGTCCCACCC‑AGTACCCCCCCCAACCCAGGACGCGATGGTCGCCGGA‑AAAAAAATAATACGGGAAGCGAAGATTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                            <  1:125248/150‑1 (MQ=255)                                         tttgtCGGCCAGTCCAAACC‑AGTCCCGCCCCCAACCCATGGCGGGATGGTCGCGGGA‑AACAGAATAATACGGGAAGAAAAGATGGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacga                                                                                          <  1:523398/149‑1 (MQ=255)                                                      tCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACg                                                                            >  2:473510/1‑151 (MQ=255)                                                       ccAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                           >  2:436511/1‑151 (MQ=255)                                                                        ccGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc                                                            >  2:34137/1‑150 (MQ=255)                                                                                  aCGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                    >  2:665051/1‑148 (MQ=255)                                                                                    gcgATGGTCGCCGGA‑AAAAAAATAATACTGAAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                  <  1:493296/148‑1 (MQ=255)                                                                                      gATGGCCGCCGGA‑AACAAAATAATACGGAAAGCGAAGTTTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                  <  1:174957/146‑1 (MQ=255)                                                                                        gtgTCCCCGGGAAACAAAAAAATGCTGGAGGAGAAATTCGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                            <  1:265996/149‑1 (MQ=255)                                                                                         ggTCGCGGGA‑AAAAAAATAATACGGGAAGCGAAGTTTGCCGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCa                                          <  1:497689/151‑1 (MQ=255)                                                                                              ccGGA‑AACAAAATAATACGGAAAGCAAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                     <  1:199828/151‑1 (MQ=255)                                                                                                    aaCAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                  <  1:159828/149‑1 (MQ=255)                                                                                                    aaCAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                                 >  2:600616/1‑150 (MQ=255)                                                                                                      cAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                              >  1:463380/1‑151 (MQ=255)                                                                                                        gAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                            >  2:148077/1‑151 (MQ=255)                                                                                                        gAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                            >  2:142706/1‑151 (MQ=255)                                                                                                             aaacggaaaggaaaGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                         <  2:635353/137‑1 (MQ=255)                                                                                                                cTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                       >  2:603213/1‑148 (MQ=255)                                                                                                                cTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                       >  2:572028/1‑148 (MQ=255)                                                                                                                cTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGAGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                       >  1:318709/1‑148 (MQ=255)                                                                                                                        gAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt            >  2:555781/1‑151 (MQ=255)                                                                                                                        gAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt            >  2:216063/1‑151 (MQ=255)                                                                                                                             ttGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa         <  1:296581/149‑1 (MQ=255)                                                                                                                               gcgggtaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa         <  1:397163/141‑1 (MQ=255)                                                                                                                                  gggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatacta  <  2:336468/151‑1 (MQ=255)                                                                                                                                        |                                                                                                                                                 CACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACC‑AGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGG‑AACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTA  >  NC_000913/3443144‑3443421

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