Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,424,111	T→A	S110S (TCT→TCA)	trkG →	Rac prophage; potassium transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,424,111	0	T	A	96.8%	213.9 / ‑1.6	63	S110S (TCT→TCA)	trkG	Rac prophage; potassium transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (31/30);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCTGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAATTTTATAGGTGGT  >  NC_000913/1424015‑1424207                                                                                                 |                                                                                                 gTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGa                                                                                               >  2:958429/1‑100 (MQ=255)  ttttGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGaa                                                                                              <  1:1230461/100‑1 (MQ=255)   tttGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAAt                                                                                             <  2:235214/100‑1 (MQ=255)   tttGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAAt                                                                                             >  1:1289578/1‑100 (MQ=255)    ttGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATa                                                                                            >  1:862399/1‑100 (MQ=255)    ttGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATa                                                                                            <  1:1340275/100‑1 (MQ=255)        aTTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATaacaa                                                                                        >  1:28506/1‑100 (MQ=255)          tttGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATaacaaca                                                                                      <  1:1059401/100‑1 (MQ=255)           ttGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATaacaacaa                                                                                     >  1:1446676/1‑100 (MQ=255)               tttCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAgg                                                                                 <  1:736564/100‑1 (MQ=255)                ttCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGa                                                                                <  2:1260127/100‑1 (MQ=255)                 tCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAg                                                                               <  1:1254958/100‑1 (MQ=255)                 tCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAg                                                                               >  2:1510156/1‑100 (MQ=255)                 tCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAg                                                                               <  1:424469/100‑1 (MQ=255)                  cTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGc                                                                              <  1:1117810/100‑1 (MQ=255)                    gTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCaa                                                                            >  2:857898/1‑100 (MQ=255)                       attaGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTg                                                                         >  2:473114/1‑100 (MQ=255)                         taGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTa                                                                       <  1:757928/100‑1 (MQ=255)                            tGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAAtt                                                                    <  1:1398025/100‑1 (MQ=255)                             gCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATtg                                                                   <  1:135860/100‑1 (MQ=255)                             gCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTCTCTGGAATAACAACAACAGGAGCAACTGTAATtg                                                                   >  2:1007479/1‑100 (MQ=255)                               aTTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATtgat                                                                 <  1:1173982/100‑1 (MQ=255)                                ttCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATtgatg                                                                >  1:352310/1‑100 (MQ=255)                                ttCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATtgatg                                                                <  1:857903/100‑1 (MQ=255)                                ttCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATtgatg                                                                <  1:333526/100‑1 (MQ=255)                                   ccTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGTGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGatgatg                                                             >  2:1407638/1‑100 (MQ=255)                                         tGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGtt                                                       <  1:1311560/100‑1 (MQ=255)                                           gATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCa                                                     >  1:700161/1‑100 (MQ=255)                                              tGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTa                                                  <  1:384586/100‑1 (MQ=255)                                                aCTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTCTCTGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTAcc                                                <  1:1007485/100‑1 (MQ=255)                                                 cTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCt                                               >  2:161457/1‑100 (MQ=255)                                                  tCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTc                                              <  2:719776/100‑1 (MQ=255)                                                  tCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTc                                              >  1:936257/1‑100 (MQ=255)                                                  tCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTc                                              <  2:936252/100‑1 (MQ=255)                                                   cAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCg                                             >  2:404248/1‑100 (MQ=255)                                                   cAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCg                                             >  2:489609/1‑100 (MQ=255)                                                   cAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCg                                             >  2:944583/1‑100 (MQ=255)                                                   cAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCg                                             >  2:1024627/1‑100 (MQ=255)                                                      aaCTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGc                                          <  2:1300076/100‑1 (MQ=255)                                                       aCTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCa                                         >  1:725231/1‑100 (MQ=255)                                                         ttAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCata                                       >  1:1039311/1‑100 (MQ=255)                                                             tttAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTg                                   <  1:958434/100‑1 (MQ=255)                                                                  cGTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTa                              >  2:1118358/1‑100 (MQ=255)                                                                   gTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTAt                             >  1:594012/1‑100 (MQ=255)                                                                   gTTCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTAt                             >  1:850923/1‑100 (MQ=255)                                                                    ttCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATc                            <  1:508181/100‑1 (MQ=255)                                                                    ttCATTGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATc                            >  2:527473/1‑100 (MQ=255)                                                                        ttGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTc                        <  2:416323/100‑1 (MQ=255)                                                                        ttGATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTc                        >  1:822111/1‑100 (MQ=255)                                                                          gATGCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTcac                      <  1:843109/100‑1 (MQ=255)                                                                          gATGCTATGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGAAAATGTTAATTCATTTCCTCAGGCATATCTGTAATAACGGTcan                      <  2:1023488/100‑2 (MQ=255)                                                                            tGCTATGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAg                    <  2:317250/100‑1 (MQ=255)                                                                             gCTCTGTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGt                   >  1:948338/1‑100 (MQ=255)                                                                               tctGTTTAAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTa                 >  1:828522/1‑100 (MQ=255)                                                                                  gTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAAt              <  1:1315519/100‑1 (MQ=255)                                                                                  gTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAAt              <  1:1341422/100‑1 (MQ=255)                                                                                  gTTTGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAAt              >  1:686315/1‑100 (MQ=255)                                                                                    ttGAAGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAAttt            <  2:1371068/100‑1 (MQ=255)                                                                                       aaGGGGTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAATTTtat         <  1:1510165/100‑1 (MQ=255)                                                                                       aaGGGGTTTCAGGAATAACAACAACAGGAGCAACGGTAATTGATGAACTTAGTTCATTACCCCGGGCATATTTGTACTATcctcaacac                    >  2:1424993/1‑81 (MQ=255)                                                                                           ggTTTCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAATTTTATAggt     >  2:1193857/1‑100 (MQ=255)                                                                                             tttCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAATTTTATAggtgg   >  1:8408/1‑100 (MQ=255)                                                                                              ttCAGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAATTTTATAggtggt  >  2:1491175/1‑100 (MQ=255)                                                                                                 |                                                                                                 GTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCTGGAATAACAACAACAGGAGCAACTGTAATTGATGATGTTAGTTCATTACCTCGGGCATATTTGTACTATCGGTCACAGTTAAATTTTATAGGTGGT  >  NC_000913/1424015‑1424207

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑