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

CAGACGCATAGAGTAACACATAAAGCGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACT  >  NC_000913/3443128‑3443426                                                                                                                                                        |                                                                                                                                                      cAGACGCATAGAGTAACACATAAAGCGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAc                                                                                                                                                      >  1:142516/1‑151 (MQ=255)       cATAGAGTAACACATAAAGCGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGc                                                                                                                                                >  1:467897/1‑151 (MQ=255)             gTAACACATAAAGCGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATAtt                                                                                                                                          >  2:881231/1‑151 (MQ=255)               aaCACATAAAGCGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTc                                                                                                                                         >  2:1148439/1‑150 (MQ=255)                         gCGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCgg                                                                                                                              >  1:156114/1‑151 (MQ=255)                          gggTTGCCCCGGGCTAAAAATACAGAAAATTTTTTGTCA‑CCCTTTCCAACCCAGCCCCCCCCCCAACCCATGGCGGGTGGGCCGCCGGAAACAAAATAATCCGGAAAGCAAAGAGTGCGGGTTTTCCCCCCGCCATATTCACTTTCAGCg                                                                                                                               <  1:347915/149‑1 (MQ=255)                          cGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGt                                                                                                                             >  2:214318/1‑151 (MQ=255)                              tG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                         >  2:555799/1‑151 (MQ=255)                              tG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                         >  2:7177/1‑151 (MQ=255)                                 cccAGGCTGAAAATACAGCAAATTTGTTGCCAGCCAGTCCCAA‑CCAGTCCCCCCCCCGAACCATGCCGGGATGGCCGCCGGAAAAAAAATAATACTGGAAGCGAAGATGGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                         <  2:430937/149‑1 (MQ=255)                                  ccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg                                                                                                                      >  1:270415/1‑151 (MQ=255)                                           aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                             >  1:57699/1‑151 (MQ=255)                                           aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                              >  1:210810/1‑150 (MQ=255)                                           aaaTACAGCGAAACTGTTGCCAGCCAGTTCCAA‑CCAGGACCGCACCCGAACCATGATGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                             >  1:214293/1‑151 (MQ=255)                                              tACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                            >  2:729697/1‑149 (MQ=255)                                                     aatttgttgTCAGCCAGTTCCAC‑CATTTCCCCCCCCCGAACCATGTCGCGGTGGTCCCCGGAAAAAGAATAATACTGGAGGCAAAGTTTGCCGGCATTTCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                                    <  2:1150035/147‑1 (MQ=255)                                                       aTTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                                 >  2:829841/1‑151 (MQ=255)                                                       aTTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                                 >  2:69889/1‑151 (MQ=255)                                                       aTTGTTGGCGGCCAGTCCCAA‑TCAGAACCCCCCCAGAACCATCATTAAGGGGCCGTGGGAAACAGAAAAAGTCGGGAAGCAAAGAGTGGGGTGCTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                                    <  2:574392/148‑1 (MQ=255)                                                           ttcccAGCCAGTCCCAC‑CCAGCCCCCCCCCCGAACCAGGACGGGGGGGCCGCGGGAAAAAAAAAAATCCTGGAAGCGAAGTTGGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                              <  1:859358/146‑1 (MQ=255)                                                                aGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                        >  1:122089/1‑151 (MQ=255)                                                                aGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                        >  2:1001539/1‑151 (MQ=255)                                                                  ccgttcccAC‑CCGGTCCCCCCCCCGACCCATGGGGGGAGGGTCGCGGGGAAAAAAAAAATACTGAAGGGGAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                        <  2:947682/143‑1 (MQ=255)                                                                  ccggtcccAC‑CCAGTCCCCCCCCCGAACAATGACGGGAGGGTCGCGGGGAACAAAAAAATGCTGGAGGCGAGGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                        <  2:436337/143‑1 (MQ=255)                                                                  ccAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                                      >  1:737059/1‑151 (MQ=255)                                                                     ttcccAA‑CAAGCCCCCCCCCCGACCGTTGCCGGGGTGGTCGCCGGAAAAATAATAATACGGGAGGAGAATTTGGGGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc                                                                                          <  1:227599/141‑1 (MQ=255)                                                                      ttCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                                   >  2:751710/1‑150 (MQ=255)                                                                       cccAA‑CAAGCCCCCCCCCCAACCCATGACGCGTGGGCCGCCGGGAACAAAATATTACTGGAAGCGAAGTTTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                                  <  2:208206/149‑1 (MQ=255)                                                                                   cccccgccGACCCAGGACGCGATGGTCGCCGGGAAAAGAATAATACTGGAAGGGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                                      <  2:564891/129‑1 (MQ=255)                                                                                    cccccccGACCCAGGACGCGAGGGTCGCGGGAAAAAAAAAAATCCTGGAAGCAAAGTTTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                                     <  2:129567/149‑1 (MQ=255)                                                                                       ccccGACCAATGACGGGTTGGTCGCCGGAAACAAAATAATACGGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAc                                                                  <  2:1134906/151‑1 (MQ=255)                                                                                           gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                               >  2:941206/1‑150 (MQ=255)                                                                                                   aCGCGAGGGCCGCCGGAAAAAAAATAATTGGGGAGGAGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                       <  1:52070/150‑1 (MQ=255)                                                                                                     gcgATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTACCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGCGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                                      >  1:830138/1‑149 (MQ=255)                                                                                                        aTGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                                 >  1:691271/1‑151 (MQ=255)                                                                                                          ggtcgcgagaaaaaaaataatATTGAAGGAAAAGATCGCCGGAATCCCCGCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc                                                <  2:189525/135‑1 (MQ=255)                                                                                                          ggccgcgGGGAAAAAAAAAATACGGGAAGCGAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCa                                               <  2:689379/148‑1 (MQ=255)                                                                                                           gTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                              >  2:588319/1‑151 (MQ=255)                                                                                                                 gggAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                                         >  2:505075/1‑150 (MQ=255)                                                                                                                    aaCAAAAAAATACTGGAAGAGAGGATGCCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                                     <  1:927616/151‑1 (MQ=255)                                                                                                                        gataaatcatgaaaGCGAAGATTGGGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                   <  2:969129/137‑1 (MQ=255)                                                                                                                            aataCTGGAAGCAAAGATGGCGGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaa                               <  1:125138/149‑1 (MQ=255)                                                                                                                             ataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                             >  1:492651/1‑150 (MQ=255)                                                                                                                             ataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                             >  2:910666/1‑150 (MQ=255)                                                                                                                                cTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                            >  1:345996/1‑148 (MQ=255)                                                                                                                                 tGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                         >  2:512056/1‑150 (MQ=255)                                                                                                                                  ggAAGAAAAGGTGGCCGGAATACCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                         <  1:178334/149‑1 (MQ=255)                                                                                                                                   gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGCAAATGTGTGCTCTGTGCAGCATAGACACGACTACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                         >  2:263530/1‑148 (MQ=255)                                                                                                                                    aagcaaatatgccggggTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                     <  1:588290/137‑1 (MQ=255)                                                                                                                                         aaGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCa                   <  2:467918/148‑1 (MQ=255)                                                                                                                                             ttGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa              <  2:868131/149‑1 (MQ=255)                                                                                                                                              tGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat           >  2:762005/1‑151 (MQ=255)                                                                                                                                               gcggggTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata          <  1:878423/148‑1 (MQ=255)                                                                                                                                               gCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat           <  2:373793/150‑1 (MQ=255)                                                                                                                                                  gggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact        >  1:66111/1‑150 (MQ=255)                                                                                                                                                  gggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact        >  2:76186/1‑150 (MQ=255)                                                                                                                                                       tACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatact  >  2:592647/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CAGACGCATAGAGTAACACATAAAGCGGTTG‑CCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACT  >  NC_000913/3443128‑3443426

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