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,3620GA83.3% 64.3 / 7.6 30H104Y (CAC→TAC) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base A (13/12);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 6.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCA  >  CP000730/1744229‑1744498
                                                                                                                                     |                                                                                                                                        
attATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGtt                                                                                                                                   <  1:201392/141‑1 (MQ=255)
 ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTa                                                                                                                                  >  1:133264/1‑141 (MQ=255)
     aaaTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAAc                                                                                                                              <  1:97708/141‑1 (MQ=255)
             tACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAc                                                                                                                      >  2:245813/1‑141 (MQ=255)
              aCCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAcc                                                                                                                     <  2:241674/141‑1 (MQ=255)
                ctactaCGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTCCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCat                                                                                                                   <  2:329716/141‑1 (MQ=255)
                  actacGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                                                                 >  2:106726/1‑141 (MQ=255)
                                tAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTc                                                                                                   <  1:224334/141‑1 (MQ=255)
                                     tttCAACTAATTTTATTGCTGATTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAAc                                                                                              >  1:244535/1‑141 (MQ=255)
                                      ttCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACg                                                                                             >  2:23606/1‑141 (MQ=255)
                                         aaCTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAAt                                                                                          <  2:167334/141‑1 (MQ=255)
                                           cTAATTTTATTGCTGCTTCAATCGTTCCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGa                                                                                        <  1:134226/141‑1 (MQ=255)
                                                                    accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCttt                                                               >  2:16776/1‑141 (MQ=255)
                                                                        ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTa                                                           >  2:281820/1‑141 (MQ=255)
                                                                        ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTTAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCTCTTCTTTTCTa                                                           >  1:284048/1‑141 (MQ=255)
                                                                                  ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaa                                                 >  2:278779/1‑141 (MQ=255)
                                                                                      atcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGcc                                             <  1:99500/141‑1 (MQ=255)
                                                                                        catcTGTAATTAACACACGTTGACCTGGGTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTTACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCaa                                           >  1:246780/1‑141 (MQ=255)
                                                                                          tcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATc                                         >  1:94935/1‑141 (MQ=255)
                                                                                          tcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATc                                         >  1:95006/1‑141 (MQ=255)
                                                                                          tcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATc                                         >  1:209835/1‑141 (MQ=255)
                                                                                                    acacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAAt                               <  2:24769/141‑1 (MQ=255)
                                                                                                     gaaataatGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATa                              <  2:34860/134‑1 (MQ=255)
                                                                                                     cacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATa                              >  2:318628/1‑141 (MQ=255)
                                                                                                       cacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCACTTCTTTTCTAACAAGTGCAAAGCCAATCCACATTATATaat                            >  2:279700/1‑140 (MQ=255)
                                                                                                               ccTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGa                    >  1:213710/1‑141 (MQ=255)
                                                                                                                   gTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGc                <  1:165918/141‑1 (MQ=255)
                                                                                                                        gctGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATg           <  2:58514/139‑1 (MQ=255)
                                                                                                                           gCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCAAACAGGTGCAAAGCCAACCCCCATTGAATAAGCTACAGGACAGCCAATGATa        >  1:78090/1‑141 (MQ=255)
                                                                                                                                 ttGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCa  >  1:39341/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                        
ATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCA  >  CP000730/1744229‑1744498

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