Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,424	0	C	A	48.0%	65.4 / 31.1	51	V37F (GTT→TTT)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base C (22/4);  new base A (0/24);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.72e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAA  >  NC_000913/2248276‑2248557                                                                                                                                                     |                                                                                                                                       ccGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGcaaca                                                                                                                                       >  1:786133/1‑150 (MQ=255)  cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGcaacaa                                                                                                                                      >  2:341135/1‑150 (MQ=255)            aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                            >  1:921367/1‑150 (MQ=255)                gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑g                                                                                                                       >  1:1034616/1‑150 (MQ=255)                    cGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATgg                                                                                                                   >  1:226976/1‑150 (MQ=255)                      gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAAc                                                                                                                <  2:716254/151‑1 (MQ=255)                      gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAAc                                                                                                                <  2:580491/151‑1 (MQ=255)                       aaaCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAAcc                                                                                                               <  1:410741/151‑1 (MQ=255)                         aCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGc                                                                                                             <  2:660900/151‑1 (MQ=255)                            ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCa                                                                                                           >  1:935662/1‑150 (MQ=255)                            ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCa                                                                                                           >  1:629919/1‑150 (MQ=255)                                         ccAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGt                                                                                             >  1:1038296/1‑151 (MQ=255)                                           atttCCCCGAGCCTGGTCATCAGAAAGAAACCCACCGGCCCAACAAGAATATACGAGAGAAATCCCCCCCCGGCCCCTCCCTGAAAAATCTTTCCGCCAGAGCAAAAAAAAAGCCCTGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCa                                                                                           <  1:891001/149‑1 (MQ=255)                                             ttCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAgg                                                                                         >  2:871081/1‑151 (MQ=255)                                               cACCGAGACTGGTCATAAGAAAGTAACCCATCAGGCCAATCAGCATATACGAGAGAAATCCCCCCCCCGGACCTGCCTGAAAAATCGTGGCGCCAGAGGAAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTg                                                                                       <  1:826969/151‑1 (MQ=255)                                                     gaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACgcgc                                                                                 >  1:144046/1‑151 (MQ=255)                                                         gttcatcaGGAAAAAACCCACCAGGCCAACCAGCATATACGAGAAAAATCCCCCGCCCGGACCTGCCTGAAAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                   <  2:951759/143‑1 (MQ=255)                                                         gttcataggaaaGAAAACCATCAGGCAAACAAAAATAGCCGAGAGCAATCCCCCGCCGGCCCCTCCCTGAAAAATCGTGGCGCCAGAGACAAAAAAAAGCCCTGCCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                             <  2:514941/141‑1 (MQ=255)                                                            catcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAt                                                                          >  1:668515/1‑151 (MQ=255)                                                             atcaGGAAGAAACCCACCGGGCAAAACAACATAAAGGAGACAAACCCCCCGCCCGGCCCTCCGTGAAAAACCGTGGCGCAAGAGAAAAAAAAAAGCCCTTTCCG‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTaa                                                                           <  2:40861/149‑1 (MQ=255)                                                                agaaaGAAAACCACCAGCCCAATCAGAATATACGAGAGCATTCCCCCGCCCGGCCCTGCCTGAAAAATCGTTGCGCAAGAGGCAAAAAAAAGCCCTGTACC‑TATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                         <  2:391990/145‑1 (MQ=255)                                                                   aaataaccctacggaccAATCAAAATATACAAGACAAATGCCCCACCCGGCCCTGCCTAAAAATTTGTTGCGCAAGATGCAA‑AAAAAGGCCCGTTCCGGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGc                                                                   <  2:284285/136‑1 (MQ=255)                                                                          ccATCGGGCAAAACACCAAATAGGAGAGCAACGCCCCCCCCGGCCCTGCTGGAGAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGcc                                                            <  2:935652/151‑1 (MQ=255)                                                                                     accAGAATATCCGAGACAAACGCCCCCCCCGGCCCGGCATAAAAAATCGTTGCCCCAGAGACAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                  <  2:835239/148‑1 (MQ=255)                                                                                     aTCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtg                                                 >  2:558990/1‑151 (MQ=255)                                                                                     aTCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtg                                                 >  2:221485/1‑151 (MQ=255)                                                                                            tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGCGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                          >  2:1002868/1‑151 (MQ=255)                                                                                               acaagacCAATCCCCCCCCCGGCCCTGCCTGAAAAACCGTGGCGCCAGGGCAAAAAAAAAGCCCTGTCCG‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                          <  2:871017/145‑1 (MQ=255)                                                                                               aCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                       >  2:340201/1‑151 (MQ=255)                                                                                                    aactgccccccGCCGGGCCCAGCCGGAGAAGGGGTGGCGCCGGAGCAAAAAAAAAGACCTGTCCT‑ATTGGACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                       <  2:904605/139‑1 (MQ=255)                                                                                                       aacgCCCCGCCCGGCCCTGCCGGAGAATTCGTTGCGCAAGGGCCAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                <  1:528128/147‑1 (MQ=255)                                                                                                       aaccccccGCCCGGCCCTGCCGGAAAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                <  1:454186/146‑1 (MQ=255)                                                                                                       aaTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                >  1:40769/1‑150 (MQ=255)                                                                                                       aaTCCCCCGCCCGGACCTGCTTGAAAAATCGTGGCGCCAGAGGAAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGa                                 <  1:626982/149‑1 (MQ=255)                                                                                                           ccccccccGGCCCTGCCAGAAAAATCGTTGCGCCAGAGACAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCa                            <  1:180473/150‑1 (MQ=255)                                                                                                           ccccGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                           >  1:772918/1‑151 (MQ=255)                                                                                                             ccgccggAACCTGCATGAAAAATCGTGGCGCAAGAGGAAAAAAAAGGCCTTGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCa                            <  2:806273/148‑1 (MQ=255)                                                                                                                   ggACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAAt                   >  1:369532/1‑151 (MQ=255)                                                                                                                    gACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                  >  2:145991/1‑151 (MQ=255)                                                                                                                    gACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                  >  1:843955/1‑151 (MQ=255)                                                                                                                      cctgcctgAAAAACCTTTGCGCCAGAGAAAAAAAAAACACCTGTACC‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                  <  1:86872/149‑1 (MQ=255)                                                                                                                      ccgccctaaaaaaTCTTGGCGCCAGAGCAAAAAAAAAGCCCTGCACC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                  <  2:9351/139‑1 (MQ=255)                                                                                                                         gcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc             >  1:882286/1‑151 (MQ=255)                                                                                                                          cctgAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc             <  1:643909/150‑1 (MQ=255)                                                                                                                              aaaaaTCGTTGCGCAAGAGGAAAAAAAAAGCCCGTTACG‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc            <  2:993172/145‑1 (MQ=255)                                                                                                                               gaAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt       >  1:663262/1‑151 (MQ=255)                                                                                                                                aactcgtcccgcgAGAGAAAAAAAAAAGCCCTGTCCG‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa      <  1:870262/142‑1 (MQ=255)                                                                                                                                aaacggTTGCGCCAGAGCAAAAAAAAAGCCCGGTACG‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt       <  2:959810/145‑1 (MQ=255)                                                                                                                                aaaTCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt       >  1:548667/1‑150 (MQ=255)                                                                                                                                aaaTCGTGGCCCCAGAGGCAAAAAAAAGCCCTGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt       <  2:946506/150‑1 (MQ=255)                                                                                                                                 aaTCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGAGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa      >  2:676818/1‑150 (MQ=255)                                                                                                                                     tgtGCGCGAGGGCCAAAAAAATGACCTGCACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaa     <  2:194536/145‑1 (MQ=255)                                                                                                                                         cgcAAGAGGAAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaaa  <  1:245874/146‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                       CCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAA  >  NC_000913/2248276‑2248557

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