Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,484,856 A→T 100% pseudogene (115/750 nt) gapC ← pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,484,8560AT100.0% 105.4 / NA 29pseudogene (115/750 nt)gapCpseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/15);  total (14/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

CGCTTTCCACGGAATATTTTTAGCCTCTTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTATGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACCCAACACCAGTCGACCGATACGACCAAAACCGTTAATACCAACTTTACTCATGGTTTTCTCCTGTCAGGAAC  >  CP009273/1484710‑1484990
                                                                                                                                                  |                                                                                                                                      
cGCTTTCCACGGAATATTTTTAGCCTCTTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTgg                                                                                                                                     >  2:176785/1‑150 (MQ=255)
 gCTTTCCACGGAATATTTTTAGCCTCTTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTggg                                                                                                                                    >  2:128837/1‑150 (MQ=255)
       cACGGAATATTTTTAGCCTCTTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGt                                                                                                                               >  1:131135/1‑149 (MQ=255)
         cGGAATATTTTTAGCCTCTTTGTCGGCGTAAACCGCGATGCTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTgag                                                                                                                            <  1:194431/150‑1 (MQ=255)
            aaTATTTTTAGCCTCTTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGAt                                                                                                                          >  2:169136/1‑149 (MQ=255)
                        ctctTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGAc                                                                                                             <  1:43738/150‑1 (MQ=255)
                        ctctTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGCAGGCGAGAATTTTTggg                                                                                                                                    >  1:2052/1‑127 (MQ=255)
                         tctTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTgggg                                                                                                                                   <  2:2052/127‑1 (MQ=255)
                                 ggCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTAt                                                                                                    <  1:176785/150‑1 (MQ=255)
                                    gTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGtt                                                                                                 <  2:14587/150‑1 (MQ=255)
                                        aCCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCtt                                                                                             >  1:57102/1‑150 (MQ=255)
                                             gATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGAc                                                                                        >  1:12844/1‑150 (MQ=255)
                                             gATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGGAGTTTGAATCATGTTTCAGCAGGGAGGCGAGAATTTTTGGGGAAGTGAGATCATTATTAGCGACAACGTCTATGTTGCTTTTGAc                                                                                        >  2:74598/1‑150 (MQ=255)
                                                    ttCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAg                                                                                  <  2:8632/149‑1 (MQ=255)
                                                    ttCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAg                                                                                  <  2:56133/149‑1 (MQ=255)
                                                      cccATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTa                                                                                <  1:35324/149‑1 (MQ=255)
                                                       ccATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAAt                                                                              <  1:128837/150‑1 (MQ=255)
                                                                    gTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCa                                                                 <  2:161603/150‑1 (MQ=255)
                                                                        aTCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTc                                                              >  2:169223/1‑149 (MQ=255)
                                                                             ccGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGAcc                                                          <  1:169136/148‑1 (MQ=255)
                                                                                tAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGACCGATAc                                                     <  2:72980/150‑1 (MQ=255)
                                                                                  aaaTCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGACCGATACGa                                                   <  2:131135/150‑1 (MQ=255)
                                                                                              cAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGACCGATACGACCAAAACCGTTa                                       <  1:193468/150‑1 (MQ=255)
                                                                                                     aTGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGACCGATACGACCAAAACCGTTAATACCaa                                >  1:182316/1‑150 (MQ=255)
                                                                                                            gTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGACCGATACGACCAAAACCGTTAATACCAACTTTACt                         >  1:109796/1‑150 (MQ=255)
                                                                                                                 ttGAATCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGACCGATACGACCAAAACCGTTAATACCAACTTTACTCATgg                    >  1:182768/1‑150 (MQ=255)
                                                                                                                     aTCATGTTTCAGCAGGTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACATCAGCCGACCGATACGACCAAAACAGTTAATACCAACTTTACTCATGGtttt                >  1:160239/1‑150 (MQ=255)
                                                                                                                                   ggTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTCGACCGATACGACCAAAACCGTTAATACCAACTTTACTCATGGTTTTCTCCTGTCAGGAAc  <  2:146557/150‑1 (MQ=255)
                                                                                                                                   ggTAGGCGAGAATTTTTGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACGCAACACCAGTAGACCGATACGACCAAAACCGTTAATACCAACTTTACTCATGGTTTTCTCCTGTCAGGAAc  >  2:134233/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                      
CGCTTTCCACGGAATATTTTTAGCCTCTTTTTCGGCGTAAACCGCGATACTTTTCCCATCAACGATAAGTGAATCTTCCGTAAAATCAACGCTCCAGGGGAATGGTCCGTAGTTTGAATCATGTTTCAGCAGGTAGGCGAGAATTTATGGGGAAGTGAGATCATTAATAGCGACAACGTCTATGTTGCTTTTGACTTCAAGTAATCGACCCAACACCAGTCGACCGATACGACCAAAACCGTTAATACCAACTTTACTCATGGTTTTCTCCTGTCAGGAAC  >  CP009273/1484710‑1484990

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