Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,943,893 C→T A186V (GCT→GTT)  csdA → cysteine sulfinate desulfinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,943,8930CT100.0% 15.5 / NA 8A186V (GCT→GTT) csdAcysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (2/6);  total (2/6)

AACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTATGGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGAC  >  NC_000913/2943656‑2944110
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
aaCCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCGGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGTTCATTCAGCCggn                                                                                                                                                                                                              <  1:52990/251‑2 (MQ=255)
                                                         ggCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACGCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGTTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGTGGCAGTGCATTTCCCCGCGGATGTTc                                                                                                                                                       <  1:161168/249‑1 (MQ=255)
                                                                                                                                                 aCTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGTTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTGCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATTCGTTTTCAGGTCACAAACTGTATGGCCCGACAGGTATCGGCGTGCTGTATg                                                                            <  1:247774/236‑1 (MQ=255)
                                                                                                                                                                                             cGAACNNTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGTTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCa                                                                                                                                                                             <  2:14376/95‑1 (MQ=255)
                                                                                                                                                                                                            gTTGCCCGGATCTGGCGCGAGCGATTACCTTTGTTCATTCAGCCGGGATGGTGGTGAAGGATGATGGTGCTCCGGGGGCAGTGCATTTCCCCGAGGATGTTCAGCAACTGGATATTGCTTTCTATGCTTTTTCAGGTCACAAACTGTATGGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGCACTGCTGGAGGGGATGTCGCCCTGGCTGGGCGGCGGCACAATGGTTCACGACGTGCGTTTTGAc  >  1:130063/1‑251 (MQ=255)
                                                                                                                                                                                                               gCCCGGATCTGGCGCGAGCGATTACCTTTGTTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCNNCGGATGNNc                                                                                                                                                       >  2:16199/1‑99 (MQ=255)
                                                                                                                                                                                                               gCCCGGATCTGCAGCGAGCGATTACCTTTGTTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCa                                                                                                                                                   <  1:16199/103‑1 (MQ=255)
                                                                                                                                                                                                                      tCTGGCGCGAGCGATTACCTTTGTTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTATGGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGggc                                   <  1:41871/208‑1 (MQ=255)
                                                                                                                                                                                                                               agcgaNTANCNTTNNTNATTCNGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCngng                                                                                                                                                                          >  2:121038/1‑64 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
AACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTATGGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGAC  >  NC_000913/2943656‑2944110

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