Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,068,728	C→A	6.4%	M27I (ATG→ATT)	argO ←	arginine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,068,728	0	C	A	6.4%	212.6 / 2.2	62	M27I (ATG→ATT)	argO	arginine transporter
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base A (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATT  >  NC_000913/3068595‑3068864                                                                                                                                      |                                                                                                                                         ccAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATcaca                                                                                                                                    <  1:1625455/140‑1 (MQ=255)  cAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGAGTCATCACaa                                                                                                                                   >  1:1649584/1‑140 (MQ=255)      gcCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGc                                                                                                                               <  1:917909/140‑1 (MQ=255)      gcCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGc                                                                                                                               <  2:1507688/140‑1 (MQ=255)      gcCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGc                                                                                                                               <  2:1507669/140‑1 (MQ=255)      gcCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGc                                                                                                                               <  1:918038/140‑1 (MQ=255)       ccAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCa                                                                                                                              >  1:580210/1‑140 (MQ=255)                      aCTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTTACGACGTATGCCCTGATTAATCACAAAAGCATTTTGTGGACCGAGc                                                                                                               <  2:551750/140‑1 (MQ=255)                      aCTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTTACGACGTATGCCCTGATTAATCACAAAAGCATTTTGTGGACCGAGc                                                                                                               <  2:551468/140‑1 (MQ=255)                         gCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTCAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGt                                                                                                            <  1:305185/140‑1 (MQ=255)                         gCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGt                                                                                                            <  2:1503904/140‑1 (MQ=255)                         gCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGt                                                                                                            <  2:1065839/140‑1 (MQ=255)                         gCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGt                                                                                                            <  2:102556/140‑1 (MQ=255)                         gCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGt                                                                                                            <  1:1281370/140‑1 (MQ=255)                         gCATCAATAACGCGCTGCCAACAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGt                                                                                                            <  1:1279352/140‑1 (MQ=255)                                     cgcTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAgcc                                                                                                >  2:1667826/1‑140 (MQ=255)                                        tGCCACCAAAAATCCCGGCGCAAATCAGGACCAACTCGCTGATAGCACAAAGTAAGGCCATCATAAGGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTGTGTGGACCGAGCGGTCGGATCAGAgccgcc                                                                                             >  2:1392122/1‑140 (MQ=255)                                        tGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAgccgcc                                                                                             >  2:1392346/1‑140 (MQ=255)                                          ccaccaAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGcccc                                                                                           <  1:822835/140‑1 (MQ=255)                                          ccaccaAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGcccc                                                                                           <  2:274691/140‑1 (MQ=255)                                          ccaccaAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGcccc                                                                                           <  1:1043227/140‑1 (MQ=255)                                          ccaccaAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGcccc                                                                                           <  1:1102193/140‑1 (MQ=255)                                                                tCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCAGCCCCAAGTGCAAGACCTTGAAAGtaa                                                                     <  2:963313/140‑1 (MQ=255)                                                                    gACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAg                                                                 >  1:1027189/1‑140 (MQ=255)                                                                         aaTCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCTCAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAac                                                            >  1:8543/1‑140 (MQ=255)                                                                            cGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACg                                                         <  1:162650/140‑1 (MQ=255)                                                                            cGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACg                                                         <  1:162785/140‑1 (MQ=255)                                                                              cTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGtt                                                       <  2:1483561/140‑1 (MQ=255)                                                                              cTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGtt                                                       >  2:169576/1‑140 (MQ=255)                                                                              cTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGtt                                                       <  2:176260/140‑1 (MQ=255)                                                                              cTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGtt                                                       <  2:176344/140‑1 (MQ=255)                                                                                         aaGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGa                                            >  2:203293/1‑140 (MQ=255)                                                                                            tAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAAtt                                         <  2:91677/140‑1 (MQ=255)                                                                                              aGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTaa                                       <  1:1667826/140‑1 (MQ=255)                                                                                                    tCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCt                                 >  1:89797/1‑140 (MQ=255)                                                                                                    tCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCt                                 >  1:89548/1‑140 (MQ=255)                                                                                                    tCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCt                                 >  1:89353/1‑140 (MQ=255)                                                                                                      aTAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGa                               >  2:747346/1‑140 (MQ=255)                                                                                                      aTAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGAATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGa                               >  2:982102/1‑140 (MQ=255)                                                                                                         aTGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGc                            <  2:1260405/140‑1 (MQ=255)                                                                                                         aTGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGc                            <  2:817363/140‑1 (MQ=255)                                                                                                            tgGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCaga                         <  1:1098463/140‑1 (MQ=255)                                                                                                             ggTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCagag                        <  1:502982/140‑1 (MQ=255)                                                                                                             ggTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCagag                        <  1:1085948/140‑1 (MQ=255)                                                                                                              gaaCTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGTATTAAGTATCTGATGCagagt                       >  2:18729/3‑139 (MQ=255)                                                                                                              gTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGACAGTAATAAGAAAACTCGTTGGATATTCCGAATGAAGTATCTGATGCagata                       >  2:1057342/1‑138 (MQ=255)                                                                                                              gTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACAGGTTGGATATTCCGACTTCAGTATCTGATGCAgaga                       >  2:1057208/1‑140 (MQ=255)                                                                                                              gTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAgaga                       >  1:1537792/1‑140 (MQ=255)                                                                                                              gTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAgaga                       >  2:1082646/1‑140 (MQ=255)                                                                                                               tACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGAc                      >  1:1381265/1‑140 (MQ=255)                                                                                                                 cTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACta                    >  1:329561/1‑140 (MQ=255)                                                                                                                    acgacgTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTatat                 >  1:538164/1‑140 (MQ=255)                                                                                                                    acgacgTATGCCCTGATTCAACACAAAAGCATATAGTGGACCGAGCGGTAGGATCATAGACGCCCCAAGTGCAAGACCATGAAAGAAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTatat                 >  1:177463/1‑140 (MQ=255)                                                                                                                      gccgTATGCCCTGATTCATCACAAAAGCATTGTGTGGACCGAGCGGTAGGATCAGAGCCGCCCCAAGTGCAAGACCTTGAAAGTACTAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATcg               <  1:979453/138‑1 (MQ=255)                                                                                                                      gacgTATGCCCTGATTCATCACAAAAGCCTTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTACTAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATcg               <  1:979454/140‑1 (MQ=255)                                                                                                                        cgTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGCAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATcgcg             <  1:1664208/140‑1 (MQ=255)                                                                                                                         gTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCgg            >  1:1120582/1‑140 (MQ=255)                                                                                                                           aTGCCCTGATTAATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAATCACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGa          >  1:1124624/1‑140 (MQ=255)                                                                                                                           aTGCCCTGATTAATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGa          >  1:701046/1‑140 (MQ=255)                                                                                                                              cccTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTg       >  2:1145716/1‑140 (MQ=255)                                                                                                                                   aTTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGttatt  >  1:1329580/1‑140 (MQ=255)                                                                                                                                   aTTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGttatt  >  1:228770/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                         CCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATT  >  NC_000913/3068595‑3068864

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