Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,464,435 T→C L40P (CTG→CCG)  h16_A1354 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,464,4350TC100.0% 76.1 / NA 21L40P (CTG→CCG) h16_A1354conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/11);  total (10/11)

GATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGC  >  AM260479/1464311‑1464556
                                                                                                                            |                                                                                                                         
gaTGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctg                                                                                                                 <  2:285393‑M1/134‑6 (MQ=255)
    ccATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacg                                                                                                             >  1:512028‑M1/1‑125 (MQ=255)
    ccATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacg                                                                                                             >  1:525485‑M1/1‑125 (MQ=255)
                 tgACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcc                                                                                                 <  1:256727‑M1/133‑22 (MQ=255)
                     cTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                            <  2:388757‑M1/134‑27 (MQ=255)
                                   ggtggGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacgg                                                                              >  1:506153‑M1/1‑94 (MQ=255)
                                           aaCCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggt                                                                      <  1:477985‑M1/134‑49 (MQ=255)
                                              cGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                   <  1:357321‑M1/134‑52 (MQ=255)
                                                          ccgccACACGAGCACCGGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaata                                                         <  1:546828‑M1/133‑62 (MQ=255)
                                                             ccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgg                                                                                                            >  1:130006‑M1/1‑69 (MQ=255)
                                                             ccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgg                                                                                                            <  2:130006‑M1/79‑11 (MQ=255)
                                                                               gCCTCTGGCATTGCGGCCATGAGCAACCGGTGCCAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                   >  2:97506‑M1/1‑51 (MQ=255)
                                                                               gCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                   <  1:97506‑M1/102‑52 (MQ=255)
                                                                                    tGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggt                                <  2:512028‑M1/132‑87 (MQ=255)
                                                                                     ggCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcga                             >  2:59726‑M1/1‑45 (MQ=255)
                                                                                       cATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggttagccgccgaatatgtggcgaatgggcccgcgcccggtcgact                           <  2:1055177‑M1/134‑92 (MQ=255)
                                                                                              gCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctga                    >  1:413323‑M1/1‑36 (MQ=255)
                                                                                                  tGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                            <  2:1013410‑M1/58‑27 (MQ=255)
                                                                                                  tGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGAt                                                                                            >  1:1013410/1‑58 (MQ=21)
                                                                                                        aCCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggc          >  1:32572‑M1/1‑26 (MQ=255)
                                                                                                                 aGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcagg  >  2:676943‑M1/1‑17 (MQ=255)
                                                                                                                            |                                                                                                                         
GATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGC  >  AM260479/1464311‑1464556

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