Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,632,437 (T)7→6 intergenic (‑71/+8) BW25113_RS08135 ← / ← BW25113_RS25620 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,632,4310T.94.4% 62.8 / ‑1.0 18intergenic (‑65/+14)BW25113_RS08135/BW25113_RS25620hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (4/13);  total (4/14)
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.01e-01

TCGCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAG  >  NZ_CP009273/1632344‑1632500
                                                                                       |                                                                     
tCGCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttt                                                                     >  1:283767/1‑90 (MQ=255)
      tttGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttttttcat                                                               <  2:122609/90‑4 (MQ=255)
         ggATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCatgat                                                           <  2:146458/90‑1 (MQ=255)
          gATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCatgatg                                                          <  1:112938/90‑1 (MQ=255)
               tCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTAcc                                                     <  2:535517/90‑1 (MQ=255)
               tCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTAcc                                                     <  2:377909/90‑1 (MQ=255)
                 aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                   <  2:182373/90‑1 (MQ=255)
                  cGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAGAAAAACTATGATTTTTTCATGATGTTACCgta                                                  <  2:484961/90‑1 (MQ=255)
                  cGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCgt                                                   <  1:614362/90‑1 (MQ=255)
                     atgatgGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgtagta                                               <  2:160741/90‑1 (MQ=255)
                     atgatgGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgtagta                                               <  2:620110/90‑1 (MQ=255)
                               tAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATc                                     <  1:148537/90‑1 (MQ=255)
                                 aCTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCa                                   <  2:432700/90‑1 (MQ=255)
                                   tCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATc                                 <  2:283767/90‑1 (MQ=255)
                                    ccAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCg                                >  1:242146/1‑90 (MQ=255)
                                         cATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAg                           >  2:57400/1‑90 (MQ=255)
                                               tGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCa                     <  2:5578/90‑1 (MQ=255)
                                                 gCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAg                   >  1:512246/1‑90 (MQ=255)
                                                            gctgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTc        >  1:493999/1‑90 (MQ=255)
                                                                  tCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAg  <  1:678123/90‑1 (MQ=255)
                                                                                       |                                                                     
TCGCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAG  >  NZ_CP009273/1632344‑1632500

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