Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,593 C→G 10.1% 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
*CP0092731,460,5930CG10.1% 838.4 / 53.5 303pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (160/102);  new base G (12/15);  total (183/120)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 350 total aligned reads displayed.

TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCA  >  CP009273/1460526‑1460660
                                                                   |                                                                   
  gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAAAGGAGGCACCGGCACTCa                                                                   >  1:32840970/1‑68 (MQ=11)
  gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa                                                                   >  1:8940294/1‑68 (MQ=255)
   aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                  >  1:32131641/1‑68 (MQ=255)
   aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGGGGCACCGGCACTCAg                                                                  >  1:17464320/1‑68 (MQ=14)
    cGGTGACAAGGCGTTCGTCAAAAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:30278202/68‑1 (MQ=14)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:23661593/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCAGa                                                                 <  1:30641936/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:2341985/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:30846993/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:28901805/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:25083398/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:24209798/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:12046285/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:553554/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:5140778/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAAAAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:9691554/68‑1 (MQ=17)
    cGGGGACAAGGCCCTCGTAAAAAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:26219390/68‑1 (MQ=14)
    cGGGGACAAGGCCATCGTCAAAAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:21571239/68‑1 (MQ=17)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGAt                                                                >  1:1849640/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                                >  1:21118890/1‑68 (MQ=255)
       tGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                               <  1:32512807/67‑1 (MQ=255)
       tGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                               <  1:836198/67‑1 (MQ=255)
        gACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTaa                                                             <  1:4098744/68‑1 (MQ=255)
        gACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTaa                                                             <  1:8873653/68‑1 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGATTAAt                                                            >  1:23279670/1‑68 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            >  1:14098527/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            <  1:17147833/67‑1 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            <  1:34120860/67‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:33876784/68‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:10923787/68‑1 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtg                                                        >  1:7570502/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGtg                                                        >  1:25255876/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        >  1:26442431/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        >  1:4865368/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCCCCACGCAGAGTAATGGtg                                                        >  1:34309081/1‑68 (MQ=18)
              gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       >  1:3268504/1‑68 (MQ=255)
              gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       >  1:13804922/1‑68 (MQ=255)
              gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       >  1:7816888/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                       <  1:22872331/67‑1 (MQ=14)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                      >  1:3572689/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                      >  1:14240527/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       <  1:27205300/67‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       <  1:15355597/67‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       <  1:5906066/67‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACTCAGATAAAGGGtaag                                                      >  1:10227611/1‑65 (MQ=11)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGGGGCACCGGCACGCAGATAAAAGGtgag                                                      >  1:13473347/1‑67 (MQ=14)
                cATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgatg                                                     <  1:16579754/68‑1 (MQ=255)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                   >  1:7742851/1‑67 (MQ=255)
                     tCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGcc                                                >  1:16003825/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             <  1:31035502/68‑1 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             >  1:257949/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             >  1:27849159/1‑68 (MQ=255)
                          aaCGAAGGCGAGAGCACTATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:13208620/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           >  1:18157821/1‑68 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:28165368/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:11600509/68‑1 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGAGGCCCCGGCACCCAGATTAATCGACAACCCGCCACGGCa                                          >  1:11500595/1‑68 (MQ=14)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         <  1:9495523/68‑1 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         <  1:24628388/68‑1 (MQ=255)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGATTAATGGTGATGACGCCACGGCaa                                         >  1:22310557/1‑67 (MQ=255)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         >  1:15849677/1‑67 (MQ=255)
                                ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGCCGCCACGGCAaaca                                      >  1:24928417/1‑67 (MQ=255)
                                       ggtccATCACAAACCGTGTCACCGGCACACAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                               <  1:18904213/64‑1 (MQ=11)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                               <  1:32413105/67‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                <  1:34493021/66‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                              >  1:29509007/1‑68 (MQ=255)
                                        cACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCa                             >  1:10346723/1‑68 (MQ=255)
                                         aCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                            >  1:26662648/1‑68 (MQ=255)
                                         aCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                            >  1:10792594/1‑68 (MQ=255)
                                         aCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                            >  1:4261357/1‑68 (MQ=255)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                       >  1:9797435/1‑67 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                    >  1:3345544/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:23726647/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:17025012/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:100085/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:15501496/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:25519923/1‑68 (MQ=40)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:19619296/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACgg                                            >  1:16770640/1‑40 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtcg                 >  1:6232036/1‑65 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:10500705/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAAAAACAGGGGAAAAACCACGGTTGa                >  1:29831203/1‑68 (MQ=14)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACCGCAAAAAACACGGGAAAAACCACTGTTGa                >  1:20134806/1‑68 (MQ=14)
                                                      gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAACCCACCGTAGAc               >  1:19707376/1‑68 (MQ=255)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg             <  1:29239097/68‑1 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg              <  1:22664320/67‑1 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg              <  1:12407100/67‑1 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg              <  1:30455831/67‑1 (MQ=40)
                                                         gcaccggcaccCAGATTAATGGTGATGACGCCACGGCTAATAACAGGGGCAAAACCACTGTTGATgga            >  1:20661460/1‑67 (MQ=14)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa          >  1:2859064/1‑68 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa          >  1:27303263/1‑68 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa           <  1:9239182/67‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        >  1:6899085/1‑68 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        <  1:20410897/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        <  1:17770497/68‑1 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt     <  1:14832640/67‑1 (MQ=40)
                                                                  cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        <  1:19276384/63‑1 (MQ=39)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc   >  1:1475538/1‑67 (MQ=37)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc   >  1:31315607/1‑67 (MQ=37)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa  <  1:14377778/68‑1 (MQ=37)
                                                                   |                                                                   
TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCA  >  CP009273/1460526‑1460660

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