Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,593,020 G→A pseudogene (5067/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,593,0200GA92.6% 78.8 / NA 28pseudogene (5067/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (9/16);  minor base C (0/2);  total (9/19)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGGAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCAAAGGTATACCCACCCTCAACAGA  >  NC_000913/1592893‑1593135
                                                                                                                               |                                                                                                                   
gAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTTGAAATCt                                                                                                                >  2:93976/1‑133 (MQ=255)
  aCTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCttt                                                                                                              <  2:311763/133‑1 (MQ=255)
            tttACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATc                                                                                                    >  1:49659/1‑133 (MQ=255)
            tttACCATCGTCACGATGGTGATGGCTGGTCAGGTAGGTTTACACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATc                                                                                                    <  2:189255/133‑1 (MQ=255)
                ccATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAgg                                                                                                <  1:447403/133‑1 (MQ=255)
                        aCGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTgc                                                                                        >  1:469732/1‑133 (MQ=255)
                          gCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTgcgc                                                                                      >  2:455603/1‑133 (MQ=255)
                              gTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTACGGTCGAAAACTTTCACACCCATCCAGGTGAGTTGCGCCtgt                                                                                  >  2:163082/1‑133 (MQ=255)
                                gATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTgg                                                                                <  2:109109/133‑1 (MQ=255)
                                gATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTgg                                                                                <  2:126685/133‑1 (MQ=255)
                                        ttCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGt                                                                        >  1:45258/1‑133 (MQ=255)
                                        ttCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTAATTGGCGCCTGTGACTGGACGt                                                                        >  1:153055/1‑133 (MQ=255)
                                           aggtaggtTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTAcc                                                                     <  1:311111/133‑1 (MQ=255)
                                             gtaggtTTTCACCCCAAGTCGCGTTTGTACATTTCCGTCGCCTTCCGTTTCAATGCCCGTTCCGTCTTTCCGGGTATGGCCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAg                                                                   >  2:300770/1‑133 (MQ=255)
                                                ggtTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTa                                                                <  1:455603/133‑1 (MQ=255)
                                                      cACCCCAACTCGCGTTTGGACACTTCCGTCCCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGc                                                          <  1:173589/133‑1 (MQ=255)
                                                       aCCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCg                                                         >  1:35729/1‑133 (MQ=255)
                                                            aaGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCt                                                    >  1:102448/1‑133 (MQ=255)
                                                                 gcgTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCt                                               <  2:469732/133‑1 (MQ=255)
                                                                                cGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTcc                                <  1:76860/133‑1 (MQ=255)
                                                                                cGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTcc                                <  1:300770/133‑1 (MQ=255)
                                                                                      cTTCCGTTTCAATGCCCGTTCCGTTTTTCCGGGTATGGTCGCAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTc                          <  2:134621/133‑1 (MQ=255)
                                                                                      cTTCCGTTTCAATGCCCGTTCCGTCTCTCCGGGTATGGTCGCAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTc                          <  2:223058/133‑1 (MQ=255)
                                                                                       ttCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCa                         <  1:96256/133‑1 (MQ=255)
                                                                                       ttCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCa                         <  2:49659/133‑1 (MQ=255)
                                                                                           gTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCAAAgg                     <  1:93976/133‑1 (MQ=255)
                                                                                                 aTGCGCGTCCCGTCTTTCGGGGAATGGTCGGAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCAAAGGTATacc               <  1:173159/133‑1 (MQ=255)
                                                                                                  tGCGCGTTCCGTCTTTCCGGGTATGGTCGAAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCAAAGGTATaccc              <  1:14261/133‑1 (MQ=255)
                                                                                                              cTTTCCGGTTATGGTGGACATCATTCACACACATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCAAAGGTATACCCACCCTCAACaga  <  1:163082/133‑1 (MQ=255)
                                                                                                                               |                                                                                                                   
GAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGGAATCTTTCACACCCATCCAGGTGATTTGCGCCTGTGGCTGGACGTACCAGGTATTCAGCGTCCCTTCGCTGCCGCTAAATGTTCCCGCTTCAAAGGTATACCCACCCTCAACAGA  >  NC_000913/1592893‑1593135

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