Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 160,041 G→C G102A (GGG→GCG)  h16_A0146 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479160,0410GC100.0% 57.7 / NA 17G102A (GGG→GCG) h16_A0146conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (9/8);  total (9/8)

GCCACCAGCCCGCAGGGCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGGGGTACGACCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACC  >  AM260479/159921‑160162
                                                                                                                        |                                                                                                                          
gCCACCAGCCCGCAGGGCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGAcccg                                                                                                                <  2:835529/132‑3 (MQ=255)
                     aCCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTa                                                                                          >  1:309373/1‑134 (MQ=255)
                                   gggCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTc                                                                             <  2:309373/133‑1 (MQ=255)
                                                gAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACccgccgccg                                                                                                           <  2:723288/90‑1 (MQ=255)
                                                gAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACccgccgccg                                                                                                           >  1:723288/1‑90 (MQ=255)
                                                                  gacgccgccGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCg                                             >  1:1060228/4‑134 (MQ=255)
                                                                      ccgccgCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTgg                                         >  2:177843/1‑134 (MQ=255)
                                                                                           gcggGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTAc                    >  2:26740/1‑134 (MQ=255)
                                                                                             gggCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATc                                               >  1:673863/1‑105 (MQ=255)
                                                                                             gggCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATc                                               <  2:673863/105‑1 (MQ=255)
                                                                                                     gggCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGcc                                                    <  1:429203/92‑1 (MQ=255)
                                                                                                     gggCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGcc                                                    >  2:429203/1‑92 (MQ=255)
                                                                                                             tCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAAcc  <  2:1060228/134‑1 (MQ=255)
                                                                                                               gATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGaa     >  2:294494/1‑129 (MQ=255)
                                                                                                               gATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGaa     <  1:294494/129‑1 (MQ=255)
                                                                                                                 tGGCGTGCGCTACGACCCGCCGCCGTCGGcgacga                                                                                                 <  2:1048747/35‑1 (MQ=12)
                                                                                                                 tGGCGTGCGCTACGACCCGCCGCCGTCGGcgacga                                                                                                 >  1:1048747/1‑35 (MQ=9)
                                                                                                                        |                                                                                                                          
GCCACCAGCCCGCAGGGCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGGGGTACGACCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACC  >  AM260479/159921‑160162

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