Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,389,186 T→C 100% N51S (AAC→AGC)  aaeX ← DUF1656 domain‑containing protein AaeX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,389,1860TC100.0% 70.6 / NA 26N51S (AAC→AGC) aaeXDUF1656 domain‑containing protein AaeX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/14);  total (12/14)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGTTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCC  >  NC_000913/3389047‑3389307
                                                                                                                                           |                                                                                                                         
ttAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGa                                                                                                                        >  1:492801/1‑143 (MQ=255)
 tAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAACCAAGCAGCAATAGAGCGCGGTGCTGaa                                                                                                                       >  1:188141/1‑143 (MQ=255)
 tAAAAATCGCGATGAAGGCCAGAATGACTAGTACGACCGTGATGGCCGTACGGGAGAATTTTCTTACGAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGaa                                                                                                                       <  2:64310/143‑1 (MQ=255)
 tAAAAAGTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGaa                                                                                                                       <  2:367439/143‑1 (MQ=255)
    aaaTTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGATTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGaacaa                                                                                                                    >  2:141497/1‑143 (MQ=255)
            ataaaGGCCAGAACGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTACAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGAt                                                                                                            <  2:259050/140‑1 (MQ=255)
            aTGAAGGGCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGAt                                                                                                            <  2:258265/143‑1 (MQ=255)
               aaGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAACACCAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGcc                                                                                                         >  1:229178/1‑143 (MQ=255)
                gggcgtgAAGGACCAAGACGAGGGTGGTGGACGTACGGGAGAATTTGCTGATTAGTGTTTTCACGTCAAGCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCa                                                                                                        <  2:196373/137‑1 (MQ=255)
                    cAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGACCGACa                                                                                                    >  1:127454/1‑143 (MQ=255)
                       aaTGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAg                                                                                                 >  1:604836/1‑143 (MQ=255)
                            cTATTACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTa                                                                                            <  1:549241/143‑1 (MQ=255)
                            cTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAGCGAGCAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCTGACAAAGTCGTa                                                                                            <  1:348294/143‑1 (MQ=255)
                            cTAATACGACCGTGATGGCCGGACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACACGAAGCAATTTAGCGCTGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTa                                                                                            <  1:593435/143‑1 (MQ=255)
                                tACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAAAAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGGCCTGGATa                                                                                        >  2:448363/1‑143 (MQ=255)
                                  cGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATGAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATAcc                                                                                      >  1:326114/1‑143 (MQ=255)
                                           gggCCGTACGGGAGAATGTTCGTATTAGTGTTTTCACTTCCACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGcac                                                                             <  1:197895/142‑1 (MQ=255)
                                                gTACGGGAGTATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTa                                                                        <  1:568709/143‑1 (MQ=255)
                                                            tttCTTATTAGTGTTTTCACATTAACCACAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCACTAGAGCGTGGTGCTGAACAACGACGGCTGCCAGAAAAAGTCGTAGATACCGAGTGGCACAAGTACCCGGCGCAcaa                                                            >  2:161088/1‑141 (MQ=255)
                                                                     aGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGaaaa                                                   >  1:592537/1‑143 (MQ=255)
                                                                       tgtTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATc                                                 <  1:620000/143‑1 (MQ=255)
                                                                       cgagctCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATc                                                 <  1:155945/138‑1 (MQ=255)
                                                                                                          atCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTTGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTgg              <  1:1201/143‑1 (MQ=255)
                                                                                                           tCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTggg             >  1:563832/1‑143 (MQ=255)
                                                                                                                   aaCAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAg     >  1:27229/1‑143 (MQ=255)
                                                                                                                      ccgccgcaACAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGccc  <  1:102834/138‑1 (MQ=255)
                                                                                                                                           |                                                                                                                         
TTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGTTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCC  >  NC_000913/3389047‑3389307

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: