Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,470,558 G→A 100% pseudogene (42/3497 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,470,5580GA97.5% 259.9 / ‑2.2 79pseudogene (42/3497 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (49/28);  total (50/29)
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.77e-01

CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAG  >  NC_000913/1470461‑1470654
                                                                                                 |                                                                                                
cAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaa                                                                                                <  1:337040/100‑1 (MQ=255)
cAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATAGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaa                                                                                                >  2:1141011/1‑100 (MQ=255)
     cacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATac                                                                                           <  2:769229/100‑1 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                          <  1:2004237/100‑1 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                          >  1:1682339/1‑100 (MQ=255)
       ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaa                                                                                         >  1:1303901/1‑100 (MQ=255)
       ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaa                                                                                         >  1:1525351/1‑100 (MQ=255)
        cGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaac                                                                                        <  1:585108/100‑1 (MQ=255)
         gAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAAGTTAGGTTTTACTTTTGTTACTGGTAAAAATacaaca                                                                                       >  1:960059/1‑100 (MQ=255)
           aaGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAgg                                                                                     <  1:1797216/100‑1 (MQ=255)
             gagaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGta                                                                                   >  2:1298183/1‑100 (MQ=255)
                 tAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATaaa                                                                               >  2:581079/1‑100 (MQ=255)
                  aaTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAAt                                                                              >  2:2117087/1‑100 (MQ=255)
                  aaTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAAt                                                                              >  1:2245665/1‑100 (MQ=255)
                       aaaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTgg                                                                         <  2:2453657/100‑1 (MQ=255)
                        aaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGc                                                                        <  2:251079/100‑1 (MQ=255)
                            ggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGa                                                                    >  1:1309068/1‑100 (MQ=255)
                               gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtct                                                                 >  1:1644120/1‑100 (MQ=255)
                               gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtct                                                                 >  1:948580/1‑100 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCg                                                               >  1:513787/1‑100 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCg                                                               <  1:2452353/100‑1 (MQ=255)
                                  gCATTGACAGGTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAGATGCAACAGGTATAAATAGTGGCACGATCTCGt                                                              <  2:1686409/100‑1 (MQ=255)
                                  gCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                              >  1:2199115/1‑100 (MQ=255)
                                  gCATTGACAGATTGAACCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                              >  2:262001/1‑100 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  1:581086/100‑1 (MQ=255)
                                       gACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                         >  1:801059/1‑100 (MQ=255)
                                       gACAGATTAGATCTACACTGTTGAAATTCATAATTTACGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  2:1801040/98‑1 (MQ=255)
                                         cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                       >  1:259590/1‑100 (MQ=255)
                                           gATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                     >  2:2061671/1‑100 (MQ=255)
                                           gATTGAATCTACACTGTTGAAATTAAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaa                                                      >  2:300317/1‑99 (MQ=255)
                                            aTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaaa                                                    >  2:2049573/1‑100 (MQ=255)
                                             ttGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAAt                                                   <  2:269594/100‑1 (MQ=255)
                                               gAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATgg                                                 <  1:2215399/100‑1 (MQ=255)
                                               gAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATgg                                                 <  1:2307077/100‑1 (MQ=255)
                                                aaTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGt                                                >  2:1287396/1‑100 (MQ=255)
                                                 aTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTa                                               <  1:597174/100‑1 (MQ=255)
                                                  tCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaa                                              <  1:106791/100‑1 (MQ=255)
                                                   cTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                             <  1:2235574/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                            >  2:2373466/1‑100 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                            >  1:247857/1‑100 (MQ=255)
                                                     acacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGa                                           >  1:91775/1‑100 (MQ=255)
                                                     acacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAAAGGTAAAGa                                           >  2:1125526/1‑100 (MQ=255)
                                                        cTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATcc                                        >  2:2006087/1‑100 (MQ=255)
                                                          gTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCgg                                      >  2:2256243/1‑100 (MQ=255)
                                                           ttGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGc                                     <  2:1781467/100‑1 (MQ=255)
                                                           ttGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTAGCGCGAACTCGTTATTACAAAATGGTAAAGATCCGGc                                     <  2:376913/100‑1 (MQ=255)
                                                            tGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCa                                    >  1:69616/1‑100 (MQ=255)
                                                             gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                   <  1:873531/100‑1 (MQ=255)
                                                             gAAATTCAGAATTTAGGTTTTGCTATTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTAAGTAGAAAATGGTAAAAAtttgacac                                   >  2:2353480/1‑93 (MQ=255)
                                                               aaTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCg                                 >  1:1674716/1‑100 (MQ=255)
                                                                aTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGt                                >  1:480774/1‑100 (MQ=255)
                                                                 ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                               >  1:750521/1‑100 (MQ=255)
                                                                 ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                               >  1:2127059/1‑100 (MQ=255)
                                                                  tCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtct                              <  1:2061685/100‑1 (MQ=255)
                                                                    aGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                            >  1:1223653/1‑100 (MQ=255)
                                                                    aGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                            >  2:2212321/1‑100 (MQ=255)
                                                                     gAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTccc                           >  1:1828063/1‑100 (MQ=255)
                                                                       aTTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                         >  1:203459/1‑100 (MQ=255)
                                                                            ggTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGttt                    >  1:1142602/1‑100 (MQ=255)
                                                                            ggTTTTGCTTTTGTTACTGGCAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGttt                    >  2:1561447/1‑100 (MQ=255)
                                                                             gTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGtttt                   >  1:482700/1‑100 (MQ=255)
                                                                             gTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGtttt                   >  2:1370670/1‑100 (MQ=255)
                                                                             gTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGtttt                   <  1:2527756/100‑1 (MQ=255)
                                                                                ttGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACt                <  1:1609551/100‑1 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGTCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg               <  2:106634/100‑1 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCCCCCATTGTTTTACTg               >  2:186798/1‑100 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             >  1:1521070/1‑100 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             >  1:1522169/1‑100 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             <  1:423700/100‑1 (MQ=255)
                                                                                    ttttGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct            <  1:262004/100‑1 (MQ=255)
                                                                                    ttttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct            >  1:1981318/1‑100 (MQ=255)
                                                                                    ttttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct            <  2:1748315/100‑1 (MQ=255)
                                                                                      ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac          >  1:2503492/1‑100 (MQ=255)
                                                                                      ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac          <  1:610624/100‑1 (MQ=255)
                                                                                      ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac          <  2:69615/100‑1 (MQ=255)
                                                                                        gTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctacta        <  1:2183824/100‑1 (MQ=255)
                                                                                        gTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctacta        >  1:1005735/1‑100 (MQ=255)
                                                                                            cTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgg    >  1:1225234/1‑100 (MQ=255)
                                                                                            cTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgg    >  1:2065313/1‑100 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag  >  1:599468/1‑100 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag  >  1:545031/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAG  >  NC_000913/1470461‑1470654

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