Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,744,362 G→A 100% H104Y (CAC→TAC)  apt ← adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,744,3620GA95.7% 60.6 / ‑5.0 22H104Y (CAC→TAC) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (13/8);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAC  >  CP000730/1744242‑1744499
                                                                                                                        |                                                                                                                                           
tACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAc                                                                                                                       <  1:108576/141‑1 (MQ=255)
 aCCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAcc                                                                                                                      >  2:38285/1‑141 (MQ=255)
             aCCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTc                                                                                                          >  2:124448/1‑141 (MQ=255)
                 ccTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATAc                                                                                                      >  1:163859/1‑141 (MQ=255)
                                     tATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTACTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTAATACTGATAACGAATGACTTCACg                                                                                  >  1:356675/1‑141 (MQ=255)
                                         gctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGCCTTCACGAGGt                                                                              >  2:160449/1‑141 (MQ=255)
                                         gctgctTCAATCGTACAACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGt                                                                              >  1:234957/1‑141 (MQ=255)
                                          ctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTa                                                                             >  2:121720/1‑141 (MQ=255)
                                                     gTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttct                                                                  >  2:23360/1‑141 (MQ=255)
                                                      tACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttctt                                                                 <  1:38285/141‑1 (MQ=255)
                                                       accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCttt                                                                >  1:361443/1‑141 (MQ=255)
                                                                 gCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGt                                                      >  2:348193/1‑141 (MQ=255)
                                                                       aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTTATACTTATAACGAATGACTTCACGAGGTAATTTCCCTGCTTTTCTAACAGGTGCAAAg                                                <  1:241345/141‑1 (MQ=255)
                                                                         atcatcTGTTATTAACACCCGTCGACCTGGTTGAATAGAAAAAT‑‑TGCATTAGATTTAAACATTTGTCCCATATTCTAGGGCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGcc                                              <  1:89384/141‑1 (MQ=255)
                                                                          tcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCa                                             <  2:98809/141‑1 (MQ=255)
                                                                           catcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCaa                                            >  1:63240/1‑141 (MQ=255)
                                                                                gTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATcccc                                       <  1:121720/141‑1 (MQ=255)
                                                                                      aaCACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTGGTACCATATTCTAGGGCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGaa                                 <  1:103798/141‑1 (MQ=255)
                                                                                              ttGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCTTTGAATAGGCTAc                         >  1:101212/1‑141 (MQ=255)
                                                                                                       tttAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGcc                <  1:380256/141‑1 (MQ=255)
                                                                                                             tGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAt          >  1:396411/1‑141 (MQ=255)
                                                                                                                  cTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGc     <  1:250323/141‑1 (MQ=255)
                                                                                                                     tgtaCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc  <  2:35025/141‑1 (MQ=255)
                                                                                                                        |                                                                                                                                           
TACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAC  >  CP000730/1744242‑1744499

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: