Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,441 2 bp→AG coding (708‑709/1287 nt) thrC → threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,4410GA100.0% 62.4 / NA 21E236E (GAG→GAAthrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/13);  total (8/13)
*NC_0009134,4420AG100.0% 60.6 / NA 21T237A (ACG→GCG) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/13);  total (8/13)

GGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGC  >  NC_000913/4211‑4684
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                    
ggCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGt                                                                                                                                                                                                                                   >  1:63603/1‑251 (MQ=255)
          gTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGACGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGcc                                                                                                                                                                                                                          >  1:227322/1‑250 (MQ=255)
                 gCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCgg                                                                                                                                                                                                                    <  2:63603/249‑1 (MQ=255)
                  cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGATGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                                                                                                                  >  1:133867/1‑250 (MQ=255)
                                  gTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTg                                                                                                                                                                                                    <  2:252302/248‑1 (MQ=255)
                                               cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTgg                                                                                                                                                                                    <  1:128419/251‑1 (MQ=255)
                                                     cGAAACTGCTGCCATCGACGGCGATTTCGATGCCTGCCAGGCGCAGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGATTAACATCAGCCGTTTGGTGGCGCAGATTTGCTACTACGTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGTAAACTTCGGCGATTTGACGGCGTGTCTGCTGGCGAAg                                                                                                                                                                               <  2:37236/250‑1 (MQ=255)
                                                            gTTGCCAACGACGGCGATTTCGAGGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCgg                                                                                                                                                                       <  1:37719/251‑1 (MQ=255)
                                                                                            ggCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGCTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                                                                                                                            >  1:267315/1‑166 (MQ=255)
                                                                                            ggCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCTCAGGAAGCGTGGAACCAGCTGGTTGTCTCGGTg                                                                                                                                                                                                                            <  2:267315/166‑1 (MQ=255)
                                                                                                                      atgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTc                                                                                                            <  1:70866/251‑1 (MQ=255)
                                                                                                                                     tGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGGTTTTTTTCTGCGGCCCCAGTGAACGATACCGTGCCACGGTTCCTGCACCACgggc                                                                                               >  1:125943/1‑247 (MQ=255)
                                                                                                                                                            atcTCGTTTAACAGCGGCCGTTTGCTGTCGCAGATGTGCTAGTACGTGGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTAGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTGTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAg                                                                        <  2:24033/247‑1 (MQ=255)
                                                                                                                                                                cGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGAc                                                                   <  1:3104/250‑1 (MQ=255)
                                                                                                                                                                         tCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTATACTGGCGAAATCACTCACGCTGCCGTGGAAAAGTTTTATTTGCGGGACCAACGTGAACAATACCGTGCCCCCATTACTGCACCACGGTCCGTGGGCACCCAAACCGACTCAGGCGCCGTTATCCAAc                                                         >  1:183573/1‑251 (MQ=255)
                                                                                                                                                                                                       actTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACt                            <  1:10787/249‑1 (MQ=255)
                                                                                                                                                                                                                    tgggcaccTGCCGCAGTAAGCGCGCAACCAGCTGTTTGTCTCGGGGCCAAGCGGAAACTTCGGCGATTTGACGGCTGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGGCAGCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCTTGTGGaa              <  2:180626/243‑1 (MQ=255)
                                                                                                                                                                                                                         aGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTTCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTTCTTCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCCAAGCGACTCAGGCGACGTTATCCAACGCGCTTGACGTGAGTCAGCCCAACAACTGGCCGCGTGTGGAagaggtg        >  1:140074/1‑248 (MQ=255)
                                                                                                                                                                                                                             gCCGCAGGAAGCGCTCAACCAGCTGGTTGTCTCGGTGCCAAGCTGAAACTTCGGCGATTTGACGGCGGGTCTGCTAGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGTTCAGTGGTCACCCAAAGCTACTCAGGCGACGTTATCCAACGCGGTGGGCGGGAGGCAGCCGAGCAACTGGCCGCGTGTTGAAGAtttgtt      <  2:111494/249‑1 (MQ=255)
                                                                                                                                                                                                                             gCCGCAGGAAGCGCTCAACCAGCTGGTTGGCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGGGAACGATACCGTGCCACGTTTCCGGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGGGGAAGAgttgtt      <  2:227322/249‑1 (MQ=255)
                                                                                                                                                                                                                               cGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGGGAAACGTTTTATTGCTGCGACCAACGTGAACGGAACCGTGCCACGTTTCCTGCCCGACGGTCAGTGGTCCCCCAAAGCGACACAGGCAACGTTATCCCACGCGATGGACGGGAGTTCACCCAACCAATGGCCGCGGGTGGAAAGGTTGTTccgc  >  2:213681/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                    
GGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGC  >  NC_000913/4211‑4684

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