Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,302,920 +C coding (1199/1251 nt) dfp → phosphopantothenoylcysteine synthetase/decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,302,9181.C94.1% 57.3 / ‑2.9 17coding (1197/1251 nt)dfpphosphopantothenoylcysteine synthetase/decarboxylase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (7/9);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCC  >  AM260479/3302795‑3303046
                                                                                                                             |                                                                                                                                
ggACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATggc                                                                                                                         >  2:914052/1‑134 (MQ=255)
       aaCGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGgcgcg                                                                                                                  <  2:965383/134‑1 (MQ=255)
                   cTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCa                                                                                                      <  2:53624/134‑1 (MQ=255)
                     gTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATg                                                                                                    >  2:339349/1‑134 (MQ=255)
                         gACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATgcgc                                                                                                 <  1:914052/133‑1 (MQ=255)
                         gACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATgcgc                                                                                                 <  1:339349/133‑1 (MQ=255)
                          aCGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCgg                                                                                               >  2:593337/1‑134 (MQ=255)
                           cgccgcCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGc                                                                                              <  1:593337/134‑1 (MQ=255)
                                    cATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTa                                                                                     >  1:120001/1‑134 (MQ=255)
                                                         cGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTcccc                                                                <  1:1144363/134‑1 (MQ=255)
                                                             aaGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAg                                                             >  1:1111471/1‑133 (MQ=255)
                                                                              ccAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACAGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGc                                            >  1:387474/1‑133 (MQ=255)
                                                                                             gATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATggcg                                                                                                                        >  1:918482/1‑42 (MQ=37)
                                                                                             gATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATggcg                                                                                                                        <  2:918482/42‑1 (MQ=37)
                                                                                                   ccAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAAt                      <  1:579725/134‑1 (MQ=255)
                                                                                                                      acacGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATgcgcgc   >  1:31300/1‑134 (MQ=255)
                                                                                                                       cacGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGcc  <  2:120001/134‑1 (MQ=255)
                                                                                                                             |                                                                                                                                
GGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCC  >  AM260479/3302795‑3303046

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