Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,466,533 T→C pseudogene (1142/2518 nt) ydbA → putative outer membrane protein, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,5330TC94.6% 107.1 / ‑1.6 37pseudogene (1142/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/16);  total (20/17)
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.52e-01

GATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAAC  >  NC_000913/1466397‑1466666
                                                                                                                                         |                                                                                                                                     
gATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCAc                                                                                                                                    <  2:590463/140‑1 (MQ=255)
          gTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAg                                                                                                                          >  1:246617/1‑140 (MQ=255)
                  aTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCg                                                                                                                             <  2:563793/129‑1 (MQ=255)
                  aTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCg                                                                                                                             >  1:563793/1‑129 (MQ=255)
                     aaCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTAt                                                                                                               >  1:344814/1‑140 (MQ=255)
                          taataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAgg                                                                                                          <  2:478672/140‑1 (MQ=255)
                            ataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGt                                                                                                                 >  2:931647/1‑131 (MQ=255)
                            ataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGt                                                                                                                 <  1:931647/131‑1 (MQ=37)
                               aCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGAc                                                                                                     <  2:642676/140‑1 (MQ=255)
                                         tATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAgg                                                                                           <  2:262596/140‑1 (MQ=18)
                                           ttCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg                                                                                         >  1:844781/1‑140 (MQ=18)
                                            tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCgg                                                                                        <  1:313313/140‑1 (MQ=21)
                                            tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCgg                                                                                        <  2:246617/140‑1 (MQ=21)
                                             cAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGctgt                                                                                       <  2:344814/140‑4 (MQ=18)
                                                gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                                    >  1:923228/1‑140 (MQ=37)
                                                   ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaa                                                                                 <  1:381807/140‑1 (MQ=255)
                                                     cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                                               >  1:715160/1‑140 (MQ=255)
                                                          tGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGaa                                                                          >  2:716215/1‑140 (MQ=255)
                                                              tGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                             <  1:721763/133‑1 (MQ=255)
                                                              tGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                             >  2:721763/1‑133 (MQ=255)
                                                                    cggcggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                                                <  2:445337/140‑1 (MQ=255)
                                                                      gcggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGcc                                                              <  2:923228/140‑1 (MQ=255)
                                                                                    aTTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTgg                                                <  1:716215/140‑1 (MQ=255)
                                                                                              ccGGCGACAGCGCAACGGGGGATAGCAAGGAGCCTCAGACGGCTATAAAACAGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg                                      <  2:670739/140‑1 (MQ=18)
                                                                                               cGGCGACAGCGCAACGGTGGATAACAATGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGc                                     >  2:512199/1‑140 (MQ=255)
                                                                                                          cAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt                          >  2:116973/1‑140 (MQ=255)
                                                                                                                  gATAACAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGc                  >  2:940218/1‑140 (MQ=255)
                                                                                                                     aaCAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAc               >  2:883480/1‑140 (MQ=255)
                                                                                                                      aCAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACg              <  1:116973/140‑1 (MQ=255)
                                                                                                                          gggCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCGGATTAACGGTGATGACGCCACGGCaa          >  2:3310/1‑140 (MQ=255)
                                                                                                                          gggCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                >  2:601390/1‑134 (MQ=255)
                                                                                                                          gggCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                <  1:601390/134‑1 (MQ=255)
                                                                                                                             gCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac        >  2:733043/1‑140 (MQ=255)
                                                                                                                             gCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac        >  1:73373/1‑140 (MQ=255)
                                                                                                                              cACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaaca       >  1:651708/1‑140 (MQ=255)
                                                                                                                                 cATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaaca    <  2:844781/140‑1 (MQ=255)
                                                                                                                                   tGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaac  >  1:559281/1‑140 (MQ=255)
                                                                                                                                   tGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaac  >  1:50020/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAAC  >  NC_000913/1466397‑1466666

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