Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,398,446	G→A	L144F (CTC→TTC)	mreD ←	cell shape determining protein MreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,398,446	0	G	A	100.0%	150.9 / NA	46	L144F (CTC→TTC)	mreD	cell shape determining protein MreD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/27);  total (19/27)

AGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGC  >  NC_000913/3398309‑3398581                                                                                                                                          |                                                                                                                                        aGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACa                                                                                                                                       <  1:613579/140‑1 (MQ=255)   gTCACGCCAAGTTGCGCAAGTATCTTCTGGCGACGCGAAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACAcc                                                                                                                                     <  1:249491/140‑1 (MQ=255)      aCGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCa                                                                                                                                  <  1:633218/140‑1 (MQ=255)           aaGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGac                                                                                                                             <  2:721079/140‑1 (MQ=255)           aaGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGac                                                                                                                             <  1:378763/140‑1 (MQ=255)                gcgcAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGactacac                                                                                                                        <  2:582053/140‑1 (MQ=255)                 cgcAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGac                                                                                                                             >  1:92770/1‑134 (MQ=255)                 cgcAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGac                                                                                                                             <  2:92770/134‑1 (MQ=255)                          cTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGaa                                                                                                              <  1:733681/140‑1 (MQ=255)                                acgacgCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTc                                                                                                        >  2:71035/1‑140 (MQ=255)                                acgacgCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTc                                                                                                        <  1:707666/140‑1 (MQ=255)                                acgacgCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGACCACAACGCCCCATTGTCTACCCTACTCCAGAACACTTc                                                                                                        >  2:325269/1‑140 (MQ=255)                                  gacgCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCg                                                                                                      <  1:736955/140‑1 (MQ=255)                                                  aGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGtt                                                                                      <  1:185635/140‑1 (MQ=255)                                                  aGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGtt                                                                                      >  1:504450/1‑140 (MQ=255)                                                     gacagaCAGAGAAGTCACAGAACCCTTTATTGCACCGCACACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAATACTTCCGGTCTGAAAGAGACGTTAAt                                                                                   <  2:211861/135‑1 (MQ=255)                                                        aTACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAag                                                                                             >  2:490960/1‑127 (MQ=255)                                                        aTACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAag                                                                                             <  1:490960/127‑1 (MQ=255)                                                             gagaAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTaaaa                                                                           <  2:386373/140‑1 (MQ=255)                                                              agaAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTaaaaa                                                                          >  2:402112/1‑140 (MQ=255)                                                              agaAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTaaaaa                                                                          >  1:709336/1‑140 (MQ=255)                                                              agaAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTaaaaa                                                                          >  2:599578/1‑140 (MQ=255)                                                              agaAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTaaaaa                                                                          >  2:638836/1‑140 (MQ=255)                                                                 aGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAActc                                                                       >  2:376715/1‑140 (MQ=255)                                                                   tCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTg                                                                     >  1:730003/1‑140 (MQ=255)                                                                     aTAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAATACTTGCGGTCTGGAAGTGAGGTTAATAACTAAAAACTCTgtc                                                                   >  1:211861/1‑138 (MQ=255)                                                                       aGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCa                                                                 <  2:730003/140‑1 (MQ=255)                                                                       aGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCa                                                                 <  1:71035/140‑1 (MQ=255)                                                                          aaCCTTTATTGCACTTCAAACTGCTGACGGACTTTGCTCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCAGCCCAGaa                                                              <  1:256394/140‑1 (MQ=255)                                                                          aaCCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaa                                                              <  2:709336/140‑1 (MQ=255)                                                                          aaCCTTTATTGCACTGCAAACTGCTGACGGAATCTTCGCGTCAGCAAGTAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaa                                                              <  2:315700/140‑1 (MQ=255)                                                                            ccTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaaaa                                                            <  1:402112/140‑1 (MQ=255)                                                                            ccTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaaaa                                                            <  1:638836/140‑1 (MQ=255)                                                                            ccTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaaaa                                                            <  1:599578/140‑1 (MQ=255)                                                                                       cTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATcc                                                 >  2:136623/1‑140 (MQ=255)                                                                                          cAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATccacc                                              >  1:446861/1‑140 (MQ=255)                                                                                          cAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAAACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATccacc                                              >  1:53903/1‑140 (MQ=255)                                                                                               tgctgACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAg                                         >  2:389944/1‑140 (MQ=255)                                                                                                       gACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACACCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCa                                 <  1:77973/140‑1 (MQ=255)                                                                                                         cTTTGCGCATCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAAc                               >  1:683675/1‑140 (MQ=255)                                                                                                                  tCAGCAAGAAAATCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGAcc                      >  2:593856/1‑140 (MQ=255)                                                                                                                      cAAGAAAATCCACGGACAGAACACCCCATTGACTACAGTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAgcg                  <  1:325269/140‑1 (MQ=255)                                                                                                                          aaaaTCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCctg              <  2:504450/140‑1 (MQ=255)                                                                                                                          aaaaTCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCctg              <  2:683675/140‑1 (MQ=255)                                                                                                                           aaaTCCACGGCCAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCctgc             >  2:151784/1‑140 (MQ=255)                                                                                                                                      cAGAACACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGc  <  1:389944/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        AGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGC  >  NC_000913/3398309‑3398581

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