Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,038,183 (T)5→4 intergenic (+78/‑236) BW25113_RS10330 → / → yeeJ tRNA‑Asn/inverse autotransporter adhesin YeeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,038,1790T.82.3% 36.6 / 8.9 17intergenic (+74/‑240)BW25113_RS10330/yeeJtRNA‑Asn/inverse autotransporter adhesin YeeJ
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base . (10/4);  total (12/5)
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.75e-01

TCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCA  >  NZ_CP009273/2038094‑2038244
                                                                                     |                                                                 
tCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGtttta                                                               <  1:427614/90‑2 (MQ=255)
   gaggagCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAAc                                                            >  2:777340/1‑90 (MQ=255)
        gCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTa                                                      >  1:796072/1‑90 (MQ=255)
        gCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGGTTTAACATCTTa                                                      >  1:612291/1‑90 (MQ=255)
          cAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAg                                                    <  2:796072/90‑1 (MQ=255)
                   aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGa                                           <  1:230903/90‑1 (MQ=255)
                      aTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGAcc                                         <  1:777340/90‑1 (MQ=255)
                        tCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCttt                                      >  1:201476/1‑90 (MQ=255)
                          gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTAt                                    <  2:526446/90‑1 (MQ=255)
                              ttttAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGtt                                <  1:256291/90‑1 (MQ=255)
                                                gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGa              >  1:527763/1‑90 (MQ=255)
                                                  ccTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAc            >  1:399981/1‑90 (MQ=255)
                                                  ccTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAc            >  1:197821/1‑90 (MQ=255)
                                                    ttAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGc          >  1:34100/1‑90 (MQ=255)
                                                    ttAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGc          >  2:400001/1‑90 (MQ=255)
                                                     tAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGc          >  2:703469/1‑90 (MQ=255)
                                                     tAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGcc         >  2:608869/1‑90 (MQ=255)
                                                            aCATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCa  >  2:418087/1‑90 (MQ=255)
                                                                                     |                                                                 
TCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCA  >  NZ_CP009273/2038094‑2038244

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