Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,316,957	Δ1 bp	coding (262/1137 nt)	phnM ←	RPnTP hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,316,956	0	C	.	100.0%	245.5 / NA	52	coding (263/1137 nt)	phnM	RPnTP hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (18/34);  total (18/34)

TCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAA  >  NC_000913/4316876‑4317032                                                                                 |                                                                             tcatcTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGc                                                                              >  2:861751/1‑81 (MQ=255)       tctcCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGc                                                                       >  1:861846/1‑81 (MQ=255)        ctcCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCg                                                                        <  1:861752/79‑1 (MQ=255)           cAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGccac                                                                    >  1:861829/1‑80 (MQ=255)           cAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGcca                                                                     <  1:861754/79‑1 (MQ=255)             gATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCAt                                                                 <  2:861755/81‑1 (MQ=255)             gATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCAt                                                                 <  1:861756/81‑1 (MQ=255)              aTTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATc                                                                <  1:861757/81‑1 (MQ=255)               ttCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCAt                                                                 >  2:861837/1‑79 (MQ=255)                tctcCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAg                                                              >  1:861825/1‑81 (MQ=255)                     aGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGcgt                                                         >  2:861852/1‑81 (MQ=255)                          aTCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGACG‑CGCTCGCCACCATGAGCGCGTCGTgg                                                    >  2:861827/1‑81 (MQ=255)                            cgccgcCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgct                                                  >  2:861845/1‑81 (MQ=255)                             gccgccGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctg                                                 >  1:861826/1‑81 (MQ=255)                             gccgccGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgct                                                  <  2:861758/80‑1 (MQ=255)                             gccgccGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgct                                                  <  1:861759/80‑1 (MQ=255)                              ccgccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctgc                                                <  2:861761/81‑1 (MQ=255)                              ccgccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctgc                                                <  1:861760/81‑1 (MQ=255)                              ccgccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctgc                                                >  2:861833/1‑81 (MQ=255)                               cgccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctgct                                               >  2:861848/1‑81 (MQ=255)                               cgccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctgct                                               >  1:861830/1‑81 (MQ=255)                                gccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctgct                                               >  2:861844/1‑80 (MQ=255)                                gccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGgctgct                                               >  1:861859/1‑80 (MQ=255)                                gccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTc                                              <  1:861762/81‑1 (MQ=255)                                 ccgTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCa                                             <  1:861763/81‑1 (MQ=255)                                   gTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATc                                           <  1:861765/81‑1 (MQ=255)                                   gTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATc                                           <  2:861764/81‑1 (MQ=255)                                      gcgcACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGcc                                        <  1:861766/81‑1 (MQ=255)                                       cgcACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCg                                       <  1:861767/81‑1 (MQ=255)                                         cACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAg                                     >  2:861841/1‑81 (MQ=255)                                               gCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTggg                                 >  2:861838/1‑79 (MQ=255)                                               gCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGc                                <  1:861769/80‑1 (MQ=255)                                               gCCAATTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCa                               <  1:861768/81‑1 (MQ=255)                                                  aaTTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGc                            >  1:861839/1‑81 (MQ=255)                                                  aaTTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGc                            >  1:861855/1‑81 (MQ=255)                                                   aTTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGcc                           <  2:861770/81‑1 (MQ=255)                                                   aTTGCCACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGcc                           <  2:861771/81‑1 (MQ=255)                                                       ccACGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTc                       <  1:861772/81‑1 (MQ=255)                                                         aCGGCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCaa                     <  1:861773/81‑1 (MQ=255)                                                           ggCATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACt                   <  1:861774/81‑1 (MQ=255)                                                             cATCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACttt                 >  1:861840/1‑81 (MQ=255)                                                              aTCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTc                <  1:861777/81‑1 (MQ=255)                                                              aTCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTc                <  2:861776/81‑1 (MQ=255)                                                              aTCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTc                <  2:861778/81‑1 (MQ=255)                                                              aTCCAGTACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTAAACTTTc                <  2:861775/81‑1 (MQ=255)                                                                    tACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGgcgc          <  2:861779/81‑1 (MQ=255)                                                                    tACGGTGGTGATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGgcgc          <  2:861780/81‑1 (MQ=255)                                                                           gtgATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTga   <  1:861782/81‑1 (MQ=255)                                                                           gtgATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTg    <  1:861785/80‑1 (MQ=255)                                                                           gtgATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTg    <  2:861783/80‑1 (MQ=255)                                                                           gtgATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTg    <  2:861784/80‑1 (MQ=255)                                                                           gtgATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTg    <  2:861781/80‑1 (MQ=255)                                                                            tgATG‑CGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTgaa  <  1:861786/81‑1 (MQ=255)                                                                                 |                                                                             TCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAA  >  NC_000913/4316876‑4317032

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