Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,123,110	G→A	P77L (CCA→CTA)	ftsN ←	cell division protein FtsN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,123,110	0	G	A	97.6%	111.8 / NA	42	P77L (CCA→CTA)	ftsN	cell division protein FtsN
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (24/17);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAG  >  NC_000913/4122970‑4123242                                                                                                                                             |                                                                                                                                     gTTCAAGCAGCTGACGTTGGTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtagt                                                                                                                                    <  2:1069229/143‑1 (MQ=255)   tCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtagtgg                                                                                                                                  <  1:1114237/143‑1 (MQ=255)     aagaagCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTa                                                                                                                                >  2:432096/1‑143 (MQ=255)            tGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACGACCGGCTGGAGGTTCTGTGGGCGCGGGGGCTCTCTGCTGGCGCCTTGCTGGCTGTTTTATGGGGCGCCAGCGTTCTCCTGGTTTTAGTGTTGGTCCGtt                                                                                                                         >  2:160906/1‑143 (MQ=255)                gttgttCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTTGCGACTGCCCAGCTCTCTAATGTAGCGCCAGCGTGTTTCTCCTTTTAGTGATAGTCCGTTTCCg                                                                                                                     >  2:1033035/1‑143 (MQ=255)                    ttCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTACCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCGGCGCGCTGCTGGTTTTAGTGGTCGTCCGTTGcccggcc                                                                                                                 >  2:833067/1‑138 (MQ=255)                    gtCTGGTGTCAGTTGCTCCGGCGTGTTCACTGCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTGTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCa                                                                                                                 <  2:1129333/142‑1 (MQ=255)                          tgtCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCt                                                                                                           <  1:1056068/143‑1 (MQ=255)                           gtCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTg                                                                                                          <  2:514701/143‑1 (MQ=255)                           gtCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTg                                                                                                          >  1:514701/1‑143 (MQ=255)                            tCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCTGTCACTTTCTgg                                                                                                         >  1:534328/1‑143 (MQ=255)                                 tGCTCTGCCGTTATCATTTCTCTACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTg                                                                                                    <  2:961710/143‑1 (MQ=255)                                           ttttCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCTtctct                                                                                          >  1:329676/1‑142 (MQ=255)                                           ttttCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCg                                                                                          >  1:485363/1‑143 (MQ=255)                                             ttcacttcacCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTGTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGa                                                                                        >  1:103760/1‑143 (MQ=255)                                                acttcacCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGCTTTAGTGGTAGTCCGTTTCCGGTCACCTTCTGGCTTTGCAGCGtctctgccca                                                                                     >  2:956250/1‑137 (MQ=255)                                                 cttcacCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGACAGCGTTCTTCTGGGTTTAGTGGTAGTCCGTTTACGGACACTGTCTGGCCTGGCAGCGTCTGGGActgt                                                                                    >  1:216503/1‑141 (MQ=255)                                                   tcacCACCGGCAGATGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACtcttc                                                                                  <  1:838041/143‑1 (MQ=255)                                                    cctctgcGGCCGAGGTTTCTGGGGGCGGACGGTCTCCCGTGTGGGGACTTTCCAGCTCTTTGGTGTGGCGGCGGCGGTCTTCTGGTTTTAGTGGTAGTCCGTGTCCGGTCACTTGCTGGCTGTGCAGCGTCTCGGACtcttct                                                                                 <  2:704909/137‑1 (MQ=255)                                                       caccGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCGTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACtcttggttg                                                                              >  1:480652/1‑138 (MQ=255)                                                         ccGGCAGAAGGTTCTGTGGGCGCACGCACGCCCGGCTGGTGTCTTTCCAGCTCTGTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTtctttctt                                                                            <  2:869830/143‑1 (MQ=255)                                                           ggccggaGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTACTGTCGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTTACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTtgt                                                                          <  2:303259/137‑1 (MQ=255)                                                                aaGGTTCTGTGGGCGCACGCACGCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGc                                                                     <  1:432096/143‑1 (MQ=255)                                                                       tgtgGTCACACGCGCTCCCGGCTGGCGTCTTTCCAGCTCTTTAATGTAGCGCCATCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGa                                                              <  2:889719/143‑1 (MQ=255)                                                                        gtgGGCGCACGCACTCTCGGCAGGCGAATTGCCAGCTCTGTAATGTAGCGCCAGCGTTCTTCCGGTTTTAGGGGTAGTCCTTTTCCGGTCACTTTCTGGCTTTGCCGCGTCACGCACTCTCACTGCGTGTGAAGCGTatcgga                                                             >  2:692280/1‑138 (MQ=255)                                                                          gggCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGt                                                           <  1:918392/143‑1 (MQ=255)                                                                           ggCGCACGCGCGCCCGTCTGGCGACTTGCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTGTAGTGCTAGTCCGTTACCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTa                                                          <  2:17678/143‑1 (MQ=255)                                                                             cgcacggcCCCCCGGCTGGCGACTTTCCAGCTCTTTGGTGTGGCGCCAGCGTTCTTCGGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGGCTCGGACTCTTCTTTCCTGTGATGCGTAATGAAGTACa                                                        <  2:399262/143‑1 (MQ=255)                                                                             cgcacgcaCTCCCGGCTGGCGACTCTCCAGCGCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACa                                                        <  2:745344/143‑1 (MQ=255)                                                                                acgcacTCCCGGCTGGCGACTTTCCAGCTCTTTATTGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGac                                                     <  1:710327/143‑1 (MQ=255)                                                                                acgcacTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTCGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGac                                                     >  1:909789/1‑143 (MQ=255)                                                                                acgcacTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGTTGCGTACTGAAGTACAgcc                                                     >  2:191829/1‑141 (MQ=255)                                                                                                    tttCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTCTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCCCGGAGTCTTCTCGCTTGTGCTGCGCAGTGAAGTGCAGACCCCCGCTATAGGGcacacgc                                 >  2:335764/1‑140 (MQ=255)                                                                                                      tCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTCCAGACAACCGATAAAGGTCACaacaac                               >  1:185596/1‑143 (MQ=255)                                                                                                      tCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAAc                               <  1:1199918/143‑1 (MQ=255)                                                                                                       cctgCTCTTGAATGTAGCTCCAGCGTGCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTGGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAAcg                              <  2:239419/140‑1 (MQ=255)                                                                                                        cAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACGACCGATGAAGGTCACAAGAAcgg                             >  1:677365/1‑143 (MQ=255)                                                                                                         aGCTCTTTAATGTAGCGCCAGCGTTCTTCNGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTCCAGGCCACCTATACAGGTCGCAAGAAcggc                            >  1:1190420/1‑143 (MQ=255)                                                                                                              tttAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAGTGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGgcag                       >  1:550399/1‑143 (MQ=255)                                                                                                                  aTGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTACTGATGTACAGACCACCGATAAAGCTCACAATACcggcgggaccagg                   >  1:924311/1‑136 (MQ=255)                                                                                                                  aTGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGACGGCAGCAAt                   >  1:1045382/1‑143 (MQ=255)                                                                                                                          ccAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGCCCACCGATGAAGGTCACACGAACGGCGGCAGCAACGGCGGCCa           >  2:399241/1‑143 (MQ=255)                                                                                                                               gTTCGTTTGGTTTTGGCGTGAGTCCGTTTCGGGTCACTTGCTGGCTTTGCGGCGTCTCGGACTCTGCTTTCTTGTGATGTGTAATGAAGTAGAGACCACCGAGAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCg      <  2:106718/143‑1 (MQ=255)                                                                                                                                   ttctGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGgag  >  1:1096939/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                     GTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAG  >  NC_000913/4122970‑4123242

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