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,0280TC96.3% 83.4 / ‑3.4 27E121G (GAA→GGA) BW25113_RS09050hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/12);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

ATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA  >  NZ_CP009273/1806947‑1807108
                                                                                 |                                                                                
atataATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                          <  2:523229/90‑1 (MQ=255)
                        gTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCa                                                  >  1:331963/1‑90 (MQ=255)
                         tCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCaa                                                 >  2:893086/1‑90 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGa                                                                  <  1:539228/71‑1 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGa                                                                  >  2:539228/1‑71 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                              <  2:151072/75‑1 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                              >  1:151072/1‑75 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTcc                                                   <  1:381468/86‑1 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTcc                                                   >  2:381468/1‑86 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                               >  1:864450/1‑90 (MQ=255)
                           cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                               >  1:964764/1‑90 (MQ=255)
                             tAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                          <  2:976869/61‑1 (MQ=255)
                             tAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                          >  1:976869/1‑61 (MQ=255)
                             tAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATaaa                                             >  2:221640/1‑90 (MQ=255)
                                    tGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGc                                                                            >  1:777803/1‑52 (MQ=255)
                                    tGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGc                                                                            <  2:777803/52‑1 (MQ=255)
                                    tGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCttta                                      <  2:864450/90‑1 (MQ=255)
                                        tAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtttattta                                  >  2:961006/1‑90 (MQ=255)
                                                  gTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTAt                        >  1:935/1‑90 (MQ=255)
                                                             tGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCAc             <  2:389925/90‑1 (MQ=255)
                                                              gATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACt            <  2:908064/90‑1 (MQ=255)
                                                                   aaTTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCAt       >  1:606716/1‑90 (MQ=255)
                                                                     ttGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATaa     <  2:606716/90‑1 (MQ=255)
                                                                      tGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATaaa    <  2:964764/90‑1 (MQ=255)
                                                                        gTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGa  >  1:759407/1‑90 (MQ=255)
                                                                        gTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATAcc                              >  2:235587/1‑62 (MQ=255)
                                                                        gTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATAcc                              <  1:235587/62‑1 (MQ=255)
                                                                                 |                                                                                
ATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA  >  NZ_CP009273/1806947‑1807108

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