Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,019,339	G→T	C60F (TGC→TTC)	ssnA →	putative aminohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,019,339	0	G	T	100.0%	205.9 / NA	59	C60F (TGC→TTC)	ssnA	putative aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/35);  total (24/35)

GGCAAAAGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCT  >  NC_000913/3019202‑3019468                                                                                                                                          |                                                                                                                                  ggCAAAAGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGtcttct                                                                                                                                 >  2:371070/1‑140 (MQ=255)    aaaaGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGc                                                                                                                              <  2:929618/140‑1 (MQ=255)     aaaGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTTATGCCAGGGATTGTCTTCTCGca                                                                                                                             <  1:885021/140‑1 (MQ=255)       aGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGcaca                                                                                                                           <  2:374586/140‑1 (MQ=255)       aGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGcaca                                                                                                                           <  1:280588/140‑1 (MQ=255)       aGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGcaca                                                                                                                           <  1:908518/140‑1 (MQ=255)              gAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCAttt                                                                                                                    <  1:620019/140‑1 (MQ=255)                aGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCAttttt                                                                                                                  >  1:394997/1‑140 (MQ=255)                   cGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACt                                                                                                               <  2:135201/140‑1 (MQ=255)                                  cGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCcgcg                                                                                                >  1:835942/1‑140 (MQ=255)                                  cGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCcgcg                                                                                                >  1:28438/1‑140 (MQ=255)                                     aaaCGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGaa                                                                                             <  2:327603/140‑1 (MQ=255)                                     aaaCGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGaa                                                                                             <  1:614820/140‑1 (MQ=255)                                           tgtgATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATgg                                                                                       >  2:110802/1‑140 (MQ=255)                                           tgtgATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATgg                                                                                       >  2:27648/1‑140 (MQ=255)                                             tgATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCa                                                                                     <  1:629797/140‑1 (MQ=255)                                                  gTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACAt                                                                                <  1:896620/140‑1 (MQ=255)                                                     gCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGc                                                                             >  2:852798/1‑140 (MQ=255)                                                        aTCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACAt                                                                                >  2:871369/1‑134 (MQ=255)                                                        aTCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACAt                                                                                <  1:871369/134‑1 (MQ=255)                                                        aTCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGcccc                                                                          <  2:1013223/140‑1 (MQ=255)                                                         tCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCt                                                                         >  2:963068/1‑140 (MQ=255)                                                             cGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGcc                                                                     <  1:852798/140‑1 (MQ=255)                                                               aTGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCg                                                                   <  2:835942/140‑1 (MQ=255)                                                               aTGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCg                                                                   <  2:631966/140‑1 (MQ=255)                                                                 gCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGa                                                                 >  1:466000/1‑140 (MQ=255)                                                                  cccTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGAt                                                                <  1:110802/140‑1 (MQ=255)                                                                  cccTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGAt                                                                <  1:446491/140‑1 (MQ=255)                                                                         gCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCAtct                                                         <  2:824072/140‑1 (MQ=255)                                                                            aCGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCaa                                                      >  2:456367/1‑140 (MQ=255)                                                                             cGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGcc                                                                     <  1:348982/124‑1 (MQ=255)                                                                             cGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGcc                                                                     >  2:348982/1‑124 (MQ=255)                                                                                  ccccGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGa                                                <  1:996202/140‑1 (MQ=255)                                                                                  ccccGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGa                                                >  2:567307/1‑140 (MQ=255)                                                                                  ccccGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGa                                                >  1:880635/1‑140 (MQ=255)                                                                                              cTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCtggt                                    >  1:902323/1‑140 (MQ=255)                                                                                                        aTGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGa                          >  1:38745/1‑140 (MQ=255)                                                                                                        aTGCATGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGa                          >  2:375542/1‑140 (MQ=255)                                                                                                             tGGCCGGATTGTGATGCCAGGGATTGTCTTCTCTCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATcgcg                     >  2:575315/1‑140 (MQ=255)                                                                                                             tGGCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATcgcg                     <  1:29755/140‑1 (MQ=255)                                                                                                              ggCCGGATTGTGATGCCAGGGATTGTCTTCTCTCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCgcgc                    <  1:575315/140‑1 (MQ=255)                                                                                                              ggCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCgcgc                    <  2:859997/140‑1 (MQ=255)                                                                                                              ggCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCgcgc                    <  1:494035/140‑1 (MQ=255)                                                                                                              ggCCGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCgcgc                    <  2:28438/140‑1 (MQ=255)                                                                                                                ccGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGccc                  >  1:705478/1‑140 (MQ=255)                                                                                                                ccGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGccc                  >  2:363725/1‑140 (MQ=255)                                                                                                                 cGGATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCt                 >  1:242333/1‑140 (MQ=255)                                                                                                                  ggATTGTGATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCtt                >  2:874095/1‑140 (MQ=255)                                                                                                                        tgATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGaa          >  2:1029851/1‑140 (MQ=255)                                                                                                                         gATGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAag         <  2:902323/140‑1 (MQ=255)                                                                                                                          aTGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCgcgc                    <  2:366588/128‑1 (MQ=255)                                                                                                                          aTGCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCgcgc                    >  1:366588/1‑128 (MQ=255)                                                                                                                            gCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGt      <  1:956810/140‑1 (MQ=255)                                                                                                                            gCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGt      <  2:466000/140‑1 (MQ=255)                                                                                                                            gCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGt      <  2:242333/140‑1 (MQ=255)                                                                                                                            gCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGt      <  1:1029851/140‑1 (MQ=255)                                                                                                                            gCCAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGt      <  1:27648/140‑1 (MQ=255)                                                                                                                              cAGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCg    <  1:897377/140‑1 (MQ=255)                                                                                                                               aGGGATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGc   <  1:1054089/140‑1 (MQ=255)                                                                                                                                gggATTGTCTTCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGct  >  1:579960/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                  GGCAAAAGTGCAGGAAGGCGTTGATATCGCCATCGAAAACGATGTGATTGTCGCTATCGGCGATGCCCTGACGCAACGCTACCCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCT  >  NC_000913/3019202‑3019468

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