Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,530,925	T→C	Q103Q (CAA→CAG)	zipA ←	FtsZ stabilizer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,530,925	0	T	C	98.2%	194.4 / NA	55	Q103Q (CAA→CAG)	zipA	FtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/25);  minor base G (1/0);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATG  >  NC_000913/2530830‑2531017                                                                                                |                                                                                             tCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAc                                                                                          >  2:445078/1‑100 (MQ=255)  cAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAcg                                                                                         >  1:332146/1‑100 (MQ=255)  cAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAcg                                                                                         >  2:387926/1‑100 (MQ=255)  cAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCGGCGCTTCAGGCGGCTGCTGGACCGGCTGAc                                                                                          >  2:400195/1‑99 (MQ=255)    ggCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAcgcg                                                                                       >  1:560623/1‑100 (MQ=255)        gCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTg                                                                                   >  1:813235/1‑100 (MQ=255)        gCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTg                                                                                   <  2:560806/100‑1 (MQ=255)         cGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTgc                                                                                  >  2:91837/1‑100 (MQ=255)           gCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTgcgc                                                                                <  2:888784/100‑1 (MQ=255)               aTAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGAc                                                                            >  1:1150685/1‑100 (MQ=255)                tAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACg                                                                           >  2:549289/1‑100 (MQ=255)                 aggcaggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGc                                                                          >  1:1157440/1‑100 (MQ=255)                 aggcaggcCGCTGCCACGGCTGCGCTGGATACGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACg                                                                           <  1:1419849/99‑1 (MQ=255)                    caggcTGCTGCACCGGCTGCGCTGGATGTGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATa                                                                       <  2:19798/100‑1 (MQ=255)                            tgcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGtt                                                               <  2:936956/100‑1 (MQ=255)                             gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTg                                                              <  1:177096/100‑1 (MQ=255)                             gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTg                                                              >  2:1303807/1‑100 (MQ=255)                             gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTg                                                              >  1:92440/1‑100 (MQ=255)                               aCCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGt                                                            <  1:445079/100‑1 (MQ=255)                                    cTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTgg                                                       >  1:1018742/1‑100 (MQ=255)                                    cTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTgg                                                       <  1:1336896/100‑1 (MQ=255)                                     tGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGgt                                                      >  2:396744/1‑100 (MQ=255)                                      gcgcTGGATGCTGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGgtg                                                     <  2:718911/100‑1 (MQ=255)                                       cgcTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGgtgt                                                    >  1:792068/1‑100 (MQ=255)                                             aTGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggc                                              <  1:1432065/100‑1 (MQ=255)                                              tGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggcg                                             <  2:1467479/100‑1 (MQ=255)                                              tGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggcg                                             <  2:1311818/100‑1 (MQ=255)                                              tGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggcg                                             <  1:370827/100‑1 (MQ=255)                                                cGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgac                                           <  1:1093581/100‑1 (MQ=255)                                                cGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgac                                           >  2:717006/1‑100 (MQ=255)                                                cGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgac                                           >  1:1502701/1‑100 (MQ=255)                                                cGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgac                                           >  1:31761/1‑100 (MQ=255)                                                  gAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacgg                                         <  1:943665/100‑1 (MQ=255)                                                     cATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacggacg                                      >  1:1000129/1‑100 (MQ=255)                                                       tGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGGTGACGCGGCTGCGCAGATGCATAAGGCGGTAGGTACTGGTGGTgcgg                                               >  2:1225789/1‑89 (MQ=255)                                                        gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                   >  1:166474/1‑100 (MQ=255)                                                         ttGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagca                                  >  2:1182933/1‑100 (MQ=255)                                                            cggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGcc                               <  1:1409009/100‑1 (MQ=255)                                                               cggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCa                            <  2:279304/100‑1 (MQ=255)                                                               cggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCa                            <  1:1565527/100‑1 (MQ=255)                                                                 gTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATg                          <  2:345614/100‑1 (MQ=255)                                                                   aCCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCt                        >  1:719987/1‑100 (MQ=255)                                                                     cTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTcc                      <  2:332146/100‑1 (MQ=255)                                                                      tGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCt                     >  1:633422/1‑100 (MQ=255)                                                                       gcgcTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTg                    >  2:864347/1‑100 (MQ=255)                                                                       gcgcTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTg                    <  2:813230/100‑1 (MQ=255)                                                                        cgcTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGa                   <  2:1074553/100‑1 (MQ=255)                                                                           ttCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCg                <  1:653691/100‑1 (MQ=255)                                                                             naGGCGGCTGCTGGACCGGCTGACGCGGCTACGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCGCCTGAGCGtt              >  2:1037834/2‑100 (MQ=255)                                                                               ggcggcTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAg            >  1:871404/1‑100 (MQ=255)                                                                                 cggcTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGcc          >  2:1050684/1‑100 (MQ=255)                                                                                      gctgGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGc     <  1:396745/100‑1 (MQ=255)                                                                                      gctgGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGc     >  2:1134746/1‑100 (MQ=255)                                                                                       ctgGACCGGCTGACGCGGCTGCCCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCa    >  1:946976/1‑100 (MQ=255)                                                                                         gggccGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGGATg  <  2:438524/97‑1 (MQ=255)                                                                                                |                                                                                             TCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATG  >  NC_000913/2530830‑2531017

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