Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,639 C→T 16.4% pseudogene (1248/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,6390CT16.4% 103.5 / 11.5 50pseudogene (1248/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (32/7);  new base T (3/4);  total (36/14)
Fisher's exact test for biased strand distribution p-value = 4.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466491‑1466788
                                                                                                                                                    |                                                                                                                                                     
ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTAACAACGAAGGCGAGAGCGCAATCACCAACGGTGGCACCGGCACGCAGATAAACg                                                                                                                                                      >  2:948807/1‑150 (MQ=255)
          cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACTCAGATCAACGGTCACGACGc                                                                                                                                            >  2:1118379/1‑150 (MQ=14)
           gcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGcc                                                                                                                                           >  1:1283216/1‑150 (MQ=255)
           gcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGcc                                                                                                                                           >  1:1283228/1‑150 (MQ=255)
                                  cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                    >  1:448906/1‑150 (MQ=255)
                                    tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                                                                                                                  >  2:352835/1‑150 (MQ=255)
                                        cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCGAACAACAGCGGCAAAACCACCGt                                                                                                              >  1:194263/1‑150 (MQ=255)
                                         gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTCAAGACGCCACCGCGAAAAACAAAGGAAAAACCACAGtt                                                                                                             >  1:61762/1‑150 (MQ=14)
                                          tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTATCGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTg                                                                                                            >  1:956548/1‑150 (MQ=255)
                                                      aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                                <  2:1283228/150‑1 (MQ=255)
                                                      aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                                <  2:1283216/150‑1 (MQ=255)
                                                           aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCGGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAAGGCAAGGATTCCAcc                                                                                           >  1:166286/1‑150 (MQ=255)
                                                            tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                          <  1:275998/150‑1 (MQ=255)
                                                                       gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAtt                                                                               >  2:1406310/1‑150 (MQ=255)
                                                                        ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAGCGGTACGGAAATTg                                                                              >  2:427090/1‑150 (MQ=255)
                                                                           aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                           >  1:200616/1‑150 (MQ=255)
                                                                             ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTGAAGACGCCACCGCAAAAAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGCAAATAGCTGGc                                                                         >  2:909619/1‑150 (MQ=17)
                                                                                   cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGAACGCAGATTAACGGTCATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATACCTGGCAATatc                                                                   >  2:527134/1‑148 (MQ=17)
                                                                                   cAGGCGGTTGTTAACAACGAAGGCGAGAGCACCATCACCAAAGGCGGCACCGGCACTCAGATAAACGGTAAAGACGCCACCGCTAAAAACAGAGGAAAAACCACAGTTGAAGGAAATGATTCCACCGGTACCAAAATAGCTGGCAATatc                                                                   >  2:527143/1‑148 (MQ=12)
                                                                                         gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAg                                                             >  2:1297656/1‑150 (MQ=255)
                                                                                         gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCACTAACGGGAAg                                                             >  2:1297663/1‑150 (MQ=255)
                                                                                         gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGAAg                                                             >  2:1297679/1‑150 (MQ=255)
                                                                                             ttAACAACGAAGGCGAGAGCGCCATCACCAACAGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                                                         >  2:301111/1‑150 (MQ=255)
                                                                                                ataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGACATGGCGGGCAATAACGGGACGGTGATTc                                                      >  2:2538/3‑150 (MQ=11)
                                                                                                  aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                                                    >  2:610480/1‑150 (MQ=255)
                                                                                                  aacGAAGGCGAGAGCGCAATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATCACGCCGCGGCAAACAACTACTTCACAACCACCGTTGAAGGCTAGGATTCCACCGGTACGGAAATTGCTGGCGATCTCGGGCGGGTGATTCAg                                                    >  2:837151/1‑150 (MQ=14)
                                                                                                     gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                                 >  1:965733/1‑150 (MQ=255)
                                                                                                     gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                                 <  1:1379651/150‑1 (MQ=255)
                                                                                                           gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACGGCAATCAACATAGGCAAAACCACCGTTGTCGGCAAGGATTCCACCGGTACGGAAATTGCTTGCAATGACGGGCAGGTGATTCTGGACGGCGAt                                           >  2:741562/1‑150 (MQ=17)
                                                                                                              cgctacATCAGCAGCGGTTGCACCGGACCGGAGATTAACTGCGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTTCTGGCAATAACGGGAATGTGATTCAGGACGGCGATCTg                                        <  1:80485/145‑1 (MQ=17)
                                                                                                                 gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     >  2:1019784/1‑150 (MQ=255)
                                                                                                                 cccATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  2:632814/149‑1 (MQ=255)
                                                                                                                 accATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  1:802248/149‑1 (MQ=255)
                                                                                                                 accATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  2:500041/149‑1 (MQ=255)
                                                                                                                   cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                                   >  1:558573/1‑150 (MQ=255)
                                                                                                                    aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                                  >  1:581916/1‑150 (MQ=255)
                                                                                                                    aTCACCAACGGAGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                                  >  1:11102/1‑150 (MQ=255)
                                                                                                                        ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                            >  2:1009601/1‑136 (MQ=255)
                                                                                                                        ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                            <  1:1009601/136‑1 (MQ=255)
                                                                                                                        ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                            <  1:1009489/136‑1 (MQ=255)
                                                                                                                        ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                            >  2:1009489/1‑136 (MQ=255)
                                                                                                                             ggTGGCACCGGCACGCAGATTACGGGTGATGCCGCCACGGCAAACAACAACGGCAAAACCCCCGTTGACGGCAAGGATGCCCCCGGTCCGGAAAGTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggc                         <  1:2538/150‑1 (MQ=17)
                                                                                                                             ggTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGAGGGCAATACCGGCAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggc                         <  1:2540/150‑1 (MQ=14)
                                                                                                                                 gCACCGGCACTCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGATGGCAATAACGGGAAGGTGATTCAGGACGGAGATCTGGATGTCAGAGGCGGCGGTc                     >  2:1323121/1‑150 (MQ=255)
                                                                                                                                 gCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                     >  2:1323136/1‑150 (MQ=255)
                                                                                                                                     cGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                 >  1:659520/1‑150 (MQ=255)
                                                                                                                                       gCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa               <  1:301111/150‑1 (MQ=255)
                                                                                                                                          cGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg            >  1:690300/1‑150 (MQ=255)
                                                                                                                                            cAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGa                                                               <  2:147449/97‑1 (MQ=255)
                                                                                                                                            cAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGa                                                               >  1:147449/1‑97 (MQ=255)
                                                                                                                                                    cGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg  >  1:900336/1‑150 (MQ=255)
                                                                                                                                                    cGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg  >  1:900339/1‑150 (MQ=255)
                                                                                                                                                    cGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg  >  1:900401/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466491‑1466788

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