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.9%	132.5 / NA	49	E572K (GAG→AAG)	priA	primosome factor N'
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (23/24);  minor base G (0/1);  total (24/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 = 9.99e-01

CAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCT  >  NC_000913/4124959‑4125233                                                                                                                                           |                                                                                                                                         cAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGGTGGTCCGGATGTTGAGTTTGCAGCACCACTTCGCCCTGCTTACCCGCACGCCCGGCACGACCGGCGCCCGGAGTGTAAAGCAGAGCGAGACGCTCTGct                                                                                                                                      >  2:626601/1‑142 (MQ=255)        gtcgtAGCCTTTTTAGAGCATCGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGaa                                                                                                                               <  2:89724/143‑1 (MQ=255)          ctgaGCCTCTATAGAGTTACGTTTGCAAGAGCGAACAATACGGATATTGCTTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGaaaa                                                                                                                             <  2:232058/140‑1 (MQ=255)            tAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATc                                                                                                                           <  2:741169/143‑1 (MQ=255)              gCCTTTATAGGGCCACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCgg                                                                                                                         <  1:452907/143‑1 (MQ=255)              gCCTTTACAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCCCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGCAAATCgg                                                                                                                         <  1:430135/143‑1 (MQ=255)              gCCTGTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCTTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCgg                                                                                                                         <  1:367690/143‑1 (MQ=255)                 tttATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGc                                                                                                                                       <  2:567134/126‑1 (MQ=255)                 tttATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGc                                                                                                                                       >  1:567134/1‑126 (MQ=255)                   tatcGAGCAACGTTTGCAACAGAGGGCGTTCCGGATGGTGCGTTTGCAGCGCCACTTCGCCCTGTTTACCCGCACGCCCGTCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGaa                                                                                                                    <  2:466675/143‑1 (MQ=255)                     taGAGCAACGTGTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGaaaa                                                                                                                  <  1:326203/143‑1 (MQ=255)                        agCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAg                                                                                                               <  1:187557/143‑1 (MQ=255)                         gCAACGTTTGCAACAGTGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTGTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGTTTTGCCGAGCGAAAATCGGCAGAAAACAgc                                                                                                              <  2:47352/143‑1 (MQ=255)                          cAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGGCCGGCACGACCGGCGACCTGGGTGTAAAGCGGACCGAACCGCTTTGCCGAGCCAAAACACGCGGACAACAgcg                                                                                                             >  1:1024863/1‑143 (MQ=255)                                ttGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGACACGCTTTGCCGAGCGAAAGTCGGCAGAAAACAGCGCGCcgt                                                                                                       >  1:464274/1‑143 (MQ=255)                                     aCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGCCCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGCGCGAAAATCGGCGGAAAACAGTGCGCcgtccccg                                                                                                  >  1:541699/1‑143 (MQ=255)                                       agagGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCAGGCACGACCGGCGCCCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAAGCGGCAGACAACAGCGCGCagtccaccgc                                                                                                >  2:821247/1‑140 (MQ=255)                                               ttCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGCGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATg                                                                                        >  1:637755/1‑143 (MQ=255)                                                        gTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCCACCAGCg                                                                               >  2:907928/1‑143 (MQ=255)                                                          gCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCAGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGGCCACGCCCAGTAATGCAACCAGCggc                                                                             >  2:451447/1‑141 (MQ=255)                                                               tGCAGCACCTCTGCGCCCGGGTTAGCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTGTGCCGAGCGAAAAGCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATc                                                                        <  2:315942/143‑1 (MQ=255)                                                                gcagcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATcc                                                                       >  1:970470/1‑143 (MQ=255)                                                                 cagctCCACGTGGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTGAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCATGTCCAGTAATGCAACCAGCGTCACATCCg                                                                      <  2:43886/143‑1 (MQ=255)                                                                 cagcaCCGCTCCGCCGTGTGTTCCCGCACGCCCGGCACGGCCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCg                                                                      <  2:174917/143‑1 (MQ=255)                                                                  agcaCCATTTCGACCAGTTTACACGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGCGCCAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCgg                                                                     <  1:869795/143‑1 (MQ=255)                                                                  agcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAACCGCGTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCgg                                                                     >  1:469330/1‑143 (MQ=255)                                                                  aaccccccTTCGCCCTGCTTACCCGCACCCCCGTCACGACCGGCGACCTGCGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCgg                                                                     <  1:127072/136‑1 (MQ=255)                                                                          gtgtcccTGTTTTCTCGCACGCCCGGCACGACCGGCGAGCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGt                                                             <  2:794351/139‑1 (MQ=255)                                                                            cGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGa                                                           <  2:328789/143‑1 (MQ=255)                                                                            cGCCCTGGTGACGCGCGCTCCCGGCACGCCGGCCGACCCGGGGGGAAGGCAGGGCGAGACGCTTGGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGGCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGa                                                           <  2:309963/143‑1 (MQ=255)                                                                              cccTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGACATAGTGAcc                                                         >  1:688404/1‑143 (MQ=255)                                                                                      aCCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCa                                                 <  1:391535/143‑1 (MQ=255)                                                                                       cccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAg                                                <  2:706706/143‑1 (MQ=255)                                                                                        ccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACCCTTTGCCGCGCGCAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCCACCAGCGTCACATCCGGGACATGGTGGCCTCTCGCCAGc                                               >  1:971583/1‑143 (MQ=255)                                                                                         cGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAAGCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCa                                              <  2:272663/143‑1 (MQ=255)                                                                                          gCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCCGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCag                                             >  1:377661/1‑142 (MQ=255)                                                                                            aCGCCCGGCAGCACCGGCGACCTGGGGGTAAAGCTGGGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCAttt                                           <  1:349281/143‑1 (MQ=255)                                                                                                  ggCACGACCGGCGACCTGGGGGTAAAGCTGAGCGAGACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAc                                     <  2:666253/143‑1 (MQ=255)                                                                                                         ccGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCACCCACATATGTGTACCAAACaa                              >  2:71058/1‑143 (MQ=255)                                                                                                            gCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCGAGGAATGCAACCAGCGCCACATCCGGGAACTTGTGACCTCTCGCCAGCAGTTGCGCaccaaccacagt                           >  2:382852/1‑139 (MQ=255)                                                                                                                   gggTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTGGTGTACCAATCAACATCcgcgcg                    >  1:7692/1‑143 (MQ=255)                                                                                                                        tAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCAGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGCTGACCTTTCGCCAGCATTTGAGTACCAATCAAAATCCGCGcgccgcc               >  2:981085/1‑143 (MQ=255)                                                                                                                           aGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCcgcg            >  2:602408/1‑143 (MQ=255)                                                                                                                            gccgAGCGAAACGCTTGGCCGAGCGAAAATCGGCAGAACACGGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGa           <  2:556809/140‑1 (MQ=255)                                                                                                                             cTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGAt          >  1:458282/1‑143 (MQ=255)                                                                                                                               gAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCACAATCCGCGCGCCGCCGCCATGt        >  1:768613/1‑143 (MQ=255)                                                                                                                                  cGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACt     >  2:969199/1‑143 (MQ=255)                                                                                                                                    aaaCGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTc   >  1:583473/1‑143 (MQ=255)                                                                                                                                     aaCGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACGTCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGCTGTAATTCt  >  1:803060/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                         CAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCT  >  NC_000913/4124959‑4125233

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: ‑