Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,125,097	C→T	E572K (GAG→AAG)	priA ←	primosome factor N'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,125,097	0	C	T	97.7%	118.6 / NA	44	E572K (GAG→AAG)	priA	primosome factor N'
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/22);  total (21/23)
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

GCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGC  >  NC_000913/4124958‑4125235                                                                                                                                            |                                                                                                                                           gCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGc                                                                                                                                         <  1:502044/143‑1 (MQ=255)     aGGCGTCGTAGCCTGTAGAGAGCAACGGTTGCAACAGAGGATGTTCCGGATGGTGCGTGTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAg                                                                                                                                     <  2:499955/143‑1 (MQ=255)            gtAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAAt                                                                                                                              >  1:562727/1‑143 (MQ=255)               gCCTTTATAGGGCAACGTTTGCAATAGAGCATGTTCCGGCTGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCgg                                                                                                                           <  2:679877/143‑1 (MQ=255)               gCCTCTATAGAGCAACGTTTGCAACCGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCgg                                                                                                                           <  1:426404/143‑1 (MQ=255)                ccTTGATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGc                                                                                                                          <  1:357382/143‑1 (MQ=255)                  tttATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAg                                                                                                                        >  1:764904/1‑143 (MQ=255)                           cAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGCAAATCGGCAGAAAACGgcg                                                                                                               >  2:171847/1‑143 (MQ=255)                           cAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGACGAGCGGAAATAGGCAGCAAACAgcg                                                                                                               >  1:499383/1‑143 (MQ=255)                           cAACGTTTGCAACAGAGGATGTGCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAgcg                                                                                                               <  1:534296/143‑1 (MQ=255)                                  tGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtc                                                                                                        <  1:81517/143‑1 (MQ=255)                                    cAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAACACCGGCAGAAACCAGCGCGCcgtcca                                                                                                      >  2:177830/1‑143 (MQ=255)                                      aCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCGCGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGACAACCGGCAGAAAACAGCGCGCcgtccacg                                                                                                    >  2:152506/1‑143 (MQ=255)                                          agGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCTGCGACCTGGGTGTACAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccacgtcca                                                                                                >  1:481872/1‑143 (MQ=255)                                            gATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGt                                                                                              >  1:369467/1‑143 (MQ=255)                                              tGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTaa                                                                                            <  1:322126/143‑1 (MQ=255)                                               gtgccGTATGGTGCGTTTGCTGCCCCACTGCGCCCTGTTCACCCCCACGCCCGGCACGACCGGCGACCTGGGTGTCAAGGTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAAt                                                                                           <  2:86008/140‑1 (MQ=255)                                               gTTCCGGATGGTGCGTTTGCAGCCCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAAt                                                                                           <  1:591306/143‑1 (MQ=255)                                                          cacGTTTGCAGCGCACCTTTGTCCTCTTACCCCGCACGCCCGGCGCGACCGTCGACCTCGGTGCAAAGCTGAGCGAAACGCTTTGCCGCGCGAAAATCTGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGt                                                                                <  2:72415/141‑1 (MQ=255)                                                           gCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAACACATCGTGCCGTCCTCGTCCCGTCCTGCAACCACCGTc                                                                               >  1:11173/1‑143 (MQ=255)                                                                 gcagcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATcc                                                                         >  1:151183/1‑143 (MQ=255)                                                                   agtgccacTTCGCCCTGTTCCCCCGCCCGCCCGGCACGAACGGCGACCTGGGTGTAACGCTGTGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCgg                                                                       <  2:45043/139‑1 (MQ=255)                                                                   agcaCCACTTCGCCCCGTTTAACCGCACGCCCGGCACGACCGGCGACCGGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAACGAAAATCGGCAGAAAACAGCGCGCCGTAGCAGACCAGTAATGCCACCAGCGTcacaacagg                                                                       >  2:212687/1‑138 (MQ=255)                                                                      accacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCGTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCCGCGTCACATCCgggga                                                                    >  2:114699/1‑141 (MQ=255)                                                                           ttCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGt                                                               <  1:623178/143‑1 (MQ=255)                                                                               cccTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGACAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGGCACATCCGGGCAAGGGTGAcc                                                           >  2:715618/1‑143 (MQ=255)                                                                               cccTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGAcc                                                           >  1:22711/1‑143 (MQ=255)                                                                                   gTTTACCCGCGCGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTc                                                       <  1:422990/143‑1 (MQ=255)                                                                                     ttACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCCACCAGCGTCACATCCGGGAAATGGTGCACTTTCGc                                                     >  2:45867/1‑143 (MQ=255)                                                                                         ccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATg                                                                 <  1:346446/127‑1 (MQ=255)                                                                                         ccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATg                                                                 >  2:346446/1‑127 (MQ=255)                                                                                             gcGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCAttt                                             <  2:369467/142‑1 (MQ=255)                                                                                               gCCCGGCACGACCGGCGACCTGGGTCTAAAGCTGAGCGAATCGCTTTGTCGGGCGAAAATCGGCAGAAAACCGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTtgt                                           <  1:756295/143‑1 (MQ=255)                                                                                                   ggCACGACCGGCGTCCTGGGTCTAAAGGTGATCGACACGCTTAGCCGAGCGAAAATCGGCAGAAAACAGCTCGCCGTCCACGTGCACTAATGCAACCAGCGGCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAc                                       <  2:45535/143‑1 (MQ=255)                                                                                                     cACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCGTTGGGGTagcg                                     >  1:150221/1‑140 (MQ=255)                                                                                                       cggccGGGGACCTGGGCGTGAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAGCAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAAt                                   <  2:679135/140‑1 (MQ=255)                                                                                                          ccGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCACCCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCaa                                <  2:170007/143‑1 (MQ=255)                                                                                                           cGGCGACCTGGGTGTAAAGCTGAGCGAACCGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGACAGCGTTTGGGGACCaagaaga                               >  1:267443/1‑138 (MQ=255)                                                                                                           cGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTGGTGTACCAATCaaa                               >  1:248735/1‑143 (MQ=255)                                                                                                                          aaaGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACGAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAGCAATCAAACTCCGCGcgccgccg                >  1:28123/1‑143 (MQ=255)                                                                                                                                  gCGTAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGTGCCGTCCACGTCCATTAGTGCAACCAGCGTCACATCAGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTAc        <  2:577131/143‑1 (MQ=255)                                                                                                                                   cGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACt       >  1:8459/1‑143 (MQ=255)                                                                                                                                     aaaCCCGCTGCCGAGCGAAAAGCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATGTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTc     <  2:684493/143‑1 (MQ=255)                                                                                                                                        cGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGc  <  1:329021/143‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           GCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGC  >  NC_000913/4124958‑4125235

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