Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,019,938	A→G	T260A (ACG→GCG)	ssnA →	putative aminohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,019,938	0	A	G	100.0%	191.5 / NA	54	T260A (ACG→GCG)	ssnA	putative aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/31);  total (23/31)

GCCGTGAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTACAACCATCACCTTAGCGACAT  >  NC_000913/3019806‑3020071                                                                                                                                     |                                                                                                                                      gCCGTGAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGt                                                                                                                                <  1:592172/140‑1 (MQ=255)   cGTGAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCg                                                                                                                              <  2:635991/140‑1 (MQ=255)      gAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTc                                                                                                                           <  1:1020540/140‑1 (MQ=255)      gAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTc                                                                                                                           <  2:983013/140‑1 (MQ=255)      gAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTc                                                                                                                           <  2:3260/140‑1 (MQ=255)      gAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTc                                                                                                                           <  2:207218/140‑1 (MQ=255)       aaaGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCa                                                                                                                          <  2:341751/140‑1 (MQ=255)        aaGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCAt                                                                                                                         <  2:220078/140‑1 (MQ=255)           ccACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATggg                                                                                                                      <  2:330132/140‑1 (MQ=255)               aGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGt                                                                                                                  >  2:496669/1‑140 (MQ=255)                     tGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCGTGGGCTGTACTTGt                                                                                                            >  1:65030/1‑140 (MQ=255)                       gTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCg                                                                                                          <  2:406633/140‑1 (MQ=255)                            cATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGa                                                                                                     <  1:496669/140‑1 (MQ=255)                                 tCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACa                                                                                                >  2:743143/1‑140 (MQ=255)                                 tCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACa                                                                                                <  2:1121598/140‑1 (MQ=255)                                  cACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTa                                                                                                                 <  2:879136/122‑1 (MQ=255)                                  cACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTa                                                                                                                 >  1:879136/1‑122 (MQ=255)                                     gCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCGTGGGCTGTACTTGTCGAAAGATGACATCAc                                                                                            <  2:65030/140‑1 (MQ=255)                                     gCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCAc                                                                                            <  2:330258/140‑1 (MQ=255)                                     gCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCAc                                                                                            <  1:557682/140‑1 (MQ=255)                                     gCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCAc                                                                                            <  2:584703/140‑1 (MQ=255)                                           gAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACt                                                                                      >  1:509076/1‑140 (MQ=255)                                                 cTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCTATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCa                                                                                >  1:146974/1‑140 (MQ=255)                                                 cTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCa                                                                                >  1:923475/1‑140 (MQ=255)                                                    tACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAgcg                                                                             <  2:279825/140‑1 (MQ=255)                                                    tACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAgcg                                                                             <  2:285020/140‑1 (MQ=255)                                                           tttCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCg                                                                      >  2:1149525/1‑140 (MQ=255)                                                           tttCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCg                                                                      >  2:845987/1‑140 (MQ=255)                                                            ttCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGt                                                                     <  1:108428/140‑1 (MQ=255)                                                              ccTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTc                                                                   >  2:462860/1‑140 (MQ=255)                                                                 aCAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTg                                                                <  1:1149525/140‑1 (MQ=255)                                                                  cAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTgg                                                               <  1:1156108/140‑1 (MQ=255)                                                                   aGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGt                                                              >  2:1152148/1‑140 (MQ=255)                                                                       accacTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCAt                                                          >  2:484106/1‑140 (MQ=255)                                                                       accacTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCAt                                                          >  1:1135042/1‑140 (MQ=255)                                                                         cacTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGc                                                                       >  1:1017602/1‑125 (MQ=255)                                                                         cacTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGc                                                                       <  2:1017602/125‑1 (MQ=255)                                                                            tGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCTATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGc                                                     <  2:146974/140‑1 (MQ=255)                                                                            tGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGc                                                     <  2:923475/140‑1 (MQ=255)                                                                            tGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGc                                                     <  1:743143/140‑1 (MQ=255)                                                                             ggTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGcc                                                    >  2:715140/1‑140 (MQ=255)                                                                                   gCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCa                                              >  1:270086/1‑140 (MQ=255)                                                                                     aaaGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCaaa                                            >  1:544933/1‑140 (MQ=255)                                                                                      aaGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAAc                                           >  2:299913/1‑140 (MQ=255)                                                                                      aaGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAAc                                           <  1:1152148/140‑1 (MQ=255)                                                                                      aaGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAAc                                           >  1:662880/1‑140 (MQ=255)                                                                                                      aCTGGCGCAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCg                           >  2:945119/1‑140 (MQ=255)                                                                                                          gcgcAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTa                       <  2:662880/140‑1 (MQ=255)                                                                                                          gcgcAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTa                       >  2:440250/1‑140 (MQ=255)                                                                                                          gcgcAATTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTa                       <  1:484106/140‑1 (MQ=255)                                                                                                               aTTCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTACAAcc                  <  1:945119/140‑1 (MQ=255)                                                                                                                 tCGATCTCATCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTACAGCCAt                <  2:1135042/140‑1 (MQ=255)                                                                                                                         aTCGACAGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTACAACCATCACCTTAg        >  2:225204/1‑140 (MQ=255)                                                                                                                               aGCAAAGCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTACAACCATCACCTTAGCGACAt  >  1:628537/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                      GCCGTGAAAGCCACAGGCCGTGGTTTGCATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCATGTCGGCTACAACCATCACCTTAGCGACAT  >  NC_000913/3019806‑3020071

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: ‑