Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,995,314 C→T pseudogene (447/447 nt) ygeN ← pseudogene, orgB family, part of T3SS PAI ETT2 remnant

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,995,3140CT100.0% 122.4 / NA 38pseudogene (447/447 nt)ygeNpseudogene, orgB family, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/19);  total (19/19)

AAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTAT  >  NC_000913/2995207‑2995447
                                                                                                            |                                                                                                                                     
aaaataaaatTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                         <  1:973339/139‑1 (MQ=255)
                                    aTTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAattatt                                                                                                    <  2:912015/108‑1 (MQ=255)
                                    aTTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAattatt                                                                                                    >  1:912015/1‑108 (MQ=255)
                                              cATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCaaaa                                                                                                              <  2:426068/88‑1 (MQ=255)
                                              cATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCaaaa                                                                                                              >  1:426068/1‑88 (MQ=255)
                                                       atTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAat                                                                                                        <  1:1225022/85‑1 (MQ=255)
                                                       atTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAat                                                                                                        >  2:1225022/1‑85 (MQ=255)
                                                                        aGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAAc                                                                >  2:344000/1‑108 (MQ=255)
                                                                        aGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAAc                                                                <  1:344000/108‑1 (MQ=255)
                                                                        aGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCg                                 >  1:307850/1‑139 (MQ=255)
                                                                             atatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                            <  2:307850/139‑1 (MQ=255)
                                                                                ttCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                         <  1:209071/59‑1 (MQ=255)
                                                                                ttCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                         >  2:209071/1‑59 (MQ=255)
                                                                                   aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTc                      >  2:111542/1‑139 (MQ=255)
                                                                                    aaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttt                                                                      >  1:750515/1‑90 (MQ=255)
                                                                                    aaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttt                                                                      <  2:750515/90‑1 (MQ=255)
                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACa                                                                                                            >  1:679177/1‑48 (MQ=255)
                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACa                                                                                                            <  2:679177/48‑1 (MQ=255)
                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCata                                                                                                >  1:1496560/1‑60 (MQ=255)
                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCata                                                                                                <  2:1496560/60‑1 (MQ=255)
                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCa                                                                                            <  2:966460/64‑1 (MQ=255)
                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCa                                                                                            >  1:966460/1‑64 (MQ=255)
                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                         >  2:1181329/1‑49 (MQ=255)
                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                         <  1:1181329/49‑1 (MQ=255)
                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAattatt                                                                                                    >  1:35728/1‑54 (MQ=255)
                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAattatt                                                                                                    <  2:35728/54‑1 (MQ=255)
                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGtt                               >  1:627847/1‑123 (MQ=255)
                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGtt                               <  2:627847/123‑1 (MQ=255)
                                                                                                 gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCa                                                                              <  2:1555317/69‑1 (MQ=255)
                                                                                                 gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCa                                                                              >  1:1555317/1‑69 (MQ=255)
                                                                                                  aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCa                                                                                                  >  2:1234686/1‑48 (MQ=255)
                                                                                                  aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCa                                                                                                  <  1:1234686/48‑1 (MQ=255)
                                                                                                  aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                            >  2:659807/1‑118 (MQ=255)
                                                                                                  aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                            <  1:659807/118‑1 (MQ=255)
                                                                                                       gggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGattat  >  2:952197/1‑139 (MQ=255)
                                                                                                       gggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGattat  <  1:952197/139‑1 (MQ=255)
                                                                                                        ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGt                                <  1:249144/108‑1 (MQ=255)
                                                                                                        ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGt                                >  2:249144/1‑108 (MQ=255)
                                                                                                            |                                                                                                                                     
AAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTAT  >  NC_000913/2995207‑2995447

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