Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,288,554 G→C pseudogene (361/2525 nt) yejO ← pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,288,5540GC94.2% 233.9 / 7.4 69pseudogene (361/2525 nt)yejOpseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (35/30);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688
                                                                                                                                       |                                                                                                                                      
gaCGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCcccc                                                                                                                                    >  2:470611/1‑139 (MQ=255)
   ggTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGcccccccc                                                                                                                                 <  1:470611/139‑1 (MQ=255)
               aTACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAAc                                                                                                                     <  2:111034/139‑1 (MQ=255)
                         ggCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  1:1054813/139‑1 (MQ=255)
                         ggCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  2:1479661/1‑139 (MQ=255)
                               cgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCt                                                                                                     >  1:1491237/1‑139 (MQ=255)
                                gTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                                    >  1:453344/1‑139 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGcccccc                                                                                                                                   <  1:1320280/91‑1 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGcccccc                                                                                                                                   >  2:1320280/1‑91 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  2:1261472/116‑1 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTGAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  1:1261472/1‑116 (MQ=255)
                                                  ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                                                                   >  2:603294/1‑122 (MQ=255)
                                                  ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                                                                   <  1:603294/122‑1 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:1023210/1‑139 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  1:572998/1‑139 (MQ=255)
                                                           ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         <  2:453344/139‑1 (MQ=255)
                                                             cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                                                                       <  2:572998/139‑1 (MQ=255)
                                                                aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaaaa                                                                     >  2:1611730/1‑138 (MQ=255)
                                                                   tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGt                                                                 >  2:1336685/1‑139 (MQ=255)
                                                                     ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:1479769/139‑1 (MQ=255)
                                                                     ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  1:1553094/1‑139 (MQ=255)
                                                                      cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGc                                                              >  2:325994/1‑139 (MQ=255)
                                                                        caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGGCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcac                                                            >  1:1157952/1‑139 (MQ=255)
                                                                          aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                                                          <  1:325994/139‑1 (MQ=255)
                                                                                    cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACAcc                                                         <  1:1738139/130‑1 (MQ=255)
                                                                                    cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACAcc                                                         >  2:1738139/1‑130 (MQ=255)
                                                                                       ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGt                                             >  1:1295588/1‑139 (MQ=255)
                                                                                       ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGt                                             >  1:222805/1‑139 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:457764/108‑1 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  1:457764/1‑108 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        >  2:422767/1‑115 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        <  1:422767/115‑1 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGcc                                >  1:179154/1‑139 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAg                                                                                                              >  2:1168435/1‑60 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAg                                                                                                              <  1:1168435/60‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  2:69557/1‑107 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  1:69557/107‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcg                               >  1:360018/1‑139 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 <  2:226025/72‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 >  1:226025/1‑72 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         <  1:1331744/60‑1 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         >  2:1331744/1‑60 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          <  2:1724885/139‑1 (MQ=255)
                                                                                                          ccctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              <  2:24120/70‑1 (MQ=255)
                                                                                                          ccctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              >  1:24120/2‑71 (MQ=255)
                                                                                                           cctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAc                         >  2:133578/1‑139 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  2:299515/56‑1 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  1:299515/1‑56 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         >  2:288841/1‑90 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         <  1:288841/90‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTACGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:862623/1‑82 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  2:1782114/54‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  1:1782114/1‑54 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  2:202366/55‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  1:202366/1‑55 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:862623/82‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCa                      <  1:133578/139‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        >  1:132235/1‑104 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        <  2:132235/104‑1 (MQ=255)
                                                                                                               ccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCATCGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:1179723/81‑1 (MQ=255)
                                                                                                               ccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:1179723/1‑81 (MQ=255)
                                                                                                               ccTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCAt                     <  2:1553094/139‑1 (MQ=255)
                                                                                                                cTTTTCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCAtt                    >  2:1019635/1‑139 (MQ=255)
                                                                                                                  tttCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGGGCGTAGTCTGGTTTGCCGCGCCACCCATTTg                  <  1:1336685/139‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:382165/1‑76 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:382165/76‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGCCCCCCCCAGCTTTTTCAACGATTAAGcgccc                                                                                                         <  1:534712/50‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGCACCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         >  2:534712/1‑50 (MQ=255)
                                                                                                                     cGATAGCGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAAc               <  1:166198/139‑1 (MQ=255)
                                                                                                                           cGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                >  1:1045363/1‑52 (MQ=255)
                                                                                                                           cGACAAAGACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                <  2:1045363/52‑1 (MQ=255)
                                                                                                                                  gACGCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:1429798/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688

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