Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,702,865	T→G	I211M (ATT→ATG)	add →	adenosine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,702,865	0	T	G	98.2%	193.8 / ‑4.5	55	I211M (ATT→ATG)	add	adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (24/30);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAGCTGGCTG  >  NC_000913/1702720‑1703010                                                                                                                                                  |                                                                                                                                                  ttACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTg                                                                                                                                               >  2:18248/1‑150 (MQ=255)    ccGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACt                                                                                                                                           >  2:309206/1‑151 (MQ=255)    ccGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACt                                                                                                                                           <  2:771203/151‑1 (MQ=255)       cACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTggg                                                                                                                                        <  2:433794/151‑1 (MQ=255)        aCTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTggg                                                                                                                                        >  1:208245/1‑150 (MQ=255)        aCTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGt                                                                                                                                       >  2:424636/1‑151 (MQ=255)         cTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTg                                                                                                                                      >  1:346568/1‑151 (MQ=255)           tGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTg                                                                                                                                      <  2:418412/149‑1 (MQ=255)             aTTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGa                                                                                                                                  <  1:620599/151‑1 (MQ=255)              tttAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCgg                                                                                                                                   >  2:377753/1‑149 (MQ=255)               ttAGCCGGTGATGAACTTGGTTTACCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAg                                                                                                                                 >  2:206401/1‑150 (MQ=255)                 aGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCg                                                                                                                               <  1:364801/150‑1 (MQ=255)                  gCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAATCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTa                                                                                                                             <  2:918168/151‑1 (MQ=255)                   ccGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTAt                                                                                                                            <  2:775338/151‑1 (MQ=255)                     ggTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTg                                                                                                                          >  1:819815/1‑151 (MQ=255)                      gTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTgg                                                                                                                         <  2:283919/151‑1 (MQ=255)                        gatgaACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGAc                                                                                                                       <  1:701870/151‑1 (MQ=255)                           gaACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATg                                                                                                                    <  1:774957/151‑1 (MQ=255)                             aCTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGc                                                                                                                  <  1:447046/151‑1 (MQ=255)                                tGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGt                                                                                                                >  2:248666/1‑150 (MQ=255)                                 ggTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTa                                                                                                               <  2:34509/150‑1 (MQ=255)                                  gTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTaaa                                                                                                             <  1:451516/151‑1 (MQ=255)                                       ccGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCa                                                                                                        <  2:531772/151‑1 (MQ=255)                                                  gTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGa                                                                                                    >  2:418480/1‑144 (MQ=255)                                                  gTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGAt                                                                                                <  2:185528/148‑1 (MQ=255)                                                        ttCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTg                                                                                       >  2:84695/1‑151 (MQ=255)                                                            cACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGa                                                                                      <  2:526460/148‑1 (MQ=255)                                                             aCTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGa                                                                                  <  2:8076/151‑1 (MQ=255)                                                                  aaCCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtc                                                                             >  1:849703/1‑151 (MQ=255)                                                                  aaCCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACCGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtc                                                                             >  2:530733/1‑151 (MQ=255)                                                                       cgcgcgTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGcc                                                                        <  2:592178/151‑1 (MQ=255)                                                                          gcgTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCg                                                                       >  2:591240/1‑149 (MQ=255)                                                                            gTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGc                                                                    >  1:499830/1‑150 (MQ=255)                                                                              gATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGcaa                                                                  <  2:912687/150‑1 (MQ=255)                                                                                  cGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAAt                                                             <  2:532117/151‑1 (MQ=255)                                                                                   gggCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAAtt                                                            >  1:119748/1‑151 (MQ=255)                                                                                           atatTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTg                                                     >  1:517218/1‑150 (MQ=255)                                                                                           atatTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGa                                                    <  2:133433/151‑1 (MQ=255)                                                                                             atTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAAt                                                  <  1:550484/151‑1 (MQ=255)                                                                                                aCCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATcc                                                <  1:281455/150‑1 (MQ=255)                                                                                                aCCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCt                                               <  2:272901/151‑1 (MQ=255)                                                                                                 ccGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTg                                              <  1:17722/151‑1 (MQ=255)                                                                                                   gTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTc                                            >  2:264432/1‑151 (MQ=255)                                                                                                     ccATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTg                                          <  2:346580/151‑1 (MQ=255)                                                                                                               gAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAat                                 <  2:477490/150‑1 (MQ=255)                                                                                                               gAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAat                                 >  2:223325/1‑150 (MQ=255)                                                                                                                aaGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAatat                               <  2:829092/151‑1 (MQ=255)                                                                                                                 aGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATAtt                              >  2:97138/1‑151 (MQ=255)                                                                                                                            cggAAAGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCAc                   >  2:66654/1‑151 (MQ=255)                                                                                                                               aaaGCATCTGGCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGt                >  1:157411/1‑151 (MQ=255)                                                                                                                                         gCAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctg      <  1:450038/151‑1 (MQ=255)                                                                                                                                          cAGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctgg     >  2:98534/1‑151 (MQ=255)                                                                                                                                           aGGCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctggc    >  2:490391/1‑151 (MQ=255)                                                                                                                                            ggCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctggct   >  2:313649/1‑151 (MQ=255)                                                                                                                                             gCGATGCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctggctg  <  2:276997/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAGCTGGCTG  >  NC_000913/1702720‑1703010

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