Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,638,571 G→A H33Y (CAT→TAT)  cspI ← cold shock‑like protein CspI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,638,5710GA100.0% 71.9 / NA 20H33Y (CAT→TAT) cspIcold shock‑like protein CspI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/13);  total (7/13)

TTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACA  >  NC_000913/1638332‑1638806
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           
ccctggaaatgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACAt                                                                                                                                                                                                                                  >  1:382218‑M2/136‑251 (MQ=255)
    ggaaatgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCttt                                                                                                                                                                                                                              >  1:255053‑M2/132‑251 (MQ=255)
                    acccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                                                                                                                                                              <  2:261519‑M2/136‑1 (MQ=255)
                                          atgtccaccctcaagaatttttattcccggcgtaacggccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACAAGAACAGCGGCAAGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                                                                                                                                        <  2:254358‑M2/158‑1 (MQ=255)
                                                   gtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaattggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGACCGTTCTCAATACCAAATTCCACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACCCATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                                                                                                                                                                >  2:134237‑M2/85‑250 (MQ=255)
                                                                  acccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaaa                                                                                                                                                                 <  1:402903‑M2/181‑1 (MQ=255)
                                                                                 cgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGtt                                                                                                                                                   >  1:31195‑M2/55‑249 (MQ=255)
                                                                                   ggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCGTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAc                                                                                                                                               <  1:411366‑M2/199‑1 (MQ=255)
                                                                                         tgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAAGTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCGTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                                                                                                                                          <  2:179901‑M2/204‑1 (MQ=255)
                                                                                                 acaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACAt                                                                                                                                                                                                                                  >  1:420645‑M2/39‑154 (MQ=255)
                                                                                                 acaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTGAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACAt                                                                                                                                                                                                                                  <  2:420645‑M2/116‑1 (MQ=255)
                                                                                                                      gcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTgg                                                                                                             <  2:149364‑M2/233‑1 (MQ=255)
                                                                                                                                                                          cAATACCAAATTCAACTTCCTGATTCGCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCGTTTGAAGAAATTAGCCCTGGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAAtct                                                          <  2:242913/249‑1 (MQ=255)
                                                                                                                                                                               ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCTAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAgg                                                   <  1:201979/251‑1 (MQ=255)
                                                                                                                                                                                       aaCTTCCTGATTCTCAGTTACTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCa                                                        >  2:433758/1‑238 (MQ=255)
                                                                                                                                                                                       aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAAGTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCa                                                        <  1:433758/238‑1 (MQ=255)
                                                                                                                                                                                               gATTCGCAGTTAATGTTTGGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCATTTGAAGAAATTAGCCCTTTGGCAGAAGGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATaaaa                                     <  2:390114/249‑1 (MQ=255)
                                                                                                                                                                                                     cAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACa                             <  2:255053/251‑1 (MQ=255)
                                                                                                                                                                                                             gTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAAGATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTTCAAATAATAAAATTAGAACAATGACTGc                     <  1:85050/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTTGCTGAGAAATAGACAAACACATCCTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACa  >  1:402327/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           
TTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACA  >  NC_000913/1638332‑1638806

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