Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,336 C→G 7.5% pseudogene (945/2513 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,466,3360CG 7.5% 167.3 / 3.7 52pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base G (2/2);  total (29/23)
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.72e-01

CGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATT  >  NC_000913/1466198‑1466483
                                                                                                                                          |                                                                                                                                                   
cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                           <  1:21050/149‑1 (MQ=255)
 ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATATGGAGTCCATGTGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                                                                                                          <  2:109357/149‑1 (MQ=255)
  gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt                                                                                                                                         <  5:52997/149‑1 (MQ=255)
   tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                        <  2:35335/149‑1 (MQ=255)
     aCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                                                                                                      <  8:225887/149‑1 (MQ=255)
       aGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTACATGGGTATCCAGATCGACGGTGACAAGGCCAACGTCAATAACGAAGGCGAGCGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgattc                                                                                                                                    >  3:334035/1‑147 (MQ=255)
                   gATTACAAGGTCAATATGAACGTAATCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGAAAAGGACATCGTCAATAACGATTGCGAGAGCACCATCACCTACGGTGGCACCGGCACCCAGTTTAATGGTGATGACGCCACGGCaaa                                                                                                                        <  8:148730/149‑1 (MQ=18)
                     tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGGACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                                                                                                      <  6:127146/149‑1 (MQ=255)
                     tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGCTGACGCCACGGCAaac                                                                                                                       >  3:40176/1‑148 (MQ=255)
                     tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                       >  1:311559/1‑148 (MQ=255)
                      aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCCCCGGTGGCACCGGCACCCCGATTAATGGTGATGACGCCCCGGCAaacaa                                                                                                                     >  8:379205/1‑149 (MQ=255)
                          aGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                 <  8:319012/149‑1 (MQ=255)
                             gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCATACAACAACGGc                                                                                                              <  6:140604/149‑1 (MQ=255)
                                 cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGCCAAGGCCATCGTCAATAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                          >  5:56332/1‑149 (MQ=255)
                                       cGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                                                                                                    >  7:333902/1‑149 (MQ=255)
                                             cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                              <  7:401213/149‑1 (MQ=255)
                                             cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGACAAAAACGAAGGCGAGAGCCCCATCACCACCGGTGGCACCGGCACCCACATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                              >  1:165950/1‑149 (MQ=17)
                                               aTCCGGAGTCCATGGGTATCCAGATCGCCGGTGACAAGGCCATCGTCAATAACGAAGGCGCGAGCACAATCACCAACGGTGGCCACGCCACCCACATTAATGGTGATGACGCCCAGCAAAACAACAACGGCAAACACACCGtgcacgt                                                                                             >  2:313285/1‑142 (MQ=18)
                                                  cGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                          <  2:311559/148‑1 (MQ=255)
                                                    gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                  <  6:234845/138‑1 (MQ=255)
                                                    gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                  >  5:234845/1‑138 (MQ=255)
                                                         cATGGGTATCCAGATTGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                                                  >  1:168239/1‑149 (MQ=255)
                                                          aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                 >  7:98381/1‑149 (MQ=255)
                                                              gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAACACCACCGTTGACGGCAAGGCTTCCACCgg                                                                              >  1:384315/1‑148 (MQ=255)
                                                               tATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                    <  6:96723/109‑1 (MQ=255)
                                                               tATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                    >  5:96723/1‑109 (MQ=255)
                                                                  ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACg                                                                          >  3:335157/1‑148 (MQ=255)
                                                                          aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                 >  5:378934/1‑149 (MQ=255)
                                                                                aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCCAACAACAACGGCAACACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaa                                                            >  5:276736/1‑148 (MQ=255)
                                                                                aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaa                                                            >  1:129095/1‑148 (MQ=255)
                                                                                       cATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                                     >  7:203911/1‑148 (MQ=255)
                                                                                         tCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                  <  5:287750/149‑1 (MQ=255)
                                                                                              aaTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTa                                              <  2:129095/148‑1 (MQ=255)
                                                                                              aaTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTa                                              <  2:168239/148‑1 (MQ=255)
                                                                                                tAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                              <  4:34392/98‑1 (MQ=9)
                                                                                                tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                              >  3:34392/1‑98 (MQ=9)
                                                                                                tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCTGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                <  5:262234/96‑1 (MQ=255)
                                                                                                tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCGGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                >  6:262234/1‑96 (MQ=255)
                                                                                                tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                            >  5:294717/1‑148 (MQ=255)
                                                                                                tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                           <  1:11870/149‑1 (MQ=255)
                                                                                                    gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACCACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa                                       >  1:52074/1‑149 (MQ=255)
                                                                                                      aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACg                                     >  7:255849/1‑149 (MQ=255)
                                                                                                       ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                                               >  5:82855/1‑122 (MQ=255)
                                                                                                       ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                                               <  6:82855/122‑1 (MQ=255)
                                                                                                       ggCGAGAGCACCATCACCAACGGTGGCAACGGCACCCAGATTAATGGTGATGACGACACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACgg                                    <  2:366044/149‑1 (MQ=255)
                                                                                                              gCACCATCACCAACGGGGGCACCGGCACCCAGATTACTGGTGATGACGACACGGAAAACAACAACGGCAAAACCAACGTTGCCGGCAAGGATTCCACCGGTACGGTAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                             <  1:313285/149‑1 (MQ=17)
                                                                                                                  cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                         >  8:72528/1‑149 (MQ=255)
                                                                                                                      aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCa                      <  8:98381/148‑1 (MQ=255)
                                                                                                                      aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                     >  7:257996/1‑149 (MQ=255)
                                                                                                                        cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcg                   >  8:142891/1‑149 (MQ=255)
                                                                                                                        cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcg                   >  2:106158/1‑149 (MQ=255)
                                                                                                                         aaCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAAGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcgg                  >  1:150882/1‑149 (MQ=255)
                                                                                                                                   ccggcaccCAGATTTATGTTCATGACGCCACGGAAAACAACAACGGCAAATCCACCTTTGTCGGCAAGGATTCCACCGTTACGGAATTCAACGGTAATAACGGAAAAGTTATTCAGGACGTCGATCTGGATGTCATCGGCGGCGGTCAc        <  4:334035/149‑1 (MQ=255)
                                                                                                                                    cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg       >  7:320309/1‑149 (MQ=255)
                                                                                                                                    cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg       <  4:40176/149‑1 (MQ=255)
                                                                                                                                      gcaccCAGATTAATGGTTATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTTACGCCAAGGATTCCACCTCTACGGAAATCAACGGTAATAAAGGAAAAGGAATTAAGGAATACGAACTGGATGTCAGAGGCCGCGGTCAcag      >  2:303092/1‑146 (MQ=17)
                                                                                                                                        accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt   <  8:333902/149‑1 (MQ=255)
                                                                                                                                         cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGt                                                              <  8:118810/89‑1 (MQ=255)
                                                                                                                                         cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGt                                                              >  7:118810/1‑89 (MQ=255)
                                                                                                                                          ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAAAGGAAATAACGGAAAAGTTATTCAGTACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt  >  4:160322/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
CGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATT  >  NC_000913/1466198‑1466483

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