Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_017635 2,292,091 +G coding (36/39 nt) ECW_RS27980 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,292,0911.G100.0% 54.3 / NA 28V12V (GTA→GTCECW_RS27980hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (9/19);  total (9/19)

AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTT  >  NC_017635/2292042‑2292134
                                                  |                                           
agagCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTtat                                              >  1:130474/1‑50 (MQ=255)
agagCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTtat                                              >  2:512085/1‑50 (MQ=255)
 gagCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTtat                                              <  1:1232153/49‑1 (MQ=255)
      gggTTCCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTaa                                        <  2:1709490/50‑1 (MQ=17)
        gTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAaga                                      <  2:809898/50‑1 (MQ=17)
            cccccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGca                                   <  1:1988563/49‑1 (MQ=17)
              cccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGcacca                                <  1:1893003/50‑1 (MQ=17)
              cccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGcacca                                >  1:2013267/1‑50 (MQ=17)
                cGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACg                              <  2:1186205/50‑1 (MQ=17)
                 gCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACga                             <  1:841354/50‑1 (MQ=17)
                   ccTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                           <  1:724434/50‑1 (MQ=17)
                     tGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatgat                         >  2:1429234/1‑50 (MQ=17)
                        tAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAg                      >  1:194484/1‑50 (MQ=17)
                          gTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTa                    <  1:709783/50‑1 (MQ=17)
                          gTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTa                    <  2:851590/50‑1 (MQ=17)
                          gTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTa                    <  2:2054033/50‑1 (MQ=17)
                               aGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCt                 <  1:183666/48‑1 (MQ=17)
                                gTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtca              <  1:794160/50‑1 (MQ=17)
                                 tAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtcat             <  1:1907547/50‑1 (MQ=17)
                                  aaGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtcat             >  2:1114226/1‑49 (MQ=17)
                                   aGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtcatca           <  1:396579/50‑1 (MQ=17)
                                   aGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtcatca           <  1:294359/50‑1 (MQ=17)
                                     cGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTcatcat          >  1:1897140/1‑49 (MQ=14)
                                     cGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATg         <  1:1812204/50‑1 (MQ=14)
                                      ggggAAGCTTATGACTAAGAGCACCACGATGATGATTAGCTTCATCATGa        <  2:990173/50‑1 (MQ=14)
                                      ggggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGa        <  1:2033675/50‑1 (MQ=14)
                                      ggggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGa        >  2:1592104/1‑50 (MQ=14)
                                         gAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGAccc     >  2:2031298/1‑50 (MQ=11)
                                          aaGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCt    >  1:610709/1‑50 (MQ=11)
                                          aaGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCt    >  2:1143472/1‑50 (MQ=11)
                                            gCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCttt  <  1:1737219/50‑1 (MQ=1)
                                                  |                                           
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTT  >  NC_017635/2292042‑2292134

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