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,4540CT99.3% 58.8 / ‑6.8 24G73D (GGC→GAC) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/14);  minor base G (1/0);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATT  >  CP000730/1744324‑1744587
                                                                                                                                  |                                                                                                                                     
aaTTAACACACGTTAACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTGCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCa                                                                                                                              >  1:161285/1‑140 (MQ=255)
      cacacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAAt                                                                                                                        <  2:65426/140‑1 (MQ=255)
         acGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACAATTGTACCATACTCTAGGTCATACTCATAACGAATGACTTCACGCGGTAATTTCCCTTCTTTTCTAAGAGGCACAAAGTCAATCCCCATTGAATaaa                                                                                                                     >  2:93440/1‑139 (MQ=255)
                    gTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAg                                                                                                          <  2:215197/140‑1 (MQ=255)
                                    gtgCTTTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                                                                          <  2:211178/140‑1 (MQ=255)
                                     tgCATTGTTAAAACATTTGTACCATAATCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg                                                                                         >  2:174334/1‑140 (MQ=255)
                                              aaaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                                                <  1:279062/140‑1 (MQ=255)
                                              aaaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                                                <  2:193806/140‑1 (MQ=255)
                                                             tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                                                 <  1:66374/140‑1 (MQ=255)
                                                                tCTAGGTCATTCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAAc                                                              <  1:263628/140‑1 (MQ=255)
                                                                   aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGCCAGCCAATGATAAAGCCACGCGCGTCAGGTCCTACAACGCTCTCAACAtc                                                           >  2:261149/1‑140 (MQ=255)
                                                                     gTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                         <  1:139249/140‑1 (MQ=255)
                                                                     gTCATACTCATAACGAATGACTTCACGAGGTAATTTACCTTCTTTTCTACCAGGTGCAAAGTCAATCTCCATTGAATAAGCTACAGGACAGCAAATAATAAAGCAACGCGATTCAGGGCTTATAACGATATCAACAtttt                                                         >  2:131293/1‑137 (MQ=255)
                                                                      tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                                                        <  1:83809/140‑1 (MQ=255)
                                                                          aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCt                                                    <  1:174334/140‑1 (MQ=255)
                                                                             cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                                                          <  1:250856/131‑1 (MQ=255)
                                                                             cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                                                          >  2:250856/1‑131 (MQ=255)
                                                                               tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTCCAACGATATCAACATCTCTGTCTTTTGc                                               >  2:289733/1‑140 (MQ=255)
                                                                                   gAATGACGTCACTAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAt                                           <  2:201087/140‑1 (MQ=255)
                                                                                      tGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                                        >  1:257050/1‑140 (MQ=255)
                                                                                          ttaaCGAGGTAATTTCTCTACTTTTCTAAAAGGTGCAAAGGCAATGCCCATTGAATAAGCTACAGGACAGCCAATGGTAAAGCCACGCGCTTCAAAACATAAAACAATATCAACATCTATGTCTTTTTCGTATTATACaa                                    >  2:166923/4‑140 (MQ=255)
                                                                                                                 ttCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCatat             <  2:161978/140‑1 (MQ=255)
                                                                                                                           tGCAAAGTCAATCCCCATGGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCAt   <  2:96084/140‑1 (MQ=255)
                                                                                                                            gCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTTTCTCTTGCATAGCCATATgagcacctctt  >  2:173724/1‑130 (MQ=255)
                                                                                                                                  |                                                                                                                                     
AATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATT  >  CP000730/1744324‑1744587

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