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

CATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443134‑3443427                                                                                                                                                |                                                                                                                                                       cATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGc                                                                                                                                                 >  1:754939/1‑151 (MQ=255)  aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                                >  1:325174/1‑151 (MQ=255)   tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                                >  2:739225/1‑150 (MQ=255)         aaCACAAAAAGCGTTTGCCCGGGCTGAAAATACAGCGAAATTTTTGTCAGCGATTCCCAACCAGTCCCGCCCCCGACCCATGACGCGATGGTCGCCGGAAACAAAAAAATCCGGGAAGCGAAGATTGCCGGGTTTCCCCCCGCCat                                                                                                                                              <  1:515661/146‑1 (MQ=255)          acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc                                                                                                                                        >  1:493089/1‑151 (MQ=255)                 aaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGc                                                                                                                                 >  2:555924/1‑151 (MQ=255)                   ggggtggccagggCTGAAAATACAGCGAATTTTTTGTCAGCCATTCCCAACCATCCCCGCCCCCAACCCATGACGCGATGGTCGCGGGGAACAAAATAATACTGGAAGCGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGc                                                                                                                                 <  1:276832/138‑1 (MQ=255)                      gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                             >  2:768253/1‑150 (MQ=255)                        tGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                          >  1:991695/1‑151 (MQ=255)                              gccagaaaaTACGCCAAATTTGTTGTCAGCCAGTCCAAACAGGTCCCGCCCCCGACCCATGCCGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAAATGGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                                        <  1:866751/141‑1 (MQ=255)                                cTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                                  >  2:608528/1‑151 (MQ=255)                                   cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                               >  1:219793/1‑151 (MQ=255)                                    aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                              >  2:344918/1‑151 (MQ=255)                                    aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                              >  1:680256/1‑151 (MQ=255)                                    aaaTACAGCAAATTTTTTGTCAGCCAGTCCCACCCAGTCCCGCCCCCAAACCATGACGGGATGGTCGCCGGAAACAAAATAATACTGAAAGCGAAGATTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                 <  2:697448/148‑1 (MQ=255)                                    aaaTACAGCAAATTTTTTGCCACCCAGTCCCAACCAGTCCCCCCCCCAAACCATGACGCGTTGGTCGCGGGAAACAAAATAATCCTGGAGGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                               <  2:225495/150‑1 (MQ=255)                                          aGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                        >  1:520585/1‑151 (MQ=255)                                              aaaTTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                                    >  2:1006641/1‑151 (MQ=255)                                              aaaTTGTTGCCCCCCAGTCCCAACCAGTACCCCCCCCGAACCTTGACGCGATGGTCGCCGGAAACAAAATAATATTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                                     <  2:932540/150‑1 (MQ=255)                                                  tgttgtCACCCAGTCCCAACCAGTCCCCCCCCCAACCCATGACGGGATGGTCGCGGGAAAAAGAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAACGTGTGCTCTGTGCAGCATAGACacg                                                                                                 <  2:499389/150‑1 (MQ=255)                                                               ttCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                                   >  2:297696/1‑151 (MQ=255)                                                                cccacccatccccccccccAAACCATGACGGGTGGGCCGCCGGGAACAAAAAAATACTGAAAGAGAAGATGGCGGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac                                                                                            <  1:825917/127‑1 (MQ=255)                                                                  cAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                                >  2:763833/1‑151 (MQ=255)                                                                         tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                          >  1:881776/1‑150 (MQ=255)                                                                               ccccAACCCATGACGCGATGGTCGCGGGGAAAAAAATAATCCTGGAAGCGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAc                                                                   <  2:580368/151‑1 (MQ=255)                                                                                    aaCCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                               >  1:629177/1‑150 (MQ=255)                                                                                    aaCCATGACGAGAGGGTCGCGGGAAAAAAAAAAATCCGGAAGGCGAAGTTTGCGGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                               <  1:557633/150‑1 (MQ=255)                                                                                      ccATGGCGCGAGGGGCGCCGGAAACAAAATAATTCGGGACGCAAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                               <  1:278054/148‑1 (MQ=255)                                                                                            cgcgATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGt                                                      >  1:75224/1‑151 (MQ=255)                                                                                             gcgATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                     >  2:264592/1‑151 (MQ=255)                                                                                               gATGGTCGCGGGAAACAGAATAATACTGGAGGCGAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                     <  2:860589/149‑1 (MQ=255)                                                                                                aTGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                     >  2:34505/1‑84 (MQ=255)                                                                                                      gcggggAAAAAAATAATACGGGAAGGAAAGTTTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                              <  1:426262/146‑1 (MQ=255)                                                                                                          ggAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                        >  2:975988/1‑151 (MQ=255)                                                                                                          gaaaCAAAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                        <  2:915242/149‑1 (MQ=255)                                                                                                           aaaaaaaataatCCTGAAGGCAAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                                       <  2:887555/145‑1 (MQ=255)                                                                                                            aaaataataatACTGAAAGCGAAGTTTGCCGTTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaa                                         <  1:45953/143‑1 (MQ=255)                                                                                                             aCAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                                      >  1:186136/1‑150 (MQ=255)                                                                                                              cAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                    >  2:777712/1‑151 (MQ=255)                                                                                                               aaaataatACTGGAAGGGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                                   <  2:438730/149‑1 (MQ=255)                                                                                                               aaaataatACTGGAAGAGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                                   <  2:442318/149‑1 (MQ=255)                                                                                                               aGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                    >  1:791060/1‑150 (MQ=255)                                                                                                                     ataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCACATTCACTTTCCGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTGGACCACGCTCAACAAATacaaca                              >  1:268721/1‑150 (MQ=255)                                                                                                                      taCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             >  2:579253/1‑150 (MQ=255)                                                                                                                        cTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             >  1:949075/1‑148 (MQ=255)                                                                                                                          ggAAGCGAAGATGGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                          <  1:696049/149‑1 (MQ=255)                                                                                                                           gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                       >  1:434554/1‑151 (MQ=255)                                                                                                                             aGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                       <  2:196239/149‑1 (MQ=255)                                                                                                                                  aGATGGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa                         <  2:754957/142‑1 (MQ=255)                                                                                                                                   gATTGCCGTGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg                 <  1:988030/149‑1 (MQ=255)                                                                                                                                     ttGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa              <  2:571250/150‑1 (MQ=255)                                                                                                                                            gATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat      <  2:814537/151‑1 (MQ=255)                                                                                                                                             aTTACCCCCGCTATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat      <  1:45609/150‑1 (MQ=255)                                                                                                                                             aTTACACCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat      >  1:407107/1‑150 (MQ=255)                                                                                                                                                aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTg  >  2:735228/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                       CATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443134‑3443427

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