Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,601,276	G→A	P271S (CCG→TCG)	ftsX ←	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,601,276	0	G	A	100.0%	162.6 / NA	42	P271S (CCG→TCG)	ftsX	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/17);  total (25/17)

ATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACG  >  NC_000913/3601129‑3601418                                                                                                                                                    |                                                                                                                                               aTGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAgcg                                                                                                                                             >  1:223243/1‑151 (MQ=255)   gATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAgcgc                                                                                                                                            >  2:567816/1‑150 (MQ=255)         ccATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGCTGACAATCGCAGCACCAGAATTTCTGACAACATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATg                                                                                                                                      >  1:88082/1‑150 (MQ=255)             tGATCTCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATcc                                                                                                                                 >  1:163565/1‑151 (MQ=255)                        ttCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCAc                                                                                                                      >  1:655158/1‑151 (MQ=255)                                aaaacCCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCg                                                                                                                         <  1:653879/139‑1 (MQ=255)                                 aaacccggccccACTTCGGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAAt                                                                                                                  <  1:644798/136‑1 (MQ=255)                                 aaacCCTGGCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCa                                                                                                                <  2:481398/148‑1 (MQ=255)                                 aaaaCCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGtt                                                                                                             >  1:334486/1‑151 (MQ=255)                                            ctctaccGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGtt                                                                                                   <  1:396720/145‑1 (MQ=255)                                              ctcccGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGt                                                                                                        <  2:463998/140‑1 (MQ=255)                                              cTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAAt                                                                                                >  2:143709/1‑151 (MQ=255)                                               ttCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTa                                                                                                  >  2:234644/1‑148 (MQ=255)                                                     aaCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTc                                                                                          >  1:488701/1‑150 (MQ=255)                                                           cgATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCg                                                                                    >  1:376908/1‑150 (MQ=255)                                                            gATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCg                                                                                    >  2:278204/1‑149 (MQ=255)                                                                  aaTCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCTCCATACAGGAACGAGCGCAGGATGAAGCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGGGCaac                                                                              <  1:173798/149‑2 (MQ=255)                                                                  aaTCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCaa                                                                               >  2:646387/1‑148 (MQ=255)                                                                     cGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATAc                                                                         <  2:16225/151‑1 (MQ=255)                                                                      gcagcaCCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACt                                                                        >  1:593162/1‑151 (MQ=255)                                                                          caCCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTc                                                                       >  1:333123/1‑148 (MQ=255)                                                                            ccAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGAca                                                                  <  1:146974/151‑2 (MQ=255)                                                                                aaTTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGcaca                                                              <  1:546272/151‑1 (MQ=255)                                                                                 aTTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCacac                                                             >  1:339583/1‑151 (MQ=255)                                                                                   ttCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTg                                                           >  1:633338/1‑151 (MQ=255)                                                                                        aCAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTAcc                                                      >  1:376083/1‑151 (MQ=255)                                                                                                    acaacaATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGaa                                          <  2:593167/151‑1 (MQ=255)                                                                                                         aaTGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACgg                                     >  2:27420/1‑151 (MQ=255)                                                                                                                      aTCCCAGCAGTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCa                         >  2:439382/1‑150 (MQ=255)                                                                                                                          cagcagTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACg                    <  1:298419/151‑1 (MQ=255)                                                                                                                           agcagTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGc                   <  2:345229/151‑1 (MQ=255)                                                                                                                              agTGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGa                <  1:381369/151‑1 (MQ=255)                                                                                                                                tGCGCCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGatc              >  1:479568/1‑151 (MQ=255)                                                                                                                                  cgcCACCATACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGatcat            <  2:389252/151‑1 (MQ=255)                                                                                                                                    ccaccaTACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCg          >  2:8557/1‑151 (MQ=255)                                                                                                                                     caccaTACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCAgcgc         <  1:228272/151‑1 (MQ=255)                                                                                                                                       ccaTACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCg       >  2:387374/1‑151 (MQ=255)                                                                                                                                       ccaTACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCg       >  2:656149/1‑151 (MQ=255)                                                                                                                                        caTACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCg       >  2:434578/1‑150 (MQ=255)                                                                                                                                         aTACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGaa     <  1:665413/151‑1 (MQ=255)                                                                                                                                          tACAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGaa     >  1:605223/1‑150 (MQ=255)                                                                                                                                            cAGGAACGAGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAAcg  <  1:278204/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               ATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGGAGTCACGGCGAGCAAAGATACTCAGACGCACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACG  >  NC_000913/3601129‑3601418

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