Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,807,028 T→C E121G (GAA→GGA)  BW25113_RS09050 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,807,0280TC92.9% 33.7 / ‑4.1 14E121G (GAA→GGA) BW25113_RS09050hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (5/8);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA  >  NZ_CP009273/1806944‑1807108
                                                                                    |                                                                                
tGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTggg                                                                             <  1:335245/90‑1 (MQ=255)
 gTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGc                                                                            <  1:196026/90‑1 (MQ=255)
   atataATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                          <  1:445428/90‑1 (MQ=255)
             gTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAg                                                                >  1:410513/1‑90 (MQ=255)
                              cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                               >  2:30992/1‑90 (MQ=255)
                                        gAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTc                                                    <  1:187180/75‑1 (MQ=255)
                                        gAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTc                                                    >  2:187180/1‑75 (MQ=255)
                                             accaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCttt                                       <  1:124918/83‑1 (MQ=255)
                                             accaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCttt                                       >  2:124918/1‑83 (MQ=255)
                                              ccaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATAc                               <  2:359925/90‑1 (MQ=255)
                                               caaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATAcc                              <  1:198386/90‑1 (MQ=255)
                                                               cTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATACAATCTTTATTTATACCTCTTATCCGAATTTCa              >  2:7228/1‑90 (MQ=255)
                                                                      aaTTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCAt       >  2:441322/1‑90 (MQ=255)
                                                                           gTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGa  <  2:202150/90‑1 (MQ=255)
                                                                                    |                                                                                
TGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA  >  NZ_CP009273/1806944‑1807108

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