Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,744,454 C→T 100% G73D (GGC→GAC)  apt ← adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,744,4540CT95.8% 63.7 / ‑3.4 24G73D (GGC→GAC) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (12/11);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTC  >  CP000730/1744318‑1744581
                                                                                                                                        |                                                                                                                               
atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTGTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAAt                                                                                                                             <  2:144557/141‑1 (MQ=255)
    gTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTCGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATcccc                                                                                                                         >  1:68983/1‑141 (MQ=255)
            cacacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTACCAGGGGCAAAGTCAATCCCCCTTGAATa                                                                                                                 >  1:225878/1‑141 (MQ=255)
                                ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGAACAGCCAATGa                                                                                             >  1:107125/1‑141 (MQ=255)
                                 tGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAt                                                                                            <  1:242669/141‑1 (MQ=255)
                                    aTCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaaa                                                                                         >  1:138983/1‑141 (MQ=255)
                                    aTCTTTGTGCATTGTTAAAACAGTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaaa                                                                                         <  1:162678/141‑1 (MQ=255)
                                     tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAg                                                                                        >  1:43674/1‑141 (MQ=255)
                                          gggCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGTGTCATTTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg                                                                                   <  1:279347/139‑1 (MQ=255)
                                            gccttGCTAGAACACTTTTCCCATATTTTAGGTCTTTCTCAGCACGCAGTTTTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcgcg                                                                                 <  1:250355/138‑1 (MQ=255)
                                             cATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                                                <  1:160260/141‑1 (MQ=255)
                                                 gTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCa                                                                            >  1:117390/1‑141 (MQ=255)
                                                    aaaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTAGCATAGTTATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGt                                                                         <  1:8182/141‑1 (MQ=255)
                                                      aaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGTTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTcc                                                                       <  1:122397/141‑1 (MQ=255)
                                                                          gttCATACTCATAACGGATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                   <  2:108220/139‑1 (MQ=255)
                                                                            tCATACTCATAACGAATGACTTAACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAACGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTg                                                 <  2:65322/141‑1 (MQ=255)
                                                                                           aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                                  <  1:175025/141‑1 (MQ=255)
                                                                                              aCTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACa                               >  2:141444/1‑141 (MQ=255)
                                                                                               cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCCACCTCTATGTCTTTTGCGTCTTCTACaa                              >  2:260324/1‑141 (MQ=255)
                                                                                               cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCATGGTGATACTGCTCCTACAACAATATCAACATCTCTGTCTTTTGCGTATTcccgat                              >  1:78452/1‑136 (MQ=255)
                                                                                                 tCACGAGGTAATTTCCCTTCTTTTCTTACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTATCTTTTGCCTATTCTCAAAtt                            >  1:238799/1‑141 (MQ=255)
                                                                                                 acACGAGGTAATTTCCCTTCTTTTCTAACTGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACGTCTCTGTCTTTTGCGTATTCTACAAtt                            >  2:48132/2‑141 (MQ=255)
                                                                                                                          aaaCAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAACGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCtt   <  2:107125/140‑1 (MQ=255)
                                                                                                                           aaCAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGTTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTCTCTGTTGCATCGACATATGCTTc  >  2:223305/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                               
ATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTC  >  CP000730/1744318‑1744581

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