Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,910,295	A→T	11.1%	D7V (GAC→GTC)	yebQ →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,910,295	0	A	T	11.1%	167.8 / 9.4	54	D7V (GAC→GTC)	yebQ	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base T (2/4);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGG  >  NC_000913/1910156‑1910441                                                                                                                                            |                                                                                                                                                   cGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaa                                                                                                                                           <  1:46590/149‑3 (MQ=255)  gCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaac                                                                                                                                          <  1:454135/149‑4 (MQ=255)    aCCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCAtt                                                                                                                                        <  3:403930/149‑1 (MQ=255)          gtggtCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCa                                                                                                                                  <  1:420771/149‑1 (MQ=255)             gtCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCg                                                                                                                               >  1:81361/1‑149 (MQ=255)                gCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGATATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATa                                                                                                                            >  4:292079/1‑149 (MQ=255)                gCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATa                                                                                                                            >  1:248988/1‑149 (MQ=255)                  aaTCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACg                                                                                                                          <  3:28523/149‑1 (MQ=255)                     ccATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTg                                                                                                                       <  2:466957/149‑1 (MQ=255)                      cATCTTTTGCCGGTTAGTTACAATTCTGCTACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCTAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGc                                                                                                                      >  3:501169/1‑149 (MQ=255)                               ccGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAAc                                                                                                             <  4:477269/149‑1 (MQ=255)                                   ttagttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCAtt                                                                                                         <  1:237102/149‑1 (MQ=255)                                   ttagttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCAtt                                                                                                         >  3:393548/1‑149 (MQ=255)                                       ttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGa                                                                                                     >  1:375833/1‑149 (MQ=255)                                         aCAATTCTGCGATATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAtt                                                                                                   <  3:292079/149‑1 (MQ=255)                                         aCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAtt                                                                                                   <  4:244702/149‑1 (MQ=255)                                          cAATTCTGCGGCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                                                                  <  4:501169/149‑1 (MQ=255)                                          cAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                                                                  >  4:287522/1‑149 (MQ=255)                                            aTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGt                                                                                                <  2:324317/149‑1 (MQ=255)                                            aTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGt                                                                                                >  3:83776/1‑149 (MQ=255)                                                    aCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGAt                                                                                        >  1:451929/1‑149 (MQ=255)                                                       tCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGc                                                                                     <  4:18400/149‑1 (MQ=255)                                                           ccGTGAATTTCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaac                                                                                                                                          <  4:417277/92‑4 (MQ=255)                                                           ccGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTc                                                                                 >  4:81392/1‑149 (MQ=255)                                                           ccGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGAAATAAAATGCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaac                                                                                                                                          >  3:417277/1‑89 (MQ=255)                                                             gTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTTGATGGCCGTCCt                                                                               >  2:292428/1‑149 (MQ=255)                                                                aaTATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaa                                                                                                                                           >  3:250121/1‑84 (MQ=255)                                                                aaTATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaa                                                                                                                                           <  4:250121/86‑3 (MQ=255)                                                                        tGCTAGAATCATACCCCTGTTGCTTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTCACGGCGCaa                                                                    <  1:512494/149‑1 (MQ=255)                                                                              aaTCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCa                                                              >  3:397360/1‑149 (MQ=255)                                                                                           ttGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGTCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGcc                                                 >  2:230982/1‑149 (MQ=255)                                                                                                tattCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCcaacaa                                            <  4:393548/149‑1 (MQ=255)                                                                                                 attCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAAt                                           <  2:58874/149‑1 (MQ=255)                                                                                                  ttCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATc                                          <  3:287522/149‑1 (MQ=255)                                                                                                    cACCAGAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGTTTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGc                                        <  4:300059/149‑1 (MQ=255)                                                                                                     aCCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGcc                                       >  3:115743/1‑149 (MQ=255)                                                                                                     aCCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGcc                                       >  2:102119/1‑149 (MQ=255)                                                                                                       cAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCAc                                     >  2:121306/1‑149 (MQ=255)                                                                                                           gATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGAc                                 <  4:83776/149‑1 (MQ=255)                                                                                                             tataAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCt                               >  3:188912/1‑149 (MQ=255)                                                                                                               taAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTc                             >  2:190586/1‑149 (MQ=255)                                                                                                               taAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTc                             >  2:555704/1‑149 (MQ=255)                                                                                                                 aaaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCAt                           <  4:192627/149‑1 (MQ=255)                                                                                                                 aaaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCAt                           <  1:555704/149‑1 (MQ=255)                                                                                                                        tATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgc                    >  3:380724/1‑149 (MQ=255)                                                                                                                        tATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgc                    >  4:415599/1‑149 (MQ=255)                                                                                                                            ccAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagc                <  2:81361/149‑1 (MQ=255)                                                                                                                              aaaaGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagcca              >  4:461471/1‑149 (MQ=255)                                                                                                                              aaaaGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagcca              >  4:498657/1‑149 (MQ=255)                                                                                                                               aaaGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGccagccag             >  4:171621/1‑149 (MQ=255)                                                                                                                                  gTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTc          >  2:371173/1‑149 (MQ=255)                                                                                                                                         ccGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGTCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTgg   <  1:230982/149‑1 (MQ=255)                                                                                                                                         ccGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTgg   <  3:498657/149‑1 (MQ=255)                                                                                                                                          cGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTggg  <  4:397360/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                   CGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGG  >  NC_000913/1910156‑1910441

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