Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,336 C→G 13.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,3360CG13.5% 136.7 / 6.9 56pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base G (5/3);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAT  >  NC_000913/1466189‑1466482
                                                                                                                                                   |                                                                                                                                                  
tGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa                                                                                                                                                  <  2:903653/150‑1 (MQ=255)
 gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCGGGAGTCCATGGGTAACCAGAACGACGGTGACAAGGCCATCGTCAATTACGAAGGCGAGCGCCCCAACACCAACGGTGGCTCCGGCACCCAg                                                                                                                                                 >  2:522561/1‑150 (MQ=17)
 gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                                                                                                 <  2:299320/150‑1 (MQ=255)
 gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATAGACGGAGACAAGGCCATAGTAAAAAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAg                                                                                                                                                 >  2:1126650/1‑150 (MQ=14)
    atCACCGGTGCCAGCGCGACGGTGGATAACAATGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAATGCCATCGTCAATAACGAAGGCGATAGCACCATCTCGAACGGTGGCACCGGCACCCAGAtt                                                                                                                                              <  2:408122/150‑1 (MQ=17)
    atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGCTCGACGGTGACAAGGCCATCGTCAAAAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCAGCtt                                                                                                                                              >  1:1125551/1‑150 (MQ=255)
     tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTa                                                                                                                                             <  1:418304/150‑1 (MQ=255)
     tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTa                                                                                                                                             >  1:410280/1‑150 (MQ=255)
      cACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGGGACAAGGCCATCGTCAAAAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                                                                                                            >  1:1321711/1‑150 (MQ=255)
        ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                          <  2:974144/150‑1 (MQ=255)
         cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                                                                                                         <  1:823963/150‑1 (MQ=255)
         cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                                                                                                         <  1:942788/150‑1 (MQ=255)
                       cGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGTGACAAGGCCATAGTCAAAAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGTGATGACGCCACgg                                                                                                                           >  1:909056/1‑150 (MQ=14)
                            gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATAGAAGGTCACAAGACCATCGTCAATTACTAAGGAGAAATCACCATAACAAAAGGTGGCACCGGCAACCATATTAAAAGTGATGACGCCACgggagac                                                                                                                      >  2:1000745/1‑145 (MQ=17)
                                aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTGAAGAACGAAGGCGAGAGCACCATCCCCAACGGTGGCACCGGCACCCAGAGTAAGGGGGCTGCCGCCACGGCAaacaaca                                                                                                                  >  1:1341749/1‑150 (MQ=17)
                                    gggCACCATGACCGTCACCGATCCTGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                                                                                                              <  1:333076/150‑1 (MQ=255)
                                     ggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                                                                                                             >  2:222168/1‑150 (MQ=255)
                                      gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCa                                                                                                            <  1:146457/150‑1 (MQ=255)
                                         ccATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGGGATGACGCCACGGGAAACAACAACGGCaaaa                                                                                                         >  2:984459/1‑150 (MQ=255)
                                          cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGCGAGCACCATCACCAACGCTGGCACCGGCCCCCAGATTACTGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                        >  1:450023/1‑150 (MQ=255)
                                          cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                        <  2:435652/150‑1 (MQ=255)
                                          cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                        >  1:609570/1‑150 (MQ=255)
                                              acagTCACCGATCAGGATTCCATGGGTATCAAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGAGAGAGCACCATCACCAACGTTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                    <  2:418745/147‑1 (MQ=18)
                                               ccGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                                                                                                   <  1:1295574/150‑1 (MQ=255)
                                               ccGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                                                                                                   <  2:1246942/150‑1 (MQ=255)
                                                 gTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                 >  2:393323/1‑150 (MQ=255)
                                                  tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGCGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                >  2:451011/1‑150 (MQ=255)
                                                    aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                                              >  2:577051/1‑150 (MQ=255)
                                                       gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                                           <  1:756216/150‑1 (MQ=255)
                                                        aTCCGGTGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                                                                                          >  2:516715/1‑150 (MQ=255)
                                                           cGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                                       <  2:1087659/150‑1 (MQ=255)
                                                                   atcggcaTCCAGCGCGACGGAGACAAGGCCATCGTCAACAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                               <  2:1341749/144‑1 (MQ=11)
                                                                         aTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACg                                                                         >  1:1032380/1‑150 (MQ=255)
                                                                              gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGACCCGGCAACCAACAACGGCAAAACCAACGAAGACGGCAAGGATTCCACCGGTACGGAAAt                                                                    >  1:5977/1‑150 (MQ=17)
                                                                                  gACGGTGACAAGGCCATCGTCAATAACGAATGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTATTTGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCATGGATTCCACCGGTACGGAAATCAAc                                                                <  1:1253372/150‑1 (MQ=255)
                                                                                        gACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaat                                                          <  2:222515/150‑1 (MQ=255)
                                                                                          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                                          >  1:110447/1‑132 (MQ=255)
                                                                                          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                                          <  2:110447/132‑1 (MQ=255)
                                                                                            aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACCCAGATTAAAGGTGATGACGCCACCGCAAACAACAACGGCAACCCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATACACGTCAATAACg                                                      >  1:365007/1‑150 (MQ=14)
                                                                                                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGAAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATACCAGGAAATAACGGGAAGGt                                               >  2:727477/1‑150 (MQ=11)
                                                                                                      cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                   <  2:1269216/79‑1 (MQ=9)
                                                                                                      cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                   >  1:1269216/1‑79 (MQ=9)
                                                                                                      cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaaca                                                                                                                     >  2:99401/1‑77 (MQ=9)
                                                                                                      cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaaca                                                                                                                     <  1:99401/77‑1 (MQ=9)
                                                                                                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTTATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                        <  1:222168/150‑1 (MQ=255)
                                                                                                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGTCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                                       >  2:1294439/1‑150 (MQ=255)
                                                                                                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                                       >  1:649487/1‑150 (MQ=255)
                                                                                                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTAGACGGCAAGGATTCCACCGGTACCGAAATCAACGGTAATAACGGAAAAGTTAGTCAgg                                       >  1:231809/1‑150 (MQ=255)
                                                                                                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa                                      <  2:410280/150‑1 (MQ=255)
                                                                                                                         cccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGAt                         <  1:522561/149‑1 (MQ=255)
                                                                                                                          ccATCACCAACTGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACTGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                        <  1:577051/150‑1 (MQ=255)
                                                                                                                              cACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACCACAACGGCAAAACCACCGTTGCCGGCAAGGATTCCACCGGTACGGCAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGCTCGGCATGTCac                    >  1:243357/1‑149 (MQ=17)
                                                                                                                              cACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                    >  1:704268/1‑150 (MQ=255)
                                                                                                                                ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                                             >  1:714651/1‑91 (MQ=37)
                                                                                                                                ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                                             <  2:714651/91‑1 (MQ=37)
                                                                                                                                ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATGAAGGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcg                  >  2:432818/1‑150 (MQ=255)
                                                                                                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc       <  2:1352979/150‑1 (MQ=255)
                                                                                                                                              ggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt    >  1:529275/1‑150 (MQ=255)
                                                                                                                                                caccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                           >  1:1113161/1‑77 (MQ=39)
                                                                                                                                                caccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                           <  2:1113161/77‑1 (MQ=39)
                                                                                                                                                caccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt  <  1:451011/150‑1 (MQ=255)
                                                                                                                                                  cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                        >  2:497798/1‑62 (MQ=39)
                                                                                                                                                  cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                        <  1:497798/62‑1 (MQ=39)
                                                                                                                                                   |                                                                                                                                                  
TGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAT  >  NC_000913/1466189‑1466482

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: