Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,425,696 C→A N59K (AAC→AAA priB → primosomal replication protein N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,425,6960CA100.0% 91.0 / NA 26N59K (AAC→AAApriBprimosomal replication protein N
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/15);  total (11/15)

ATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGAC  >  NC_000913/4425461‑4425904
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                
atTTCGCAAACGAAACCGCTGATGATGCTGCAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTcac                                                                                                                                                                                                   >  1:184504/1‑251 (MQ=255)
     gCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCGCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATTCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGtat                                                                                                                                                                                              <  1:12062/251‑1 (MQ=255)
                    gatgatGCTTAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAATGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCGTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTc                                                                                                                                                                                <  2:168980/250‑1 (MQ=255)
                        atgCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATa                                                                                                                                                                           <  2:4061/251‑1 (MQ=255)
                                gCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTc                                                                                                                                                                    >  1:53240/1‑250 (MQ=255)
                                                          gatgaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCa                                                                                                                                          <  1:104388/250‑1 (MQ=255)
                                                                cAACCGTCTGGTGTTGTCCGGCACCGTGTTCAGGTCTCCCCTTCGAAAGTTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTATCGGAAACGAAAAACAGGCCATTTCTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACg                                                                                                                                   <  2:218320/251‑1 (MQ=255)
                                                                        tGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGc                                                                                                                           >  1:21798/1‑251 (MQ=255)
                                                                               gtCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATgg                                                                                                                    >  2:95502/1‑251 (MQ=255)
                                                                               gtCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATgg                                                                                                                    >  1:215267/1‑251 (MQ=255)
                                                                                 ccGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCGGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCTTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtt                                                                                                                  <  1:130030/251‑1 (MQ=255)
                                                                                 ccGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGACCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCACGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTAATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtt                                                                                                                  >  2:193658/1‑251 (MQ=255)
                                                                                  cGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTTAGCAAAATGGttt                                                                                                                 >  2:169700/1‑251 (MQ=255)
                                                                                   ggCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtttt                                                                                                                <  1:94129/251‑1 (MQ=255)
                                                                                           gtgCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGcc                                                                                                         <  2:81586/250‑1 (MQ=255)
                                                                                             gCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAg                                                                                                      <  2:57087/251‑1 (MQ=255)
                                                                                                gggCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGc                                                                                                     >  2:86845/1‑249 (MQ=255)
                                                                                                 ggCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGGTCATTTCATTCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGc                                                                                                     >  2:4746/1‑248 (MQ=255)
                                                                                                         ttCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGattgaattg                                                                                          >  2:183985/1‑251 (MQ=255)
                                                                                                                   cAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTg                                                                                <  1:47999/251‑1 (MQ=255)
                                                                                                                   cAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTg                                                                                <  1:47996/251‑1 (MQ=255)
                                                                                                                       ccATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTGGCATAACCGTTCAGGGGTTGATTTCATGCCGCAAGGCAAAGAGCGGACTGATCAAAATGGTTTTGCATTCCGAGCAGATTGAATTGATAGATTCTGgaga                                                                            <  1:15401/251‑1 (MQ=255)
                                                                                                                            aGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGAGGGAATTGATAGATTCTGGAGACTAGc                                                                       <  2:215267/251‑1 (MQ=255)
                                                                                                                                            cAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTAtt                                                       >  2:100461/1‑251 (MQ=255)
                                                                                                                                                                              aaGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                     <  1:95502/251‑1 (MQ=255)
                                                                                                                                                                                                 ggCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGAc  <  1:173202/251‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                
ATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGAC  >  NC_000913/4425461‑4425904

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: