Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,056,947 +G intergenic (+211/‑15) h16_A0964 → / → h16_A0965 conserved hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,056,9451.G100.0% 108.1 / NA 29intergenic (+209/‑17)h16_A0964/h16_A0965conserved hypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (15/14);  total (15/14)

AGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTG  >  AM260479/1056815‑1057057
                                                                                                                                   |                                                                                                                
aGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAgg                                                                                                                <  2:403878/133‑1 (MQ=255)
        ttGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGccgcccgc                                                                                                        >  1:473007/1‑134 (MQ=255)
          gCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGccgcccgc                                                                                                        >  2:687427/1‑132 (MQ=255)
          gCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGccgcccgc                                                                                                        >  1:986189/1‑132 (MQ=255)
                    cacaAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              <  2:473007/132‑1 (MQ=255)
                    cacaAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              <  1:687427/132‑1 (MQ=255)
                      caAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcgcg                                                                                          >  1:283850/1‑134 (MQ=255)
                      caAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcgcg                                                                                          <  2:598187/134‑1 (MQ=255)
                                ccgccgGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCAc                                                                                >  1:652952/1‑134 (MQ=255)
                                    cggcggCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAgg                                                                            <  2:283850/134‑1 (MQ=255)
                                     ggcggcGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              <  1:275130/115‑1 (MQ=255)
                                     ggcggcGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              >  2:275130/1‑115 (MQ=255)
                                              ggCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCaa                                                                  <  2:986189/134‑1 (MQ=255)
                                                        gCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAgggc                                                                                                              <  1:405712/79‑3 (MQ=255)
                                                        gCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAgggc                                                                                                              >  2:405712/1‑77 (MQ=255)
                                                              tggCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              <  1:46592/90‑1 (MQ=255)
                                                              tggCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              >  2:46592/1‑90 (MQ=255)
                                                                    cacGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGaa                                            <  1:949574/134‑1 (MQ=255)
                                                                     acGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGccgcccgc                                                                                                        <  1:928257/73‑1 (MQ=255)
                                                                     acGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGccgcccgc                                                                                                        >  2:928257/1‑73 (MQ=255)
                                                                           ggaggGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAActgctg                                      >  2:205831/1‑133 (MQ=255)
                                                                              gggcgggcAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGaaa                                  <  2:652952/134‑1 (MQ=255)
                                                                                                       aCTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCg         >  1:1127171/1‑134 (MQ=255)
                                                                                                            cGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTgg    >  2:625699/1‑134 (MQ=255)
                                                                                                                  ccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              >  2:420280/1‑38 (MQ=25)
                                                                                                                  ccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                              <  1:420280/38‑1 (MQ=25)
                                                                                                                  ccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTg  >  2:571604/1‑130 (MQ=255)
                                                                                                                  ccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTg  <  1:571604/130‑1 (MQ=255)
                                                                                                                         gAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTc                                               >  2:160309/1‑78 (MQ=255)
                                                                                                                         gAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTc                                               <  1:931349/78‑1 (MQ=255)
                                                                                                                         gAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTc                                               >  2:931349/1‑78 (MQ=255)
                                                                                                                         gAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTc                                               <  1:160309/78‑1 (MQ=255)
                                                                                                                                   |                                                                                                                
AGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTG  >  AM260479/1056815‑1057057

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