Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,972,184 +C intergenic (‑136/‑395) h16_A3682 ← / → h16_A3683 putative methyltransferase/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,972,1821.C100.0% 192.0 / NA 49intergenic (‑134/‑397)h16_A3682/h16_A3683putative methyltransferase/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (25/24);  total (25/24)

CGTGCTCGTCGCCTGTTTGCCTGAGTGGCCCGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCACGCGGATTGGAT  >  AM260479/3972058‑3972299
                                                                                                                             |                                                                                                                     
cGTGCTCGTCGCCTGTTTGCCTGAGTGGCCCGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGtt                                                                                                                <  1:999347/133‑1 (MQ=255)
 gTGCTCGTCGCCTGTTTGCCTGAGTGGCCCGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGc                                                                                                              <  1:1084300/134‑1 (MQ=255)
    cTCGTCGCCTGTTTGCCTGAGTGGCCCGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           <  2:245957/134‑1 (MQ=255)
        tcgCCTGTTTGCCTGAGTGGCCCGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           <  2:994242/130‑1 (MQ=255)
        tcgCCTGTTTGCCTGAGTGGCCCGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           >  1:994242/1‑130 (MQ=255)
                             ccGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           <  2:438911/109‑1 (MQ=255)
                             ccGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           >  1:438911/1‑109 (MQ=255)
                                  gtgCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           <  1:249759/104‑1 (MQ=255)
                                  gtgCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           >  2:249759/1‑104 (MQ=255)
                                                         ggggATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGt                                                       <  2:993375/133‑1 (MQ=255)
                                                                     ggTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGa                                          >  2:1002259/1‑134 (MQ=255)
                                                                                  cacGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGc                                                                                                              >  1:400824/1‑53 (MQ=255)
                                                                                  cacGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGc                                                                                                              <  2:400824/53‑1 (MQ=255)
                                                                                   acGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTcccc                                                   <  1:286239/111‑1 (MQ=255)
                                                                                   acGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTcccc                                                   >  2:286239/1‑111 (MQ=255)
                                                                                          cGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTg                                                                                                         <  1:1117092/50‑1 (MQ=255)
                                                                                          cGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTg                                                                                                         >  2:1117092/1‑50 (MQ=255)
                                                                                                 ccccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCa                                                                           >  2:791574/1‑73 (MQ=255)
                                                                                                 ccccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCa                                                                           <  1:791574/73‑1 (MQ=255)
                                                                                                 ccccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCa                                                                           <  1:27352/73‑1 (MQ=255)
                                                                                                 ccccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCa                                                                           >  2:27352/1‑73 (MQ=255)
                                                                                                 ccccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTcccc                                                   >  1:1059309/1‑97 (MQ=255)
                                                                                                 ccccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTcccc                                                   <  2:1059309/97‑1 (MQ=255)
                                                                                                 ccccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCa              >  1:175325/1‑134 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGa                                                                                                            >  2:485979/1‑38 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGa                                                                                                            <  1:485979/38‑1 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           <  2:780261/39‑1 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGAt                                                                                                           >  1:780261/1‑39 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCcgcg                                                                     <  1:158547/77‑1 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCcgcg                                                                     >  2:606412/1‑77 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCcgcg                                                                     >  2:158547/1‑77 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCcgcg                                                                     <  1:606412/77‑1 (MQ=255)
                                                                                                   ccTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCAcg            >  1:344139/1‑134 (MQ=255)
                                                                                                    cTTTCCGGCGACCTCACCGTATACGCCCTGGTTGc                                                                                                              <  1:31345/35‑1 (MQ=255)
                                                                                                    cTTTCCGGCGACCTCACCGTATACGCCCTGGTTGc                                                                                                              >  2:31345/1‑35 (MQ=255)
                                                                                                    cTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGAt                                         <  2:563362/104‑1 (MQ=255)
                                                                                                    cTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGAt                                         >  1:563362/1‑104 (MQ=255)
                                                                                                    cTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATa                                        <  1:868481/105‑1 (MQ=255)
                                                                                                    cTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATa                                        >  2:868481/1‑105 (MQ=255)
                                                                                                    cTTTCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCAcgc           >  1:36521/1‑134 (MQ=255)
                                                                                                      ttCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACg                                                        <  1:886111/87‑1 (MQ=255)
                                                                                                      ttCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACg                                                        >  2:886111/1‑87 (MQ=255)
                                                                                                      ttCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACtg                        <  1:333552/119‑1 (MQ=255)
                                                                                                      ttCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACtg                        >  2:333552/1‑119 (MQ=255)
                                                                                                      ttCCGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCAcgcg          >  2:852308/1‑133 (MQ=255)
                                                                                                        ccGGCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCACGCGGAt       >  1:611487/1‑134 (MQ=255)
                                                                                                           gCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTcc                                                     >  2:526107/1‑85 (MQ=255)
                                                                                                           gCGACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTcc                                                     <  1:526107/85‑1 (MQ=255)
                                                                                                             gACCTCACCGTATACGCCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCACGCGGATTGGAt  <  2:611487/134‑1 (MQ=255)
                                                                                                                             |                                                                                                                     
CGTGCTCGTCGCCTGTTTGCCTGAGTGGCCCGCTGTGCCGGCCGTTCCGGCGCGATCGGGGATGATACCGGTGCAACCGTGACACGTATACGTAACACCCCTTTCCGGCGACCTCACCGTATACGCCTGGTTGCGATTGACACGTATACGGCACCCCATTTCCGCCTCAGCCGCGACGCCCGTTTACGTCCCCGCGTGTGGATAACCCCGCCTCGAACTGTGGACAAGCACGCGGATTGGAT  >  AM260479/3972058‑3972299

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: