Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,187,862 G→C pseudogene (431/2434 nt) yqiG → putative outer membrane usher protein YqiG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,187,8620GC96.0% 89.3 / ‑3.0 25pseudogene (431/2434 nt)yqiGputative outer membrane usher protein YqiG
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (12/12);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATT  >  NC_000913/3187737‑3187991
                                                                                                                             |                                                                                                                                 
aagaaaaaCAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCtt                                                                                                                     <  1:334404/140‑1 (MQ=255)
             tACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTtat                                                                                                        >  1:294206/1‑140 (MQ=255)
                    tATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCctc                                                                                                 <  1:58946/140‑1 (MQ=255)
                          gATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTa                                                                                           <  2:294206/140‑1 (MQ=255)
                              aaGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGca                                                                                       <  1:334795/140‑1 (MQ=255)
                                  gAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGcaccac                                                                                   >  1:171886/1‑140 (MQ=255)
                                                    atCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGc                                                                 <  2:354827/140‑1 (MQ=255)
                                                         cACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAat                                                            >  1:456590/1‑140 (MQ=255)
                                                           ccTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAatat                                                          >  1:139058/1‑140 (MQ=255)
                                                              cAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAg                                                       <  1:75722/140‑1 (MQ=255)
                                                              cAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAg                                                       >  2:432433/1‑140 (MQ=255)
                                                              cAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAg                                                       >  2:298160/1‑140 (MQ=255)
                                                              cAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAg                                                       >  1:63631/1‑140 (MQ=255)
                                                                        tAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTATTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATggg                                             >  2:428572/1‑140 (MQ=255)
                                                                           aaTACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAg                                          >  1:130878/1‑140 (MQ=255)
                                                                                   gATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGt                                  <  2:44246/140‑1 (MQ=255)
                                                                                   gATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGt                                  <  1:298160/140‑1 (MQ=255)
                                                                                    aTGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGtt                                 <  2:359860/140‑1 (MQ=255)
                                                                                     tGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTa                                >  2:81095/1‑140 (MQ=255)
                                                                                              gACGCCTCCTGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTgggg                       >  2:423119/1‑140 (MQ=255)
                                                                                                   ctcctGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTgg                  <  1:423119/140‑1 (MQ=255)
                                                                                                      ctGAACTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGgcgc               <  2:456590/140‑1 (MQ=255)
                                                                                                           cTTTGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCg          <  1:189837/140‑1 (MQ=255)
                                                                                                             ttGGAATCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTa        >  2:58495/1‑140 (MQ=255)
                                                                                                                   tCATGGTATTCCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGAtt  <  2:63631/140‑1 (MQ=255)
                                                                                                                             |                                                                                                                                 
AAGAAAAACAAATTACAACTTATCTCGATAAAGGGAAAATGCAATTATCCATATCTGCACCTCAGGCATGGTTAAAATACAAAGATGCAAACTGGACGCCTCCTGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATT  >  NC_000913/3187737‑3187991

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