Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,583,454 G→A intergenic (+36/‑29) yrhA → / → insA pseudogene, interrupted by IS1E/IS1 repressor TnpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,583,4540GA100.0% 93.7 / NA 31intergenic (+36/‑29)yrhA/insApseudogene, interrupted by IS1E/IS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

GTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAT  >  NC_000913/3583323‑3583563
                                                                                                                                   |                                                                                                             
gTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTatga                                                                                                        <  1:693422/139‑1 (MQ=255)
                              ttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCt                                                                          >  2:565430/1‑139 (MQ=255)
                                      cccTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTAt                                                                  <  1:565430/139‑1 (MQ=255)
                                          gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                                                                      >  2:602959/1‑131 (MQ=255)
                                          gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                                                                      <  1:602959/131‑1 (MQ=255)
                                                cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTccc                                                        >  1:892808/1‑139 (MQ=255)
                                                cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTccc                                                        >  1:74431/1‑139 (MQ=255)
                                                     aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTg                                                                                        <  1:150489/102‑1 (MQ=255)
                                                     aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTg                                                   >  2:911503/1‑139 (MQ=255)
                                                     aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTg                                                   >  1:576453/1‑139 (MQ=255)
                                                     aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTg                                                   >  1:633935/1‑139 (MQ=255)
                                                     aCGGATAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCAAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTg                                                                                        >  2:150489/1‑102 (MQ=255)
                                                       ggTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGtt                                                 >  1:140948/1‑139 (MQ=255)
                                                                   gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                                                           >  2:95233/1‑117 (MQ=255)
                                                                   gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                                                           <  1:95233/117‑1 (MQ=255)
                                                                                ttCAATATTTACTTATGACATTAAAGGTGATGCTGTCAACTTACTGATTTAATGTatgat                                                                                                       >  2:788974/1‑60 (MQ=25)
                                                                                ttCAATATTTACTTATGACATTAAAGGTGATGCTGTCAACTTACTGATTTAATGTatgat                                                                                                       <  1:788974/60‑1 (MQ=255)
                                                                                  cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATgg                                                                                                     >  1:937006/1‑60 (MQ=255)
                                                                                  cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATgg                                                                                                     <  2:937006/60‑1 (MQ=255)
                                                                                         tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAg               <  2:74431/139‑1 (MQ=255)
                                                                                          aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                                             <  1:706543/76‑1 (MQ=255)
                                                                                          aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                                             >  2:706543/1‑76 (MQ=255)
                                                                                          aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                                         >  2:215066/1‑80 (MQ=255)
                                                                                          aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                                         <  1:215066/80‑1 (MQ=255)
                                                                                            ttATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc            <  2:576453/139‑1 (MQ=255)
                                                                                            ttATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc            <  1:649174/139‑1 (MQ=255)
                                                                                                 aCATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgccg       >  1:483125/1‑139 (MQ=25)
                                                                                                     tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTg                                                                                             <  1:1013125/49‑1 (MQ=255)
                                                                                                     tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTg                                                                                             >  2:1013125/1‑49 (MQ=255)
                                                                                                     tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACa   >  2:887199/1‑139 (MQ=17)
                                                                                                      aaaGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAt  <  1:887199/139‑1 (MQ=14)
                                                                                                                                   |                                                                                                             
GTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAT  >  NC_000913/3583323‑3583563

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