Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,573,528 A→G pseudogene (1115/1503 nt) yjiT → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,573,5280AG97.3% 252.3 / ‑2.2 77pseudogene (1115/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (37/36);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

AAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAA  >  NC_000913/4573432‑4573620
                                                                                                |                                                                                            
aaaccaaaGTTAGATTTGGCAGCCAAAATCCATTGTCAGGTTTGTCGTTAGTGGCTCTTGCTGGAGGGATGATTGTTTGGAGCATTAAACTTGATGGCAg                                                                                           <  1:1643513/94‑1 (MQ=255)
  aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTg                                                                                         <  1:1815187/100‑1 (MQ=255)
    cGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaa                                                                                       <  1:937430/100‑1 (MQ=255)
     gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                      >  1:665928/1‑100 (MQ=255)
     gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                      >  1:1261759/1‑100 (MQ=255)
        gTGAGATTTGGGAGGCTAAAACCATCGGTATGTTTGGCGTTAGTGGCTCGTGCTGGAGGGATGATTATTGGGTGCATTAAACTTGATGCCAGTGaaagt                                                                                    <  1:792854/99‑3 (MQ=255)
        gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                   >  2:482939/1‑100 (MQ=255)
        gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                   >  1:1429322/1‑100 (MQ=255)
          gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                 <  2:1261770/100‑1 (MQ=255)
           agaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                 <  2:677867/99‑1 (MQ=255)
           agaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCa                                                                                <  1:615954/100‑1 (MQ=255)
            gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                               >  1:1267239/1‑100 (MQ=255)
            gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                               >  1:441761/1‑100 (MQ=255)
             aTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                               <  2:441763/99‑1 (MQ=255)
              tttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTg                                                                             <  2:1515234/100‑1 (MQ=255)
               ttGGCAGGCTAAATCCATGGGCAGGTTTGTCGTTCGTGGCTCTTGCTGGAGGGGTGATTGTTGGGAGCATTAAACTTGATGGCAGTGAGATTGca                                                                                 <  1:1790820/95‑2 (MQ=255)
                 ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                          <  1:526162/100‑1 (MQ=255)
                 ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                          <  1:1045029/100‑1 (MQ=255)
                 ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                          >  2:725366/1‑100 (MQ=255)
                   caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAACTTGCCATTGgaga                                                                         >  2:1126867/1‑99 (MQ=255)
                      gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                     <  1:262941/100‑1 (MQ=255)
                       caAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGc                                                                    <  2:1467362/100‑1 (MQ=255)
                         aaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCg                                                                  <  1:962451/100‑1 (MQ=255)
                          aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGt                                                                 >  1:51833/1‑100 (MQ=255)
                           aTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGtt                                                                <  1:727428/100‑1 (MQ=255)
                            tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTa                                                               <  1:1267988/100‑1 (MQ=255)
                              cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAAc                                                             <  2:1626951/100‑1 (MQ=255)
                               aTCGGCAGGTTTGTCGTTAGTGGTTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                            <  2:1055766/100‑1 (MQ=255)
                                  ggCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTc                                                         >  1:1405619/1‑100 (MQ=255)
                                    cAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCAt                                                       <  2:830917/100‑1 (MQ=255)
                                      ggTTTGTCGTTAGTGGCTCGTGCTGGGGGGATGATTGTTGGGAGCATTAAACTTGATTTCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATc                                                      >  2:989368/1‑99 (MQ=255)
                                        tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGtt                                                   <  1:1421306/100‑1 (MQ=255)
                                         ttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTg                                                  <  2:19793/100‑1 (MQ=255)
                                          tGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGa                                                 <  1:1197777/100‑1 (MQ=255)
                                            tCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATg                                               >  1:710121/1‑100 (MQ=255)
                                             cGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                              >  1:1809934/1‑100 (MQ=255)
                                              gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                             >  1:218860/1‑100 (MQ=255)
                                                 aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                          >  1:914311/1‑100 (MQ=255)
                                                 aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                          <  2:710125/100‑1 (MQ=255)
                                                 aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                          <  2:606950/100‑1 (MQ=255)
                                                  gTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATc                                         >  1:604357/1‑100 (MQ=255)
                                                  gTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATc                                         <  1:542702/100‑1 (MQ=255)
                                                    ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCaa                                       >  1:37962/1‑100 (MQ=255)
                                                    ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTAAGGTGCCGTTAACCTTCATCGTTGATGCAGATCaa                                       >  2:626207/1‑100 (MQ=255)
                                                    ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCaa                                       >  1:873012/1‑100 (MQ=255)
                                                     gCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAAt                                      <  2:604361/100‑1 (MQ=255)
                                                      cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAAt                                      >  2:1554122/1‑99 (MQ=255)
                                                      cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                     >  2:1048472/1‑100 (MQ=255)
                                                        cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                   <  1:482984/100‑1 (MQ=255)
                                                         gTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCt                                  <  2:1267597/100‑1 (MQ=255)
                                                          ttcTGGAGTTATTATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTg                                 <  2:715231/98‑1 (MQ=255)
                                                          tGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                     >  2:882195/1‑96 (MQ=255)
                                                          tGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                   <  2:1086519/98‑1 (MQ=255)
                                                           cttGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGt                                <  1:1381459/98‑1 (MQ=255)
                                                            cTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                               <  2:742977/100‑1 (MQ=255)
                                                             tGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtag                              >  2:1381469/1‑98 (MQ=255)
                                                              ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                             >  1:1210189/1‑100 (MQ=255)
                                                              ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                             >  1:515405/1‑100 (MQ=255)
                                                              ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                             >  1:1287312/1‑100 (MQ=255)
                                                               gcgggATGATTGTTGGGAGCATTAAACATGATAGCAGTGAAATTGCGCTTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAAATCAATGGCTGTTACa                            >  1:252812/3‑100 (MQ=255)
                                                               gAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACa                            >  2:548849/1‑100 (MQ=255)
                                                                 gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                          >  1:319613/1‑100 (MQ=255)
                                                                 gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                          >  1:111118/1‑100 (MQ=255)
                                                                 gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                          >  1:1019199/1‑100 (MQ=255)
                                                                  ggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAggg                         <  2:1071654/100‑1 (MQ=255)
                                                                        ttgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGc                   <  1:898119/100‑1 (MQ=255)
                                                                           ttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttc                >  1:1124269/1‑100 (MQ=255)
                                                                            tgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttct               >  1:495180/1‑100 (MQ=255)
                                                                             gggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt              <  2:165357/100‑1 (MQ=255)
                                                                                 gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg          <  1:593827/100‑1 (MQ=255)
                                                                                  cATTAAACTTGATGGCAGTGAAATTGCCAGTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAGCAATGGCTGTTACAGGGACAAGCTTCTTTCa           >  1:1415313/1‑98 (MQ=255)
                                                                                   aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg        >  2:402018/1‑100 (MQ=255)
                                                                                   aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg        >  1:1343962/1‑100 (MQ=255)
                                                                                   aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCAAAAATACAAAACAAAAACTGTTACAAAGACAAGCTTCTTACAgtg        >  2:82286/1‑100 (MQ=255)
                                                                                      aaaCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGc     >  2:1185687/1‑100 (MQ=255)
                                                                                       aaCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGc     <  2:1082145/99‑1 (MQ=255)
                                                                                         cTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGaa  >  1:1644515/1‑100 (MQ=255)
                                                                                                |                                                                                            
AAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAA  >  NC_000913/4573432‑4573620

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