Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,020,348	C→A	G396G (GGC→GGA)	ssnA →	putative aminohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,020,348	0	C	A	100.0%	163.7 / NA	47	G396G (GGC→GGA)	ssnA	putative aminohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/29);  total (18/29)

ATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGA  >  NC_000913/3020214‑3020476                                                                                                                                       |                                                                                                                                 atgaGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAgg                                                                                                                             >  1:178833/1‑140 (MQ=255)        gCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGtt                                                                                                                      <  2:61672/140‑1 (MQ=255)         cAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTc                                                                                                                     >  1:356467/1‑140 (MQ=255)          aaCTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTca                                                                                                                    <  2:1148595/140‑1 (MQ=255)                ggCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGt                                                                                                              >  1:723470/1‑140 (MQ=255)                 gcgcGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTg                                                                                                             <  2:849252/140‑1 (MQ=255)                 gcgcGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTg                                                                                                             <  2:382391/140‑1 (MQ=255)                  cgcgAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGa                                                                                                            <  2:449212/140‑1 (MQ=255)                      aaaTTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGt                                                                                                        <  2:603948/140‑1 (MQ=255)                      aaaTTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGt                                                                                                        <  2:1091698/140‑1 (MQ=255)                       aaTTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTc                                                                                                       <  1:153687/140‑1 (MQ=255)                         tttGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCaa                                                                                                     <  1:1081821/140‑1 (MQ=255)                                ttCTGTAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACACCTCGTGGACGCCGCTGCTGGCAGACAATATCGCCGGGCTTATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTTAATGGtgtg                                                                                              <  2:835302/140‑1 (MQ=255)                                ttCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGtgtg                                                                                              <  1:1145600/140‑1 (MQ=255)                                ttCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGtgtg                                                                                              <  2:550064/140‑1 (MQ=255)                                ttCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGtgtg                                                                                              <  2:303076/140‑1 (MQ=255)                                  cTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGTTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGAt                                                                                            >  2:89935/1‑140 (MQ=255)                                  cTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGAt                                                                                            <  2:586168/140‑1 (MQ=255)                                  cTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGAt                                                                                            <  2:1101053/140‑1 (MQ=255)                                  cTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGAt                                                                                            <  2:1377222/140‑1 (MQ=255)                                    ggAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATgg                                                                                          <  2:156441/140‑1 (MQ=255)                                    ggAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATgg                                                                                          <  2:128241/140‑1 (MQ=255)                                          cGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATg                                                                                    >  1:1102106/1‑140 (MQ=255)                                            gTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGaa                                                                                  >  1:249336/1‑140 (MQ=255)                                            gTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGaa                                                                                  >  2:1303024/1‑140 (MQ=255)                                            gTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGaa                                                                                  >  2:1202201/1‑140 (MQ=255)                                            gTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGaa                                                                                  >  1:282673/1‑140 (MQ=255)                                                  aaGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGTTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGt                                                                            <  1:89935/140‑1 (MQ=255)                                                  aaGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGt                                                                            <  2:1143591/140‑1 (MQ=255)                                                    gCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGt                                                                            >  1:1317405/1‑138 (MQ=255)                                                    gCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGt                                                                            <  2:1317405/138‑1 (MQ=255)                                                        aTTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGttt                                                                      <  1:395732/140‑1 (MQ=255)                                                             aCCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACtt                                                                 >  2:1395176/1‑140 (MQ=255)                                                                  ttGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGAtt                                                            >  1:165654/1‑140 (MQ=255)                                                                  ttGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGAtt                                                            >  1:768127/1‑140 (MQ=255)                                                                       aTTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGAt                                                       <  2:723470/140‑1 (MQ=255)                                                                         tACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTc                                                     <  2:768127/140‑1 (MQ=255)                                                                                                 gCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGc                             <  1:1303024/140‑1 (MQ=255)                                                                                                   aGACAATATCGCCGGGCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCg                           <  2:249336/140‑1 (MQ=255)                                                                                                              ccGGGCATAGCGCTTTCGGTATGGAAACAGGCAGCGTTTACAGCGAGATGGTCACTGGCGCGATGGACTATGAAGACCGCCAATTTAACTTCGCTTGCGATTCCATTTATGCACAAGCCAGCAAAGTCGCTCCCAGTAtg                >  2:1123269/1‑140 (MQ=255)                                                                                                                gggCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTAtgtg              <  1:1328615/140‑1 (MQ=255)                                                                                                                 ggCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTgg             >  2:1325835/1‑140 (MQ=255)                                                                                                                 ggCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGACAGTATGTgg             >  1:194220/1‑140 (MQ=255)                                                                                                                  gCATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGc            >  1:105803/1‑140 (MQ=255)                                                                                                                   cATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGGCGGTCAGTTTAGCTTCGATTGCGATTCCATTTATGGACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCg           <  2:783350/140‑1 (MQ=255)                                                                                                                   cATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCg           <  2:1175690/140‑1 (MQ=255)                                                                                                                   cATATCGCTTTCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCg           <  2:282673/140‑1 (MQ=255)                                                                                                                           tttCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCggatgg   >  1:810617/1‑140 (MQ=255)                                                                                                                            ttCGGTATGGGATCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCggatgga  >  1:560/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                 ATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGA  >  NC_000913/3020214‑3020476

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