Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,998,635 +C intergenic (‑221/+527) h16_A3705 ← / ← h16_A3706 transcriptional regulator, LysR‑family/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,998,6331.C100.0% 110.2 / NA 29intergenic (‑219/+529)h16_A3705/h16_A3706transcriptional regulator, LysR‑family/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (16/13);  total (16/13)

ATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCCC  >  AM260479/3998521‑3998761
                                                                                                                 |                                                                                                                                
aTGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCAccc                                                                                                              >  2:65710/1‑134 (MQ=255)
                 aTCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCg                                                                                             >  1:710066/1‑134 (MQ=255)
                 aTCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCg                                                                                             >  1:395459/1‑134 (MQ=255)
                                        gTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCa                                                                                                                 <  2:506954/91‑1 (MQ=255)
                                        gTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCa                                                                                                                 >  1:506954/1‑91 (MQ=255)
                                        gTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGa                                                                       <  2:719956/133‑1 (MQ=255)
                                           tttCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACgg                                                                   >  1:474620/1‑134 (MQ=255)
                                               gtcgAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATAAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAg                                                               >  1:1158979/1‑134 (MQ=255)
                                                  gAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACGGCTACGCCCTGTCGGCTGTGGCCCCCCTCGCCGCCCCCCTGCCGCCCGACCCGGCCTTTCCCAAAAACGGCAAGTTc                                                            >  2:74971/1‑134 (MQ=255)
                                                  gAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGtt                                                             <  1:455099/133‑1 (MQ=255)
                                                  gAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTc                                                            <  1:807682/134‑1 (MQ=255)
                                                            cTCAGGCATCCATTCAATAGAGAAAGCCTGAATCATTAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCa                                                  <  1:74971/134‑1 (MQ=255)
                                                                      cATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGgcgc                                        <  1:930513/134‑1 (MQ=255)
                                                                                agaAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGcc                              <  1:65710/134‑1 (MQ=255)
                                                                                 gaAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTcgccgc                                                                                                       >  1:590350/1‑60 (MQ=255)
                                                                                 gaAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTcgccgc                                                                                                       <  2:590350/60‑1 (MQ=255)
                                                                                 gaAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCa                                                  >  1:64196/1‑113 (MQ=255)
                                                                                 gaAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCa                                                  <  2:64196/113‑1 (MQ=255)
                                                                                   aaGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCa                                                  >  1:241499/1‑111 (MQ=255)
                                                                                   aaGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCa                                                  <  2:241499/111‑1 (MQ=255)
                                                                                   aaGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAg                           >  1:220473/1‑134 (MQ=255)
                                                                                      ccTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTc                        >  1:128749/1‑134 (MQ=255)
                                                                                         gAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGAt                     >  2:476961/1‑134 (MQ=255)
                                                                                              aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGtt                                                             <  1:997753/89‑1 (MQ=255)
                                                                                              aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGtt                                                             >  2:997753/1‑89 (MQ=255)
                                                                                              aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTc                        >  2:341694/1‑126 (MQ=255)
                                                                                              aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTc                        <  1:341694/126‑1 (MQ=255)
                                                                                                           gCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGcc   <  1:476961/134‑1 (MQ=255)
                                                                                                            cTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGccc  >  2:104208/1‑134 (MQ=255)
                                                                                                                 |                                                                                                                                
ATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCCC  >  AM260479/3998521‑3998761

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