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,0450TC96.4% 261.3 / 1.3 85pseudogene (339/1527 nt)BW25113_RS01915outer membrane autotransporter barrel domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (43/38);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

ACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTAC  >  NZ_CP009273/385962‑386126
                                                                                   |                                                                                 
aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             <  2:814427/90‑1 (MQ=255)
aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             <  2:6102/90‑1 (MQ=255)
aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             <  2:49817/90‑1 (MQ=255)
aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             <  2:1551922/90‑1 (MQ=255)
aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             <  1:43105/90‑1 (MQ=255)
aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             <  1:356924/90‑1 (MQ=255)
 cGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTa                                                                            >  2:972332/1‑90 (MQ=255)
    tctcAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACtt                                                                         >  2:99494/1‑90 (MQ=255)
            tattatTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctg                                                                 >  1:490062/1‑90 (MQ=255)
             attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                >  1:884256/1‑90 (MQ=255)
             attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                >  2:810708/1‑90 (MQ=255)
             attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                >  2:417238/1‑90 (MQ=255)
             attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                >  2:389518/1‑90 (MQ=255)
             attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                >  1:643368/1‑90 (MQ=255)
             attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                >  1:564107/1‑90 (MQ=255)
             attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                >  1:1453277/1‑90 (MQ=255)
              ttattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTTACGGCTCGATTACTTCTATGCTGcc                                                               >  2:898224/1‑90 (MQ=255)
                attTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGTCTg                                                             <  2:218821/90‑1 (MQ=255)
                attTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTg                                                             <  1:1017170/90‑1 (MQ=255)
                  ttGAACCATGATTATGCAGACAGTTCTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGc                                                           <  2:933405/90‑1 (MQ=255)
                    gAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGgggg                                                         >  1:811646/1‑90 (MQ=255)
                    gAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCgg                                                         <  1:898224/90‑1 (MQ=255)
                     aaCCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGt                                                        >  2:1429357/1‑90 (MQ=255)
                         aTGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGttatt                                                    >  1:343317/1‑90 (MQ=255)
                          tGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             <  1:471226/64‑1 (MQ=255)
                          tGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                             >  2:471226/1‑64 (MQ=255)
                           gATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTtattat                                                  >  1:379384/1‑90 (MQ=255)
                              tATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGAt                                               <  2:43442/90‑1 (MQ=255)
                               aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATc                                              >  2:1173799/1‑90 (MQ=255)
                               aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATc                                              >  1:1389794/1‑90 (MQ=255)
                               aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATc                                              <  1:588200/90‑1 (MQ=255)
                               aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATc                                              >  1:175286/1‑90 (MQ=255)
                                   aGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGttt                                          <  2:643368/90‑1 (MQ=255)
                                   aGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGttt                                          <  1:99494/90‑1 (MQ=255)
                                    gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGtttt                                         >  2:1245868/1‑90 (MQ=255)
                                    gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGtttt                                         >  1:229068/1‑90 (MQ=255)
                                      cAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGa                                       >  1:1027778/1‑90 (MQ=255)
                                       aGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGG‑‑C‑‑‑TCTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                <  1:1403491/90‑1 (MQ=255)
                                         tACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGc                                    <  1:663595/90‑1 (MQ=255)
                                         tACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGc                                    <  2:1453277/90‑1 (MQ=255)
                                         tACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGc                                    <  2:1437410/90‑1 (MQ=255)
                                         tACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGc                                    >  1:1228704/1‑90 (MQ=255)
                                         tACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGc                                    <  2:564107/90‑1 (MQ=255)
                                             gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTTTTATGATCGTTTTGGTGCGGAt                                >  1:1394898/1‑90 (MQ=255)
                                             gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                <  2:343317/90‑1 (MQ=255)
                                             gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                <  2:912209/90‑1 (MQ=255)
                                             gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                <  1:1555275/90‑1 (MQ=255)
                                             gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                <  1:1518285/90‑1 (MQ=255)
                                             gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                <  1:1405564/90‑1 (MQ=255)
                                                aaTCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGt                             >  2:107793/1‑90 (MQ=255)
                                                aaTCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGt                             >  1:58445/1‑90 (MQ=255)
                                                  tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGGTTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaa                           >  1:1391114/1‑90 (MQ=255)
                                                  tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAg                                  <  2:1155407/83‑1 (MQ=255)
                                                  tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAg                                  >  1:1155407/1‑83 (MQ=255)
                                                  tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaa                           >  1:838485/1‑90 (MQ=255)
                                                  tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaa                           >  2:703385/1‑90 (MQ=255)
                                                     gcctGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaataa                        >  1:114396/4‑90 (MQ=255)
                                                     gCTTGATATCAGTAATTCAGTGATTCACGGCTCGGTTAGTTCTATGCTGCCTGGGGGTTTTTATGGTCGTTTTGATGCAGATGGtaataa                        >  1:1345462/1‑90 (MQ=255)
                                                     gCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTTCTTTTATGCTGCCTGGCGGTTATTTTGATCGTTTTGATGCAGATGGtaataa                        >  2:963286/1‑90 (MQ=255)
                                                     gCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaataa                        >  1:195437/1‑90 (MQ=255)
                                                            atCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGt                 <  1:1041094/90‑1 (MQ=255)
                                                              cAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTgg               <  2:379384/90‑1 (MQ=255)
                                                               aGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGa              >  1:787514/1‑90 (MQ=255)
                                                               aGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgc                                                                   >  1:126392/1‑37 (MQ=38)
                                                               aGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgc                                                                   <  2:126392/37‑1 (MQ=38)
                                                               aGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                <  2:168278/72‑1 (MQ=255)
                                                               aGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                >  1:168278/1‑72 (MQ=255)
                                                                  aaTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGatat           >  1:328822/1‑90 (MQ=255)
                                                                   aTTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  1:1096254/90‑1 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGt                 <  2:552453/83‑1 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGt                 >  1:552453/1‑83 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  2:570870/90‑1 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  1:44138/90‑1 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  1:963286/90‑1 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  2:756803/90‑1 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  1:703385/90‑1 (MQ=255)
                                                                   aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  2:1389794/90‑1 (MQ=255)
                                                                        gTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTTTGATCGTTTTTATGCAGATGGTAATAATCTGGGGGGATATGAtttt     >  1:52205/1‑90 (MQ=255)
                                                                        gTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGAtttt     >  2:433423/1‑90 (MQ=255)
                                                                          gATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGata            >  2:578418/1‑81 (MQ=255)
                                                                          gATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGata            <  1:578418/81‑1 (MQ=255)
                                                                          gATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTa   >  2:1254127/1‑90 (MQ=255)
                                                                           aTTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTac  <  2:221447/90‑1 (MQ=255)
                                                                             tCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          <  2:577668/80‑1 (MQ=255)
                                                                             tCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg          >  1:577668/1‑80 (MQ=255)
                                                                                   |                                                                                 
ACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTAC  >  NZ_CP009273/385962‑386126

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