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.4%	89.3 / 13.6	54	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (26/4);  new base C (0/24);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAA  >  NC_000913/3443127‑3443422                                                                                                                                                        |                                                                                                                                                  gCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTa                                                                                                                                                   <  2:198446/151‑1 (MQ=255)   aGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcc                                                                                                                                                 >  1:717731/1‑151 (MQ=255)    gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAccc                                                                                                                                                >  1:122661/1‑151 (MQ=255)    gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAccc                                                                                                                                                >  1:172897/1‑151 (MQ=255)     aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTACCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                               <  1:657305/151‑1 (MQ=255)      cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAccccc                                                                                                                                              <  1:875331/151‑1 (MQ=255)      cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAccccc                                                                                                                                              <  1:537741/151‑1 (MQ=255)          tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGc                                                                                                                                            >  2:773932/1‑149 (MQ=255)           agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCat                                                                                                                                         >  1:106422/1‑151 (MQ=255)                  cacaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACt                                                                                                                                  >  1:325451/1‑151 (MQ=255)                  cacaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACt                                                                                                                                  >  1:480125/1‑151 (MQ=255)                       aaaGCGGTTGCCCAGGCTGAAAATACAGCGAATTTGTTGTCGGCCATTCCCAACCAGTCCCGCCCCCGACCCATGA‑CGCGATGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACtt                                                                                                                                 <  1:387932/147‑1 (MQ=255)                        aaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGc                                                                                                                            >  1:506060/1‑151 (MQ=255)                                  ccgggccacAAATAAAGAAAAATTGTTGTCACCCAG‑TCCCACCAACTCCCCCCCCCAACCAAGACCGGGATGGTCCCGGGAAACAAAATAATCCTGGAAGGAAAGTTTGCGGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                                   <  2:454421/142‑1 (MQ=255)                                    aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                >  1:254857/1‑151 (MQ=255)                                    aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                >  2:479802/1‑151 (MQ=255)                                     ggCTGAAAATACAGCAAATTTGTTGTCAGCCAGTCCCACCCAGTCCCCCCCCCGAACCATGA‑CGCGATGGTCGCCGGAAACAGAATAATACTGGAAGCGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                <  1:1000029/150‑1 (MQ=255)                                          aaaaTACAGAGAATTTGTTGCCAGCCAGTTCCAACCAGCCCCCCCCCCAACCCATGA‑CGAGAGGGCCGCCGCGAACAAAAGAATACTGGAGGCGAGGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                           <  2:12036/149‑1 (MQ=255)                                           aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                         >  1:999286/1‑151 (MQ=255)                                           aaaTACAGCAAATTTTTTGTCAGCCATTCCCACCCAGTCCCCCCCCCGACCCAGGC‑CGCGAGGGGCGCCGGGAACAAAAAAATACGGGAAGCAAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                           <  1:586447/149‑1 (MQ=255)                                               aCAGCGAAATGGTTGCCCGCCAGTACCACCCAGTCCCCCCCCCGACCCATGA‑CGGGATGGTCGCCGGAAACAAAATAATCCTGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcag                                                                                                      <  1:529178/150‑1 (MQ=255)                                               aCAGCAAATTTGTTGTCAGCCAGTCCCACCCAGTCCCGCCCCCGACCCATGA‑CGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATGGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcag                                                                                                      <  1:144891/150‑1 (MQ=255)                                                           ttgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                         >  1:365448/1‑151 (MQ=255)                                                               cAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                     >  2:367280/1‑151 (MQ=255)                                                                aGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                    >  2:684222/1‑151 (MQ=255)                                                                          aaCCAGCCCCCCCCCCAAACCATGA‑CGCGATGGTCGCCGGAAAAAAAAAAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACg                                                                             <  1:451483/148‑1 (MQ=255)                                                                           aCCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                         >  1:352728/1‑151 (MQ=255)                                                                            ccAGTCCCCCCCCCAACCCAGGA‑CGCGATGGCCGCGGGAAAAAAAAAAATGCGGGAAGCAAAAATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                        <  1:556685/151‑1 (MQ=255)                                                                              aGTCCCCCCCCCGAACCATGA‑CGGGTGGGTCGCCGGAAACAAAATAATCCTGGAAGCAAAGATGGCCGGGTTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                                 <  2:65015/140‑1 (MQ=255)                                                                              aGTACCGCCCCAGGACCATGA‑CGCGATGGTCGCCGGGGACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                                      >  2:363268/1‑151 (MQ=255)                                                                                ccccccccccGAACCATGC‑CGGGAGGTCCCCCGGAAACAAAAAAATACAGGAGAGGAAGTTGGCGGAAATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                                    <  2:44832/146‑1 (MQ=255)                                                                                 aCCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                     >  2:469491/1‑149 (MQ=255)                                                                                  ccccccccGACCCATGG‑CGCGAGGGGCCCGGGAAACAAAATAATACGGGAGGAAAAGATTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                                    <  1:340960/146‑1 (MQ=255)                                                                                         cGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                           >  2:552679/1‑151 (MQ=255)                                                                                           aaCCATGA‑GGCGAGGGCCGCCGGGAAAAAAATAATACTGGAAGCGAAGATTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACca                                                            <  2:979033/148‑1 (MQ=255)                                                                                                       gagggTCGCGGGGAACAAAATAATGTTGGAAGCGAGGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                               <  1:336337/147‑1 (MQ=255)                                                                                                          ggTCGCGGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                     <  1:78961/141‑1 (MQ=255)                                                                                                          ggTCGCCGGAAAAAGAATAATACTGGAAGCAAAGTTTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCa                                           <  2:222838/151‑1 (MQ=255)                                                                                                                 gggAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                                     >  1:589686/1‑150 (MQ=255)                                                                                                                 ggaaaCAAAATATTACTGGAAGCGAGGTTTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                          <  1:791496/142‑1 (MQ=255)                                                                                                                 gaaaaCAAAATAATACTGGAGGGGAAGATGGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                        <  2:505983/144‑1 (MQ=255)                                                                                                                   gAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                   >  1:361105/1‑150 (MQ=255)                                                                                                                    aaCAGAATAATACTGGAAGCGAAGATTGCCGGGATGACCCTCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCGGCATAGACACGACGGCCTTGCTGACGTTTCGCGTAGTTTACCACAATTCGGCGATGACCACGCTcaa                                    >  2:773099/1‑148 (MQ=255)                                                                                                                       aaaataatACTGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                              <  1:328744/149‑1 (MQ=255)                                                                                                                       aaaaaaataCGGGAAGGGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                              <  1:479877/146‑1 (MQ=255)                                                                                                                         aaaaataCGGAAAGGGAAGTTTGCCGGATTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                            <  2:105729/148‑1 (MQ=255)                                                                                                                          ataataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                            >  1:59872/1‑150 (MQ=255)                                                                                                                          ataataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                            >  2:697052/1‑150 (MQ=255)                                                                                                                            aatcCTGGAAGCGAAGATGGCGGGGTTACCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                          <  2:899075/150‑1 (MQ=255)                                                                                                                               aCTGGAAGCGAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg                      <  1:231268/151‑1 (MQ=255)                                                                                                                                cTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                        >  1:607209/1‑148 (MQ=255)                                                                                                                                       cGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCAc              >  2:423191/1‑151 (MQ=255)                                                                                                                                           gATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc           >  2:9840/1‑150 (MQ=255)                                                                                                                                                cggggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata      <  1:78858/148‑1 (MQ=255)                                                                                                                                                   ggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaa  <  1:430535/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  GCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAA  >  NC_000913/3443127‑3443422

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