Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	344,674	T→C	Y167H (TAC→CAC)	yahL →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	344,674	0	T	C	96.5%	195.4 / ‑0.5	57	coding (499/816 nt)	yahL	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (26/29);  minor base . (1/1);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GTATATGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTA  >  NC_000913/344545‑344802                                                                                                                                  |                                                                                                                                 gTATATGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGc                                                                                                                         >  1:153403/1‑139 (MQ=255)    tatGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCtt                                                                                                                      >  1:1209114/1‑139 (MQ=255)      tGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTg                                                                                                                    >  2:543124/1‑139 (MQ=255)      tGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGA‑‑ACCATATGCCTTTGCt                                                                                                                  >  2:682039/1‑139 (MQ=255)            ggATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTGCCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTgaga                                                                                                              >  1:1644089/1‑139 (MQ=255)                  tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAGCACCATATGCCTTTGCTGAGAATTTTc                                                                                                        >  2:954001/1‑139 (MQ=255)                  tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTc                                                                                                        >  1:1366763/1‑139 (MQ=255)                  tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTc                                                                                                        >  2:1444781/1‑139 (MQ=255)                    cGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTccc                                                                                                      <  2:581496/139‑1 (MQ=255)                      tttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCa                                                                                                                              <  1:381048/113‑1 (MQ=255)                      tttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCa                                                                                                                              >  2:381048/1‑113 (MQ=255)                       ttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTaa                                                                                                   >  2:1331468/1‑139 (MQ=255)                       ttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACACATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTaa                                                                                                   >  1:883481/1‑139 (MQ=255)                            gATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                              <  2:1098021/139‑1 (MQ=255)                            gATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                              <  2:555254/139‑1 (MQ=255)                            gATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGA‑‑ACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            <  1:682039/139‑1 (MQ=255)                             atatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaa                                                                                             <  2:883481/139‑1 (MQ=255)                              tatTAAAGCAAACGTGCAGGCTGGACGCTTTGCCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            <  2:1644089/139‑1 (MQ=255)                              tatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            <  1:543124/139‑1 (MQ=255)                              tatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            <  1:1263341/139‑1 (MQ=255)                              tatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            <  2:1001149/139‑1 (MQ=255)                                      cAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTAt                                                                                    <  1:1331468/139‑1 (MQ=255)                                        aaCGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            <  1:203391/129‑1 (MQ=255)                                        aaCGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            >  2:203391/1‑129 (MQ=255)                                        aaCGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGc                                                                                  >  2:270982/1‑139 (MQ=255)                                          cGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAt                                                                                <  1:1444781/139‑1 (MQ=255)                                             gCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                    >  1:1226857/1‑116 (MQ=255)                                             gCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                    <  2:1226857/116‑1 (MQ=255)                                                      aCGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATgg                                                                    >  2:300153/1‑139 (MQ=255)                                                           ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            >  1:1643159/1‑110 (MQ=255)                                                           ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                            <  2:1643159/110‑1 (MQ=255)                                                               cGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaa                                                           >  1:783332/1‑139 (MQ=255)                                                                              aGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                    <  1:1265607/83‑1 (MQ=255)                                                                              aGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                    >  2:1265607/1‑83 (MQ=255)                                                                                  cAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgaga                                        >  2:236608/1‑139 (MQ=255)                                                                                      ttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                    <  1:300153/139‑1 (MQ=255)                                                                                        aTCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcg                                  <  2:1366763/139‑1 (MQ=255)                                                                                        aTCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcg                                  <  1:288730/139‑1 (MQ=255)                                                                                         tCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                 <  2:153403/139‑1 (MQ=255)                                                                                          cGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                <  2:141529/139‑1 (MQ=255)                                                                                              tCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                            >  2:492185/1‑139 (MQ=255)                                                                                               cAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                           >  1:1337149/1‑139 (MQ=255)                                                                                               cAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                           <  1:492185/139‑1 (MQ=255)                                                                                                  tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTccc                                                                                                      <  2:681996/61‑1 (MQ=255)                                                                                                  tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTccc                                                                                                      >  1:681996/1‑61 (MQ=255)                                                                                                  tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTaa                        >  1:654695/1‑139 (MQ=255)                                                                                                   gCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTcc                                             >  1:50448/1‑117 (MQ=255)                                                                                                   gCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTcc                                             <  2:50448/117‑1 (MQ=255)                                                                                                           aTTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCt               >  2:1370891/1‑139 (MQ=255)                                                                                                                aaaTGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTgcgc          <  2:1209114/139‑1 (MQ=255)                                                                                                                    gTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATgg      <  2:1337149/139‑1 (MQ=255)                                                                                                                     tATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGa     >  1:586719/1‑139 (MQ=255)                                                                                                                      aTTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGat    <  2:654695/139‑1 (MQ=255)                                                                                                                      aTTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGat    <  2:783332/139‑1 (MQ=255)                                                                                                                       ttcttcTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa         >  1:1573873/1‑133 (MQ=255)                                                                                                                       ttcttcTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa         <  2:1573873/133‑1 (MQ=255)                                                                                                                        tcttctTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatta  <  1:270982/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                 GTATATGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTA  >  NC_000913/344545‑344802

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