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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,744,362	0	G	A	86.8%	135.2 / 12.4	56	H104Y (CAC→TAC)	apt	adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base A (21/27);  total (24/32)
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.66e-01

AATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCG  >  CP000730/1744228‑1744502                                                                                                                                        |                                                                                                                                             aattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGt                                                                                                                                        >  1:268908/1‑141 (MQ=255) aattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGt                                                                                                                                        <  1:609352/141‑1 (MQ=255) aattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGt                                                                                                                                        <  2:159973/141‑1 (MQ=255)  attATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGtt                                                                                                                                       >  2:93673/1‑141 (MQ=255)   ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTa                                                                                                                                      >  1:624932/1‑141 (MQ=255)    tataAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTaa                                                                                                                                     <  2:179078/141‑1 (MQ=255)      taAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTaaaa                                                                                                                                   <  1:371726/141‑1 (MQ=255)       aaaTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAAc                                                                                                                                  <  1:30578/141‑1 (MQ=255)       aaaTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAAc                                                                                                                                  <  2:463829/141‑1 (MQ=255)       aaaTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAAc                                                                                                                                  <  1:485274/141‑1 (MQ=255)           gCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACAttt                                                                                                                              >  1:696641/1‑141 (MQ=255)             aaTACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGt                                                                                                                            >  1:232108/1‑141 (MQ=255)                 ccTACTAGTATGTCGCCTAATTTTTCAACTAATTGTATTGCTGCTGCAGGCGTACCTCCGTTATCTAATAAATCGTCTGTAATTAACACTCGTTGACCTTGTTTCATTGCA‑TCTTTGTACATTGTTTTAACATTTGTACCa                                                                                                                        <  2:582821/141‑1 (MQ=255)                      tacGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTACCATATTc                                                                                                                   <  2:657342/141‑1 (MQ=255)                      tacGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTc                                                                                                                   <  2:647254/141‑1 (MQ=255)                      tacGATACCGCCTAATATTTCAACGAATTTTTTTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTTCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTc                                                                                                                   <  2:392163/141‑1 (MQ=255)                       acGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCt                                                                                                                  >  1:542539/1‑141 (MQ=255)                                      ttttCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATaa                                                                                                   <  1:57886/141‑1 (MQ=255)                                         tCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGa                                                                                                <  2:33238/141‑1 (MQ=255)                                          cAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGaa                                                                                               <  2:764295/141‑1 (MQ=255)                                            aCTAATTTTATTGCTGCGTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAAGTAACACACGTTGACCGGAGTTAATGTCA‑TCGGTTGACATTTTTTAAAGAGTTGTACCATATTCTAGGTCATACTCATAACGAAtt                                                                                             >  1:393866/1‑140 (MQ=255)                                                     aTTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGa                                                                                    <  2:424245/141‑1 (MQ=255)                                                       tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAgg                                                                                  <  2:731242/141‑1 (MQ=255)                                                       tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAgg                                                                                  <  1:479831/141‑1 (MQ=255)                                                         ctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTAcc                                                                                                                         >  1:498335/1‑100 (MQ=255)                                                         ctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTAcc                                                                                                                         <  2:498335/100‑1 (MQ=255)                                                            ctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAtt                                                                             <  1:215555/141‑1 (MQ=255)                                                               cAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                                                                          <  1:42346/141‑1 (MQ=255)                                                                caTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGNCATNCTCATAACGAATGACTTCACGAGGTAATTTccc                                                                         <  1:633511/140‑1 (MQ=255)                                                                aaTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTccc                                                                         <  1:93673/141‑1 (MQ=255)                                                                aaTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTccc                                                                         >  2:179558/1‑141 (MQ=255)                                                                  tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATATCGTATTACTTCACGAGGTAATTTCCctt                                                                       >  2:231871/1‑141 (MQ=255)                                                                  tCGTACCACCAGTAGCTAATAAATCATCTGTAATCAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCctt                                                                       <  1:726244/141‑1 (MQ=255)                                                                   cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttc                                                                      >  1:492608/1‑141 (MQ=255)                                                                     tACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCcttctt                                                                    >  1:478568/1‑141 (MQ=255)                                                                      accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCttt                                                                   >  2:31751/1‑141 (MQ=255)                                                                       ccaccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt                                                                  <  2:721606/141‑1 (MQ=255)                                                                        caccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTc                                                                 >  2:470395/1‑141 (MQ=255)                                                                         accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCt                                                                >  1:112307/1‑141 (MQ=255)                                                                          ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTa                                                               <  1:525314/141‑1 (MQ=255)                                                                              tAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg                                                           <  2:269702/141‑1 (MQ=255)                                                                               aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTCATTTCCCTTCTTTTCTAACAgg                                                          <  2:313483/141‑1 (MQ=255)                                                                                 cTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTg                                                        >  2:483289/1‑141 (MQ=255)                                                                                    ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACACATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaa                                                     >  1:637364/1‑141 (MQ=255)                                                                                         tcatcTGTAATTAACACACGTTGATCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCa                                                <  1:232912/141‑1 (MQ=255)                                                                                         tcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAACCATTTGTACCATATTCTATGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGTTGCAAAGCCa                                                >  1:267211/1‑141 (MQ=255)                                                                                         tcatcTGTAATTAACACACGTCGACCTCGTGTAATTCCACGCGTTGTTGGGT‑TTAAGATTGTGGGACCTGATGCTGGGGCTGACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCGTTTCTAACAGGTGCAAAGCCa                                                <  1:704003/141‑1 (MQ=255)                                                                                          catcTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCaa                                               >  1:442950/1‑141 (MQ=255)                                                                                           atcTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAAt                                              >  1:603761/1‑141 (MQ=255)                                                                                           atcTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAAt                                              <  2:739512/141‑1 (MQ=255)                                                                                            tcTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATc                                             >  2:129684/1‑141 (MQ=255)                                                                                                tAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCa                                         <  1:511448/141‑1 (MQ=255)                                                                                                           cGTTGACCTGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCt                              >  2:444620/1‑141 (MQ=255)                                                                                                                   tGGTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACa                      >  1:514864/1‑141 (MQ=255)                                                                                                                     gTTTAATTGCA‑TCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGc                    <  1:475307/141‑1 (MQ=255)                                                                                                                                 tCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAg         >  1:461582/1‑141 (MQ=255)                                                                                                                                        gCATTGTNAAAACATTTGTAACACACTCTAGGTCATACTCATAAAGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAGAAGGTGCAAAGCCCCTCCCCATTGAATAAGTTATTGGCCAGCCGATGATAAAGCCAcgcg  <  2:589164/141‑1 (MQ=255)                                                                                                                                        |                                                                                                                                             AATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCA‑TCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCG  >  CP000730/1744228‑1744502

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