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% 159.8 / NA 56N51S (AAC→AGC) aaeXDUF1656 domain‑containing protein AaeX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/29);  total (27/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGTTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAACACCACGATAAC  >  NC_000913/3389047‑3389322
                                                                                                                                           |                                                                                                                                        
ttAAAAATTGCGATGAAGTGCAGAATGACTAATACGACCTTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGa                                                                                                                                       <  2:12289/143‑1 (MQ=255)
 tAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGaa                                                                                                                                      >  1:646873/1‑143 (MQ=255)
    aaaTTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTGTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGaacaa                                                                                                                                   >  2:441520/1‑143 (MQ=255)
         gccaTGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCg                                                                                                                              <  2:298568/140‑1 (MQ=255)
           gATGAAGGCCTGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCCGCACTCGAGCGCGGTGCTGCACTCCGCCggt                                                                                                                            >  1:199069/1‑142 (MQ=255)
               aaGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAAAACAACAACCAGGAATGGGGTGCGGTGCTGAACAACGCCGGATGcc                                                                                                                        >  1:411237/1‑143 (MQ=255)
                 ggCCAGAATGACTACTGCGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCCACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAg                                                                                                                      <  1:322448/143‑1 (MQ=255)
                   ccAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGAc                                                                                                                    <  1:642989/143‑1 (MQ=255)
                         tGACTAATGCGACCTTGATGGCCGTACGGGAGCATTTTCTTATTAGTGTTTTCACTTCAACCTCACACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTc                                                                                                              <  2:510677/143‑1 (MQ=255)
                         tGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTc                                                                                                              >  1:585832/1‑143 (MQ=255)
                               aTACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGAt                                                                                                        >  1:141914/1‑143 (MQ=255)
                                  cGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAACAGAGCGCGGTGCTGAACAACGCAGGATGCCAGACAAAGTCGAAGAAAcc                                                                                                     >  2:249224/1‑143 (MQ=255)
                                      cGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTATCACGTCAAGCTCAAGCGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCAGATGCCAGACAAAGTCGTAGATACCTGtt                                                                                                 <  2:850349/143‑1 (MQ=255)
                                      cGTGATGGCCGGACGGGAGGAGTTTGTAATGAGTGCGTGCTCCTCTACCGCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGtt                                                                                                 <  1:133809/143‑1 (MQ=255)
                                       gTGATGGCCGTGCGCGAGAATTTGCGTAGGAGTCTTTTCACTTCAACCCCGCGCCAACAGTCCCGATATCAAATAAAACAAGCAGCAATACAGCCCCGTGCTGAACAACGCCGCATGCCAGACAAAGTCGTAGATACCTGTTg                                                                                                <  2:760210/143‑1 (MQ=255)
                                         gATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGc                                                                                              <  1:542498/143‑1 (MQ=255)
                                          aTGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATCAAACAAGCAGCAATAGCGCGCGGTGCTGAACAACGCCGGAGCCCAGACAAAGTCGTAGATACCAGTTGGca                                                                                             >  2:499241/1‑143 (MQ=255)
                                            ggCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATCTCAAATAAAACAAGCCGCAATAGAACGCGGTGCCGCACAACGCCGGATACCAGCAAAAGTCGTAGACAACGGTTGGccca                                                                                           >  2:158176/1‑143 (MQ=255)
                                            ggCCGTACGGGAGAATTTTCTTATTAGTGTTAACTCTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCAGAACAACGCCGAAGGCGCGAAAAAGTCATAGATGCGTGATGGcaca                                                                                           >  2:72028/1‑143 (MQ=255)
                                               cGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGAAAAAGCACTAGATACCTGTTGGcagacgt                                                                                        >  2:265941/1‑138 (MQ=255)
                                                  aCGGGTGAATTTTCTTATTATTGTTTTCACTTCAGCCTAGAACGAACAGGCGCGGGATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTAcc                                                                                     <  1:66115/143‑1 (MQ=255)
                                                     ggAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCgg                                                                                  >  1:891789/1‑143 (MQ=255)
                                                        gaATTTTCTTATTAGTGTTCTCACTTCAACCTCAAACGAGCAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGgcgc                                                                               <  2:865950/143‑1 (MQ=255)
                                                           ttatCTTAATACTGGTTTCACTTCAACCTCAAACGAACATTCGCGATATCAACTAATACAAGCAGCAATAGAGCGCGCTACTGAACAACGCAGGATGCCAGACACAGTCATCGATGCCTGTTGGCACACGTAGCCGGCGCAcc                                                                            >  2:642386/4‑143 (MQ=255)
                                                            tttCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCCGCAATAGAGCGCGGTGCTGAACGACGCCGGCTGCCTGACAAAGTCGTAGATACCTGTTGGCACAGCGAGCCGGCGCAcca                                                                           >  1:268788/1‑143 (MQ=255)
                                                            tttCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTAGCAAAAGAAACCGACGCAcca                                                                           >  1:840839/1‑143 (MQ=255)
                                                                   tgattgtTTTCACTTCAACCTGAACCGAACATTCGCGATATCAACTAAAACAAGCAGCAATAGAGCGAGGTACTGAACAACGCCGGAGGCTATACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGaa                                                                    <  2:349393/139‑1 (MQ=255)
                                                                       tgtTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATc                                                                <  1:688763/143‑1 (MQ=255)
                                                                         ttttCACGTCCACCTCAGACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGc                                                              <  2:265037/143‑1 (MQ=255)
                                                                         ctttAACTCCAACCCCAAACGAACAGTCGCGATATCAAATCAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTCCCCGGCGCACCAGCCAGAAAATCGc                                                              <  1:969662/142‑1 (MQ=255)
                                                                               cgtCAACCTCACACGAACAGCCGCGACATCACATAATACAGTGGGCTGTAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGa                                                        <  1:929502/141‑1 (MQ=255)
                                                                                tccAACCCCATACGAACCGCCGCGATATCAAATAAAACAAGCAGCAGTAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACCAAGTCGTATATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGAt                                                       <  2:581094/141‑1 (MQ=255)
                                                                                 tCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCCGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTCCCCGGCGCACCAGCCAGAAGATCGCCACTGCTa                                                      >  2:138179/1‑143 (MQ=255)
                                                                                   aaCCTCAAACTAACATTCGCGATATCAAATAAAACAAGGAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATaaa                                                    <  1:501489/143‑1 (MQ=255)
                                                                                    aCCTCAATGGCACACGCGAGATTTCAAATAAACCAAGCAGCAACTGAGCGGGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGGAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATaaaa                                                   <  2:325256/143‑1 (MQ=255)
                                                                                    aCCTCAAACGGCCAGGCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATaaaa                                                   <  2:643793/143‑1 (MQ=255)
                                                                                    aCCTCAAACGAACAGGCGCAACATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCTGGCGCACCAGCCAGAAAATCGCCAGTGATaaaa                                                   <  2:719098/143‑1 (MQ=255)
                                                                                      cTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGAATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGACGCACCAGCCCGAAAATCGCCAGTGAGAAAAGc                                                 >  2:651013/1‑143 (MQ=255)
                                                                                       tCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAGGCa                                                >  1:319962/1‑143 (MQ=255)
                                                                                          aaCGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGCTAAAAGCAAtt                                             >  1:885434/1‑143 (MQ=255)
                                                                                                     gcgATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGCTACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATACAAGCAATTCAAAAACTATc                                  >  2:115225/1‑143 (MQ=255)
                                                                                                        atatCAAATAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGGAAATCGCCAGTGATAAAAGCAATTCAAAAAATAACGGt                               >  2:318072/1‑143 (MQ=255)
                                                                                                          atCAAATAAAGCAGGCAGGAGTAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTgg                             <  2:585832/143‑1 (MQ=255)
                                                                                                             aaaTACGACAAGCAGCAATAGAGCGCGGTGCGGCACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGAACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGaa                          <  2:878323/143‑1 (MQ=255)
                                                                                                             aaaTAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCAAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATGAAAGCAATTCAAAAAATATCGGTGGGaa                          >  1:324109/1‑143 (MQ=255)
                                                                                                              aaTAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGTGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAg                         <  1:441520/143‑1 (MQ=255)
                                                                                                               aTAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACCAAGTCGTAGATACCTGTTGGCACAAGAAACAGGCGCACAAACCAGAAAAACGCAAGTGATAAAAGCAATTCAAAAAATCAGGGTGTGAAgg                        >  1:130004/1‑143 (MQ=255)
                                                                                                               aTAAAACAAGCAGCAATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAACTCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAgg                        >  1:855413/1‑143 (MQ=255)
                                                                                                                  aaccacgcagcaATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGAAATCTGTTGGCAGAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGGGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACa                     <  2:914772/137‑1 (MQ=255)
                                                                                                                         cagcaATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGCAAATCGCCCGTGATAAAAGCAATCCAAAAAATACCAAAGACAAGGACAACCCaaa              >  1:408370/1‑143 (MQ=255)
                                                                                                                         cagcaATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTGGTCACAAGAACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCaaa              <  1:225748/143‑1 (MQ=255)
                                                                                                                           gcaATAGAGCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCcacaca            >  2:511231/1‑143 (MQ=255)
                                                                                                                              aTAGAGTGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAAcacca         <  2:191465/143‑1 (MQ=255)
                                                                                                                              aTAGAGCGCGGTGCTGAACAACGCCGGATGCCGGGCAAAGACGGAGAGACCTGTTGGCACAAGTACCCGGAGCACCAGCCAGAAAATCGCCGGTGAGAAAAGCAATGCAAAAACTATCGATGGGAAGGCGAGCCCAAacaaga         >  2:515561/1‑140 (MQ=255)
                                                                                                                                cgagCGCGGTGCTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAACACCACg       <  1:702974/142‑1 (MQ=255)
                                                                                                                                     gcgGTGCTGAACAACGCCGGATGCCAGACAAAGGCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAACACCACGATAAc  <  1:650995/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
TTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACAAGCAGCAATAGAGCGCGGTGTTGAACAACGCCGGATGCCAGACAAAGTCGTAGATACCTGTTGGCACAAGTACCCGGCGCACCAGCCAGAAAATCGCCAGTGATAAAAGCAATTCAAAAAATATCGGTGGGAAGGACAGCCCAAACACCACGATAAC  >  NC_000913/3389047‑3389322

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