Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,253,193 +C coding (243/258 nt) h16_A2071 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,253,1911.C100.0% 65.9 / NA 18T81P (ACC→CCC) h16_A2071Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (8/10);  total (8/10)

GCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACA  >  AM260479/2253061‑2253281
                                                                                                                                   |                                                                                          
gCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAAccc                                                                                         <  2:133820/134‑2 (MQ=255)
                            cACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                <  2:487705/132‑1 (MQ=255)
                                    cACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGcc                                                                      >  2:469143/1‑118 (MQ=255)
                                    cACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGcc                                                                      <  1:469143/118‑1 (MQ=255)
                                        gATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGGGATCGCGCAGTACCgaga                                                  >  1:642447/1‑134 (MQ=255)
                                                   ccACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                >  2:246697/1‑109 (MQ=255)
                                                   ccACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                <  1:246697/109‑1 (MQ=255)
                                                   ccACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGccc                                       >  1:469213/1‑134 (MQ=255)
                                                                      cTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAg                    <  1:809148/134‑1 (MQ=255)
                                                                         tgCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGc                   <  2:469213/132‑1 (MQ=255)
                                                                            aTTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAAccc                                                                                         >  1:42310/1‑57 (MQ=255)
                                                                            aTTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAAccc                                                                                         <  2:42310/58‑2 (MQ=255)
                                                                              ttCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGAt                                                                 >  1:453570/1‑81 (MQ=255)
                                                                              ttCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGAt                                                                 <  2:453570/81‑1 (MQ=255)
                                                                                    cGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGt                >  1:360536/1‑124 (MQ=255)
                                                                                    cGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGt                <  2:360536/124‑1 (MQ=255)
                                                                                         acaacaacaGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAaca  <  1:151633/133‑1 (MQ=255)
                                                                                                       ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                >  1:55577/1‑57 (MQ=255)
                                                                                                       ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                <  2:55577/57‑1 (MQ=255)
                                                                                                                    cGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCgaga                                                  <  2:994216/58‑1 (MQ=255)
                                                                                                                    cGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCgaga                                                  >  1:994216/1‑58 (MQ=39)
                                                                                                                                   |                                                                                          
GCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACA  >  AM260479/2253061‑2253281

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