Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,253,193 +C coding (243/258 nt) h16_A2071 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,253,1911.C86.7% 45.7 / 0.7 15coding (241/258 nt)h16_A2071Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (5/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGG  >  AM260479/2253058‑2253251
                                                                                                                                      |                                                             
cTGGCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCaaa                                                                >  2:91125/1‑134 (MQ=255)
           cgcgGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCg                                                    >  2:885541/1‑134 (MQ=255)
             cgGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCt                                                  <  1:885541/134‑1 (MQ=255)
                               cACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                  <  2:893849/132‑1 (MQ=255)
                                                     gCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGcc          <  2:206504/134‑1 (MQ=255)
                                                             gcgcAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAg   <  1:91125/133‑1 (MQ=255)
                                                                                 ttCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                  >  1:670568/1‑82 (MQ=255)
                                                                                 ttCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                  <  2:670568/82‑1 (MQ=255)
                                                                                          tAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACACTGCGCCCGCTGGTAGCGGCCGAGATCGCGCa                             <  1:692461/78‑1 (MQ=255)
                                                                                          tAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACACTGCGCCCGCTGGTAGCGGCCGAGATCGCGCa                             >  2:692461/1‑78 (MQ=255)
                                                                                                          ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCgaga                                    <  1:315205/55‑1 (MQ=255)
                                                                                                          ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCgaga                                    >  2:315205/1‑55 (MQ=255)
                                                                                                               tCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCa                             <  1:485192/57‑1 (MQ=39)
                                                                                                               tCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCa                             >  2:485192/1‑57 (MQ=39)
                                                                                                                  gCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAgg  >  1:296542/1‑81 (MQ=255)
                                                                                                                  gCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAgg  <  2:296542/81‑1 (MQ=255)
                                                                                                                                      |                                                             
CTGGCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGG  >  AM260479/2253058‑2253251

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