Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	104,426	C→T	L149F (CTT→TTT)	ftsA →	cell division protein FtsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	104,426	0	C	T	100.0%	125.4 / NA	33	L149F (CTT→TTT)	ftsA	cell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/20);  total (13/20)

GATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCAAAGCG  >  NC_000913/104348‑104503                                                                               |                                                                              gATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGAttt                                                                             >  2:22702/1‑81 (MQ=255)      gCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCggg                                                                        >  2:22703/1‑81 (MQ=255)           gtgtGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTg                                                                    <  1:22646/80‑1 (MQ=255)           gtgtGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGc                                                                   <  1:22645/81‑1 (MQ=255)             gtgCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCg                                                                  >  1:22746/1‑80 (MQ=255)                   cATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAg                                                           <  1:22647/81‑1 (MQ=255)                       tgATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCa                                                        <  2:22648/80‑1 (MQ=255)                         aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCaaaa                                                     <  1:22649/81‑1 (MQ=255)                         aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCaaaa                                                     <  2:22650/81‑1 (MQ=255)                            ccGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTg                                                  <  2:22651/81‑1 (MQ=255)                              gCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCa                                                >  1:22707/1‑81 (MQ=255)                                aaGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAcc                                              <  1:22652/81‑1 (MQ=255)                                 agagTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCt                                             <  2:22653/81‑1 (MQ=255)                                  gagTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                                            <  2:22654/81‑1 (MQ=255)                                     tATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATc                                         <  2:22655/81‑1 (MQ=255)                                     tATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATc                                         <  2:22656/81‑1 (MQ=255)                                         cGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACAt                                     >  2:22774/1‑81 (MQ=255)                                            ttGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGt                                   >  2:22693/1‑80 (MQ=255)                                            ttGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                                  <  1:22657/81‑1 (MQ=255)                                            ttGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                                  <  2:22658/81‑1 (MQ=255)                                              gACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTca                                 >  1:22695/1‑80 (MQ=255)                                              gACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTca                                 >  2:22698/1‑80 (MQ=255)                                              gACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTca                                 <  1:22659/80‑1 (MQ=255)                                                     aGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGata                         >  2:22709/1‑81 (MQ=255)                                                       gAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATg                       >  2:564859/1‑81 (MQ=255)                                                        aaGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATgg                      <  1:22660/81‑1 (MQ=255)                                                        aaGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATgg                      >  2:22708/1‑81 (MQ=255)                                                          gggATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGc                     >  2:22700/1‑80 (MQ=255)                                                                aaGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAAc              <  1:22661/81‑1 (MQ=255)                                                                  gAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACAt            >  2:22726/1‑81 (MQ=255)                                                                    aTCCGGTAGGATTTTCGGGAGTGAGGATGCAGGCAAAAGTGCACCTAATCACATGTCACAACGATATGGCGAAAAACATCg          <  1:22662/81‑1 (MQ=255)                                                                       cGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCa       <  1:22663/81‑1 (MQ=255)                                                                          tAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCAAAg    <  2:22664/81‑1 (MQ=255)                                                                            ggATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCAAAGCg  >  1:22723/1‑81 (MQ=255)                                                                               |                                                                              GATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCAAAGCG  >  NC_000913/104348‑104503

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