Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,187,697 (A)6→7 pseudogene (266/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,6911.A95.2% 67.9 / ‑3.2 21pseudogene (260/2434 nt)yqiGputative outer membrane usher protein YqiG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (10/10);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATC  >  NC_000913/3187612‑3187760
                                                                                |                                                                      
ggTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTg                                                                         <  2:531425/80‑1 (MQ=255)
       aTTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGaaaaaaac                                                                <  2:531426/81‑3 (MQ=255)
          gATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGaaaaaaactc                                                              <  2:531428/80‑5 (MQ=255)
          gATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGaaaaaaactc                                                              <  2:531427/80‑5 (MQ=255)
                   aCTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCt                                                    >  2:531453/1‑81 (MQ=255)
                      gCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGAGGATTTTTTTGAAAAAACCTCACCAAACCCTGAt                                                 >  1:531457/1‑81 (MQ=255)
                        aGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCg                                               <  1:531430/81‑1 (MQ=255)
                           ttGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATg                                            >  1:531463/1‑81 (MQ=255)
                                 gACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGgttgt                                       >  1:531462/1‑80 (MQ=255)
                                  acacCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGgttgtt                                      <  1:531431/80‑1 (MQ=255)
                                    acCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTAcc                                   >  1:531449/1‑81 (MQ=255)
                                      cAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCt                                 <  2:531431/81‑1 (MQ=255)
                                         agCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATc                              >  2:531454/1‑81 (MQ=255)
                                          gCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATc                              >  2:531468/1‑80 (MQ=255)
                                            aaTTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCgaa                           <  1:531434/81‑1 (MQ=255)
                                            aaTTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCga                            >  2:531457/1‑80 (MQ=255)
                                            aaTTAACATTATTAGGATTTACTGATGAATTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCgaaa                          <  2:531432/81‑1 (MQ=255)
                                             aTTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCgaaa                          <  2:531434/81‑1 (MQ=255)
                                                 aCATTATTAGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCgaaaaag                       <  2:531435/80‑1 (MQ=255)
                                                    ttattaGGATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCgaaaaagaaaa                   <  1:531437/81‑1 (MQ=255)
                                                          ggATTTACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAAt             <  1:531438/81‑1 (MQ=255)
                                                               tACTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACaa        >  1:531456/1‑81 (MQ=255)
                                                                 cTGATGAATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACt      >  2:531460/1‑81 (MQ=255)
                                                                    atgaATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTAt   >  1:531459/1‑81 (MQ=255)
                                                                     tgaATTTATTGAAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTAtc  <  1:531439/81‑1 (MQ=255)
                                                                                |                                                                      
GGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATGGTTGTTACCCTATCGAAAAAGAAAAACAAATTACAACTTATC  >  NC_000913/3187612‑3187760

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