Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,251,287 T→C pseudogene (327/1323 nt) yjbI → pseudogene, SopA‑related, pentapeptide repeats‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,251,2870TC96.8% 209.6 / ‑2.1 62pseudogene (327/1323 nt)yjbIpseudogene, SopA‑related, pentapeptide repeats‑containing
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (34/26);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAA  >  NC_000913/4251192‑4251381
                                                                                               |                                                                                              
gtgATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTc                                                                                            >  2:1279577/1‑100 (MQ=255)
gtgATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTc                                                                                            >  1:212056/1‑100 (MQ=255)
 tgATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTcg                                                                                           <  1:57749/100‑1 (MQ=255)
  gATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTcgc                                                                                          >  2:1270693/1‑100 (MQ=255)
   aTTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTcgcg                                                                                         <  2:1183186/100‑1 (MQ=255)
   aTTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTcgcg                                                                                         <  2:888325/100‑1 (MQ=255)
     ttttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGaa                                                                                       >  1:1266921/1‑100 (MQ=255)
      tttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAAt                                                                                      >  1:1170238/1‑100 (MQ=255)
      tttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAAt                                                                                      <  1:1274040/100‑1 (MQ=255)
       ttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATg                                                                                     <  2:1291163/100‑1 (MQ=255)
        tATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGt                                                                                    >  2:1104009/1‑100 (MQ=255)
         aTCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTc                                                                                   >  1:569098/1‑100 (MQ=255)
          tCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCa                                                                                  >  2:72280/1‑100 (MQ=255)
          tCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCa                                                                                  >  2:414889/1‑100 (MQ=255)
           cAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAg                                                                                 <  1:998396/100‑1 (MQ=255)
           cAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAg                                                                                 <  1:749101/100‑1 (MQ=255)
           cAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAg                                                                                 <  2:1128638/100‑1 (MQ=255)
                 gAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGa                                                                           <  1:465929/100‑1 (MQ=255)
                  aGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGaa                                                                          <  1:972191/100‑1 (MQ=255)
                   gtggtgAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGaaa                                                                         <  1:901632/100‑1 (MQ=255)
                     ggtgAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAAtt                                                                       >  2:768257/1‑100 (MQ=255)
                      gtgAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTc                                                                      >  1:208360/1‑100 (MQ=255)
                        gAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCaa                                                                    >  1:446736/1‑100 (MQ=255)
                             tgtCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTc                                                               <  1:152424/100‑1 (MQ=255)
                              gtCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCg                                                              >  1:621630/1‑100 (MQ=255)
                              gtCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCg                                                              >  1:451099/1‑100 (MQ=255)
                                 ctgttCAGCAAGAATATTTTTTATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTc                                                               <  2:1170599/93‑1 (MQ=255)
                                   tttCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAAt                                                         <  1:1104019/100‑1 (MQ=255)
                                   tttCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAAt                                                         >  1:1063668/1‑100 (MQ=255)
                                      taGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGt                                                      <  2:719587/99‑1 (MQ=255)
                                      cAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGt                                                      >  2:385543/1‑100 (MQ=255)
                                       aGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCACAGTAACTTTCGCGAATGTCGGATTAGAAATTCAAACTTCGATAATTGta                                                     >  1:719592/1‑99 (MQ=255)
                                        gCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTa                                                    <  1:1197924/100‑1 (MQ=255)
                                        gCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTa                                                    <  1:446640/100‑1 (MQ=255)
                                         cAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTAt                                                   >  2:488313/1‑100 (MQ=255)
                                          aaaaaTATTCTTGATACATTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTa                                                    <  2:1186575/98‑1 (MQ=255)
                                             aaTATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCttt                                               >  1:520045/1‑100 (MQ=255)
                                              atatTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTa                                              <  1:1279588/100‑1 (MQ=255)
                                               tatTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTAc                                             <  2:849154/100‑1 (MQ=255)
                                                     ttaTACAGTTGACTTCTCGCAGAGTACAATTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCACACTTCGATAATTGTTATCTTTACGCTTCg                                       >  2:458444/3‑100 (MQ=255)
                                                     tGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCg                                       <  2:826224/100‑1 (MQ=255)
                                                          cAGTTGACTTCTCTCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACtt                                  >  1:882146/1‑100 (MQ=255)
                                                               gACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACca                             >  2:956636/1‑100 (MQ=255)
                                                                aCTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGCTTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcag                            >  1:1105347/1‑100 (MQ=255)
                                                                aCTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcag                            <  1:768262/100‑1 (MQ=255)
                                                                 cTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcaga                           >  2:573749/1‑100 (MQ=255)
                                                                   tctcGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagc                         >  2:641006/1‑100 (MQ=255)
                                                                   tctcGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagc                         <  2:718241/100‑1 (MQ=255)
                                                                      cGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagcaga                      >  1:1169758/1‑100 (MQ=255)
                                                                       gCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagcagag                     >  1:860462/1‑100 (MQ=255)
                                                                        cAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGt                    >  1:1160028/1‑100 (MQ=255)
                                                                             tCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctg               <  1:385547/100‑1 (MQ=255)
                                                                              cAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgt              >  2:1047035/1‑100 (MQ=255)
                                                                              cAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgt              >  1:236892/1‑100 (MQ=255)
                                                                                aCTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgtct            <  2:791753/100‑1 (MQ=255)
                                                                                 cTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgtctg           <  1:381808/100‑1 (MQ=255)
                                                                                     gTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCaa       >  1:554768/1‑100 (MQ=255)
                                                                                       cATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAg     >  1:401789/1‑100 (MQ=255)
                                                                                       cATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAg     <  1:929254/100‑1 (MQ=255)
                                                                                        aTAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGa    >  1:1173126/1‑100 (MQ=255)
                                                                                         tAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGaa   >  1:777325/1‑100 (MQ=255)
                                                                                          aGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGaaa  >  1:136034/1‑100 (MQ=255)
                                                                                               |                                                                                              
GTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAA  >  NC_000913/4251192‑4251381

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