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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,744,3620GA97.5% 114.2 / ‑3.7 41H104Y (CAC→TAC) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/23);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

ATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCA  >  CP000730/1744229‑1744498
                                                                                                                                     |                                                                                                                                        
attATAACTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGtt                                                                                                                                   <  1:706859/141‑1 (MQ=255)
 ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTTTACATTGTTa                                                                                                                                  >  1:65935/1‑141 (MQ=255)
 ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTa                                                                                                                                  <  1:154677/141‑1 (MQ=255)
     aaaTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAAc                                                                                                                              <  1:279734/141‑1 (MQ=255)
           aaTACCCACTACGATCCCGCCTAATTTTCCAACTAATTTTATTGCTGCCTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGt                                                                                                                        <  1:537415/141‑1 (MQ=255)
            aTACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTa                                                                                                                       >  1:241769/1‑141 (MQ=255)
               ccTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAcct                                                                                                                    >  1:577005/1‑140 (MQ=255)
                  actacGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                                                                 >  1:49761/1‑141 (MQ=255)
                            cGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACAACCGGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCTTTTTCTAGGTCAt                                                                                                       >  1:223217/1‑141 (MQ=255)
                                tAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTTTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGATCATACTc                                                                                                   >  1:112742/1‑141 (MQ=255)
                                   tttttCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCCTCTGGGCGTAACTTACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATa                                                                                                <  1:912023/141‑1 (MQ=255)
                                     tttCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAAc                                                                                              >  1:382944/1‑141 (MQ=255)
                                                tttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATGTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCa                                                                                   <  1:661793/141‑1 (MQ=255)
                                                                 cGTACCACCAGTAGATAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAATATATATTTACAAGGGTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttc                                                                  >  1:484847/1‑141 (MQ=255)
                                                                     ccaccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt                                                              <  1:109078/141‑1 (MQ=255)
                                                                     ccaccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt                                                              <  1:851262/141‑1 (MQ=255)
                                                                                 aataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCa                                                  <  1:612688/141‑1 (MQ=255)
                                                                                    aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAg                                               <  1:845460/141‑1 (MQ=255)
                                                                                       tcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCa                                            <  1:523519/141‑1 (MQ=255)
                                                                                         atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAAt                                          >  1:549395/1‑141 (MQ=255)
                                                                                         atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAAt                                          <  1:140220/141‑1 (MQ=255)
                                                                                         atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAAt                                          >  1:267198/1‑141 (MQ=255)
                                                                                              tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCTATTGACTGCTAACAGGTGCAAAGCCAATCCCCa                                     <  1:919700/141‑1 (MQ=255)
                                                                                              tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCa                                     >  1:113628/1‑141 (MQ=255)
                                                                                                 ttAACACACGTTGACCTGGTTTAATTGCATCTTTTTACATAGATCCCCCGTTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTg                                  <  1:436949/141‑1 (MQ=255)
                                                                                                  tAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGa                                 <  1:824606/141‑1 (MQ=255)
                                                                                                  tAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGa                                 <  1:722808/141‑1 (MQ=255)
                                                                                                    acatacGTTGACCTGGTTTGATTGCATGTTTGTACATTGTGAAAACATTTGGACCATATTCGAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAAt                               <  1:879943/141‑1 (MQ=255)
                                                                                                    acacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAAt                               <  1:612846/141‑1 (MQ=255)
                                                                                                    acacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAAt                               <  1:285947/141‑1 (MQ=255)
                                                                                                          ggtGACCTGGTTCAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGTTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTa                         <  1:805713/139‑1 (MQ=255)
                                                                                                           ttGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAc                        <  1:65096/141‑1 (MQ=255)
                                                                                                           gtcacgtGGTTTCATTGCATCTTTGTCCATTGTTAACACATTTTCACCATCTTCTAGGTCATCCTCATAACGAATGACTTCCCGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAc                        <  1:206661/135‑1 (MQ=255)
                                                                                                            tGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACa                       >  1:193297/1‑141 (MQ=255)
                                                                                                                cTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACCGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGAc                   >  1:741000/1‑141 (MQ=255)
                                                                                                                cTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGAc                   >  1:833240/1‑141 (MQ=255)
                                                                                                                cTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGAc                   >  1:740986/1‑141 (MQ=255)
                                                                                                                    tttAATCGCATCCTTGCACATTGTTCCAACATTTGTACCATATTCTAGGTCCTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGcc               <  1:192576/141‑1 (MQ=255)
                                                                                                                       aaTTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAt            <  1:806811/141‑1 (MQ=255)
                                                                                                                            cATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACCAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCTCATTGAATAAGCTACAGGACAGCCAATGATaa       >  1:618896/1‑141 (MQ=255)
                                                                                                                                 ttGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCa  <  1:376517/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                        
ATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCA  >  CP000730/1744229‑1744498

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