Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 386,045 T→C pseudogene (339/1527 nt) BW25113_RS01915 → outer membrane autotransporter barrel domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273386,0450TC87.8% 119.4 / 10.6 49pseudogene (339/1527 nt)BW25113_RS01915outer membrane autotransporter barrel domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (22/21);  total (25/24)
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.77e-01

CACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACTG  >  NZ_CP009273/385961‑386130
                                                                                    |                                                                                     
cACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGTTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAt                                                                                  >  2:30455/1‑90 (MQ=255)
 aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                                 <  1:394279/90‑1 (MQ=255)
 aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGCAATTCAGTGATTCACGGCTCGAtt                                                                                 <  1:661472/90‑1 (MQ=255)
           ccTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgc                                                                       <  2:817117/90‑1 (MQ=255)
              attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTAtgctgc                                                                    >  2:787220/1‑90 (MQ=255)
              attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                    >  2:108724/1‑90 (MQ=255)
                   ttGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGc                                                               <  1:54516/90‑1 (MQ=255)
                     gAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGGTTACTTCTATGCTGCCTGgggg                                                             >  2:434650/1‑90 (MQ=255)
                      aaCCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGt                                                            >  2:374521/1‑90 (MQ=255)
                      aaCCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGt                                                            >  1:399649/1‑90 (MQ=255)
                      aaCCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGt                                                            <  2:797038/90‑1 (MQ=255)
                       acagtGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGtt                                                           <  1:612123/86‑1 (MQ=255)
                        ccATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGtta                                                          >  1:399831/1‑90 (MQ=255)
                           tGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGttatta                                                       <  2:49603/90‑1 (MQ=255)
                                aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGTTTACTTCTATGCTGCCCGGCGGTTATTTTGATc                                                  >  2:527040/1‑90 (MQ=255)
                                aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTTTATGCTGCCTGGCGGTTATTATGATc                                                  >  1:95726/1‑90 (MQ=255)
                                aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATc                                                  >  1:430650/1‑90 (MQ=255)
                                  gCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGt                                                <  2:132490/90‑1 (MQ=255)
                                    aGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGttt                                              <  2:443289/90‑1 (MQ=255)
                                     gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGtttt                                             >  2:751690/1‑90 (MQ=255)
                                     gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGGTTACTTCTATGGTGCCTGGGGGGTATTATGGTCgctgt                                             >  2:612123/1‑86 (MQ=255)
                                     gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGtttt                                             >  2:552201/1‑90 (MQ=255)
                                       cAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGa                                           >  1:74768/1‑90 (MQ=255)
                                         gTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATg                                         >  2:69304/1‑90 (MQ=255)
                                          tACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACTGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGc                                        <  1:751690/90‑1 (MQ=255)
                                          tACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGc                                        <  1:552201/90‑1 (MQ=255)
                                            cTGTTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAg                                      <  1:30455/90‑1 (MQ=255)
                                              gCTAATCAGCTTGATATCAGTAATTCAGTGTTTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                    <  2:579596/90‑1 (MQ=255)
                                              gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                    <  1:467499/90‑1 (MQ=255)
                                              gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                    <  2:38099/90‑1 (MQ=255)
                                              gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                    <  2:122585/90‑1 (MQ=255)
                                              gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                    <  1:69304/90‑1 (MQ=255)
                                                   tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTTTGATCGTTTTGATGCAGATGGtaa                               >  2:718562/1‑90 (MQ=255)
                                                   tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTTTGATCGTTTTGATGCAGATGGtaa                               >  2:482477/1‑90 (MQ=255)
                                                   tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaa                               >  1:210598/1‑90 (MQ=255)
                                                   tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaa                               >  2:729954/1‑90 (MQ=255)
                                                     aGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTTTGATCGTTTTGATGCCGATGGtaata                             >  2:157722/1‑90 (MQ=255)
                                                      gCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaataa                            >  2:422661/1‑90 (MQ=255)
                                                      gCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaataa                            >  1:733850/1‑90 (MQ=255)
                                                                aGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGa                  >  1:176958/1‑90 (MQ=255)
                                                                    aTTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg              <  2:275858/90‑1 (MQ=255)
                                                                    aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg              <  2:355917/90‑1 (MQ=255)
                                                                    aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg              <  1:527040/90‑1 (MQ=255)
                                                                           gATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTa       >  1:339456/1‑90 (MQ=255)
                                                                            aTTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTac      <  1:422661/90‑1 (MQ=255)
                                                                            aTTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTac      <  1:37409/90‑1 (MQ=255)
                                                                               cACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACt   <  2:430650/90‑1 (MQ=255)
                                                                               cACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACt   <  1:108724/90‑1 (MQ=255)
                                                                                aCGGCTCGATTACTTCTATGGTGCCTGGCGGGTATTATGATCGTTTTTATGCAGATGGTAATAATCTGGGGGGGTATGATTTTTACACTg  >  2:742438/1‑90 (MQ=255)
                                                                                    |                                                                                     
CACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACTG  >  NZ_CP009273/385961‑386130

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