Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,592,137 A→G pseudogene (207/5422 nt) BW25113_RS07915 ← ECSE_1600 family autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,592,1370AG89.8% 125.9 / 8.4 49pseudogene (207/5422 nt)BW25113_RS07915ECSE_1600 family autotransporter
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (27/17);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCGTTGTC  >  NZ_CP009273/1592052‑1592224
                                                                                     |                                                                                       
cagccaAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCAt                                                                                     >  1:413072/1‑90 (MQ=255)
     aaaaTATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCa                                                                                <  1:308321/90‑1 (MQ=255)
      aaaTATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCaa                                                                               <  2:682547/90‑1 (MQ=255)
                cACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACt                                                                     <  2:1248283/90‑1 (MQ=255)
                  cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg                                                                   >  1:1191018/1‑90 (MQ=255)
                  cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg                                                                   >  1:1289646/1‑90 (MQ=255)
                  cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg                                                                   >  2:1329021/1‑90 (MQ=255)
                  cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg                                                                   >  1:1417508/1‑90 (MQ=255)
                  cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg                                                                   >  2:1040886/1‑90 (MQ=255)
                       gTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCa                                                              <  2:1346511/90‑1 (MQ=255)
                         gCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGa                                                            >  1:324723/1‑90 (MQ=255)
                          ccATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGAc                                                           <  2:849362/90‑1 (MQ=255)
                            aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc                                                         >  2:1201522/1‑90 (MQ=255)
                            aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc                                                         >  1:630803/1‑90 (MQ=255)
                            aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc                                                         >  1:817273/1‑90 (MQ=255)
                            aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc                                                         >  1:852657/1‑90 (MQ=255)
                            aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc                                                         >  2:909598/1‑90 (MQ=255)
                            aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc                                                         <  2:532113/90‑1 (MQ=255)
                            aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGc                                                         >  1:453401/1‑90 (MQ=255)
                              aCCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTg                                                       <  1:1307297/90‑1 (MQ=255)
                               ccAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGa                                                      <  2:662626/90‑1 (MQ=255)
                                 aGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCaa                                                                               >  2:1435986/1‑63 (MQ=255)
                                 aGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCaa                                                                               <  1:1435986/63‑1 (MQ=255)
                                      aCCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAAcc                                                   <  1:70325/86‑1 (MQ=255)
                                      aCCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAAcc                                                   >  2:70325/1‑86 (MQ=255)
                                           gTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGc                                          <  2:166438/90‑1 (MQ=255)
                                            tAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGcc                                         <  2:26832/90‑1 (MQ=255)
                                                catcatAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCa                                     >  1:396984/1‑90 (MQ=255)
                                                   catAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCaaaa                                  >  1:696220/1‑90 (MQ=255)
                                                        gCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACAccc                             <  1:794306/90‑1 (MQ=255)
                                                           aCATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc                          >  1:838681/1‑90 (MQ=255)
                                                           aCATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc                          >  2:888380/1‑90 (MQ=255)
                                                           aCATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc                          >  2:45649/1‑90 (MQ=255)
                                                            cATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCg                         >  2:699348/1‑90 (MQ=255)
                                                             aTCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGt                        <  1:543399/90‑1 (MQ=255)
                                                                aGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCaa                     <  1:1129840/90‑1 (MQ=255)
                                                                 gTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCa                                                              >  2:225707/1‑48 (MQ=255)
                                                                 gTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCa                                                              <  1:225707/48‑1 (MQ=255)
                                                                 gTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAAt                    >  2:201409/1‑90 (MQ=255)
                                                                     cAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTAc                <  1:485179/90‑1 (MQ=255)
                                                                      aaTATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACt               >  2:323268/1‑90 (MQ=255)
                                                                      aaTATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACt               >  1:1149860/1‑90 (MQ=255)
                                                                           tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt          >  2:1504533/1‑90 (MQ=255)
                                                                           tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt          >  2:1316811/1‑90 (MQ=255)
                                                                           tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt          >  1:161369/1‑90 (MQ=255)
                                                                           tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt          >  2:959092/1‑90 (MQ=255)
                                                                            ggTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAttt         >  1:359127/1‑90 (MQ=255)
                                                                              tAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCg       <  1:1040886/90‑1 (MQ=255)
                                                                                   tGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCGTTGTc  <  2:1462041/90‑1 (MQ=255)
                                                                                     |                                                                                       
CAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCGTTGTC  >  NZ_CP009273/1592052‑1592224

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