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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,125,0970CT100.0% 106.6 / ‑6.7 41E572K (GAG→AAG) priAprimosome factor N'
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/16);  minor base G (1/2);  total (23/18)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG  >  NC_000913/4124958‑4125234
                                                                                                                                           |                                                                                                                                         
gCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGc                                                                                                                                        <  2:739372/143‑1 (MQ=255)
gCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGc                                                                                                                                        >  2:311859/1‑143 (MQ=255)
     ggCGTCGTAGCCTTTATAGAGCAACGTATGCAACAGAGGATGTTCCGGCTGGTGCGTTTGCAGCGCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGc                                                                                                                                   <  2:435738/143‑1 (MQ=255)
            tAGCCGTTATAGAGAAACGGTTTCAACAGAAGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCAGTGTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGCAAAGCTGAGCGAAACGCCTTGCCGGGCGAACATc                                                                                                                            >  2:623724/1‑143 (MQ=255)
              gCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCgg                                                                                                                          <  2:756245/143‑1 (MQ=255)
                     gagagCAACGTTGGCAACAGAGGATGTTCCGGATGGTGCGTTTGCCGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAGAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGaaaa                                                                                                                   <  2:66168/142‑1 (MQ=255)
                       gagCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACa                                                                                                                 >  2:321308/1‑143 (MQ=255)
                       gagCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAAGAGGCGGaaacta                                                                                                                 >  2:807259/1‑140 (MQ=255)
                        agCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAg                                                                                                                >  1:92386/1‑143 (MQ=255)
                                    aaCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGCCCTGGGTGTACAGCTGTGCGAAACGCTTTGGCGAGCGGCAATCGGCAGAAAACAGCGCGCcggccac                                                                                                    >  2:700773/1‑143 (MQ=255)
                                      cAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGGACGCCCGGGACGACCGGCGACCCGGGTGTAAGGCTGAGCGAAACGCTGTGGCGAGGGAAAGACGGCAGAGAACAGGGCGCcgtgagcga                                                                                                  >  2:545030/1‑137 (MQ=255)
                                                tCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGACAACAGCGCGCCGTCCACGTCCAGTAATGc                                                                                        >  1:729941/1‑143 (MQ=255)
                                                             tttGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTcaca                                                                           <  1:527938/143‑1 (MQ=255)
                                                                gcagcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCCCGTCCAGTCACGCAACCAGCGTCACATcc                                                                        >  1:807876/1‑143 (MQ=255)
                                                                 cagcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCCGAAAACAGCGCGCCGTCCACGTCCAGTAACGCTACCAGCGGCACATGCg                                                                       >  2:486731/1‑143 (MQ=255)
                                                                 cagcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCg                                                                       <  1:614208/143‑1 (MQ=255)
                                                                   gcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCCGCGTCACATCCggg                                                                     >  1:22312/1‑143 (MQ=255)
                                                                    caccacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCCGAGCGAAACCCCTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAAGGCAACCAGCGTCCCATCCGGGa                                                                    >  2:111391/1‑143 (MQ=255)
                                                                       cacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCAGGTCGCGTCATGCAACCTGCGTCACGTCgggggtat                                                                 >  1:250676/1‑139 (MQ=255)
                                                                       cacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGaaac                                                                 >  1:779547/1‑142 (MQ=255)
                                                                               ccTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGCAATGGCGACCt                                                         >  2:91138/1‑143 (MQ=255)
                                                                                 tggcgccccGCACGCCCGGCACGACCGGCGACCTGGGTGTAACGCTGCGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCTTCACATCCGGGAAATGGTGACCttt                                                       <  1:465500/137‑1 (MQ=255)
                                                                                    ttgcccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGc                                                    <  1:15249/140‑1 (MQ=255)
                                                                                      aCCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCa                                                  <  1:621326/143‑1 (MQ=255)
                                                                                              gCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTtgt                                          <  2:79679/143‑1 (MQ=255)
                                                                                              gCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGCGACCTTTCGCCAGCATAtgg                                          >  1:912677/1‑142 (MQ=255)
                                                                                                  ggCGCGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAc                                      <  2:68677/143‑1 (MQ=255)
                                                                                                  ggCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTGTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAc                                      <  1:654348/143‑1 (MQ=255)
                                                                                                       gACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACAGCCGGGAAATGGTGACCTGTCGCCAGCATTCGTGTCCCAATc                                 >  1:94623/1‑143 (MQ=255)
                                                                                                              gACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGACACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAACAATCAACATcc                          >  1:374829/1‑143 (MQ=255)
                                                                                                               aCCTGCGTGCAAAGGTAAGCTAAACTCTTTGCCGAGCGAAAATCGGCACAAAACAGCGCGCCGTCCACGTCTAGTAATTCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCcg                         <  2:502180/143‑1 (MQ=255)
                                                                                                                ccTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTCATGCAACCAGCGTCACATATCGGAGATGGTGACCTTTCGCCAGCATTTGTGTGCCAATCAGAATCcgc                        >  1:893680/1‑143 (MQ=255)
                                                                                                                ccTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCcgc                        <  2:733218/143‑1 (MQ=255)
                                                                                                                  tGGGGGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCcgcgc                      <  2:92386/143‑1 (MQ=255)
                                                                                                                    ggTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCgcgcgc                    <  1:293254/143‑1 (MQ=255)
                                                                                                                       gtAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGcgccgc                 >  1:806228/1‑143 (MQ=255)
                                                                                                                        tAAAGCTGAGCGAAACGATTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTCGTGTACCAAGCAAAATCCGCGcgccgcc                >  2:337066/1‑143 (MQ=255)
                                                                                                                            gCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGa            <  2:612190/143‑1 (MQ=255)
                                                                                                                            gCTCAGCGAAACGCTGTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGa            <  1:621626/143‑1 (MQ=255)
                                                                                                                                   gAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACGTCCGGGGAATGGTGACGTTTCGCCACCCTTTGTGTACCCGTCGCTCTCCGCGCGCCCCGGCGACGTACtt     >  1:278636/1‑143 (MQ=255)
                                                                                                                                      aCGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATTGTTACCTTTCGCCAGCGTTTGTGTGCCACTCCAAATCCGCGCGCCGCCGCGATGTACTTCTg  >  1:256269/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                         
GCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG  >  NC_000913/4124958‑4125234

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: