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,0280TC100.0% 134.1 / NA 36E121G (GAA→GGA) BW25113_RS09050hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/17);  total (19/17)

GATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGATTATTA  >  NZ_CP009273/1806942‑1807114
                                                                                      |                                                                                      
gATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTg                                                                                     <  1:395896/90‑1 (MQ=255)
gATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTg                                                                                     <  1:704/90‑1 (MQ=255)
     atataATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                                <  2:2734530/90‑1 (MQ=255)
             ttGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGa                                                                        >  2:2941447/1‑90 (MQ=255)
               gTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAg                                                                      >  2:2026239/1‑90 (MQ=255)
               gTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAg                                                                      >  2:1116854/1‑90 (MQ=255)
                 tCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  2:2509588/90‑1 (MQ=255)
                 tCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  1:992380/90‑1 (MQ=255)
                 tCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  1:926837/90‑1 (MQ=255)
                  cATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGa                                                                   <  1:2928023/90‑1 (MQ=255)
                    tAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACa                                                                 >  2:2735565/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:2055300/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:2566603/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:2862925/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  1:1752829/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:3002263/1‑90 (MQ=255)
                                       aGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:1241337/90‑1 (MQ=255)
                                       aGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:861470/90‑1 (MQ=255)
                                        gTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              >  1:2448512/1‑89 (MQ=255)
                                        gTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:2448512/89‑1 (MQ=255)
                                          gAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGaaa                                                                              >  1:1790050/1‑55 (MQ=255)
                                          gAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGaaa                                                                              <  2:1790050/55‑1 (MQ=255)
                                               accaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                                >  1:1199149/1‑48 (MQ=255)
                                               accaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                                <  2:1199149/48‑1 (MQ=255)
                                                  aaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:2586994/79‑1 (MQ=255)
                                                  aaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              >  1:2586994/1‑79 (MQ=255)
                                                  aaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACct                                   >  2:1282950/1‑90 (MQ=255)
                                                           tttATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGa                          <  2:1430890/90‑1 (MQ=255)
                                                                 cTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  2:2588371/42‑1 (MQ=255)
                                                                 cTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    >  1:2588371/1‑42 (MQ=255)
                                                                 cTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCa                    >  2:1300858/1‑90 (MQ=255)
                                                                  tGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCAc                   >  1:1602610/1‑90 (MQ=255)
                                                                   gATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCAt                                                            >  1:2920675/1‑48 (MQ=255)
                                                                   gATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCAt                                                            <  2:2920675/48‑1 (MQ=255)
                                                                   gATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACt                  >  2:2480338/1‑90 (MQ=255)
                                                                           tGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATaaa          <  1:1033961/90‑1 (MQ=255)
                                                                                   tttcCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGattatta  <  2:2388449/90‑1 (MQ=255)
                                                                                      |                                                                                      
GATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGATTATTA  >  NZ_CP009273/1806942‑1807114

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