Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,844,413	C→T	S150S (AGC→AGT)	adeD →	adenine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,844,413	0	C	T	96.8%	99.0 / NA	31	S150S (AGC→AGT)	adeD	adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/15);  minor base A (1/0);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAT  >  NC_000913/3844188‑3844655                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   ggATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCt                                                                                                                                                                                                                           <  2:170020/250‑1 (MQ=255)                aaCGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTg                                                                                                                                                                                                            <  2:194608/251‑1 (MQ=255)                  cggcggTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGcc                                                                                                                                                                                                          <  2:47222/251‑1 (MQ=255)                    gaggTGCCCGGGTTTATTGATGCTCGCCTGCATATTGAATCCAGCATGATGACGAAGGTCACTTTTGACACCGCTACAATGCCGCGCGGCATGACGACCGTTATTTGCGACAATCATGAAGTCGTCAACGTGATGGTCGAAGCAGGATTCGCCTAGTTTGCCCGATGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTAGGTTCTTGCGTACTCGCCCTGGCAGGCTGAGATGCTAACGGttcc                                                                                                                                                                                                          >  2:370913/3‑246 (MQ=255)                                   aTTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAAc                                                                                                                                                                                          <  2:2914/250‑1 (MQ=255)                                             cccttcATATTGAATCCAGCATGATGACGCCGGTCCCTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATTAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGcc                                                                                                                                                                               <  2:177124/246‑1 (MQ=255)                                                         aatccagaaggatgcagCCGGTAACTTTTTAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACGCTCATGAATTCGTCAACGTGATGGGCGAGGCCGGCTTCGCCTGATTTGCCCGCTGTGCCGAACAGGCACGGCAAAACCAGTACTTACAGGTCAGATCTTGCGTACCCGCCCTGGCAGGCTGCGATGCTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCac                                                                                                                                                                   >  2:250908/17‑250 (MQ=255)                                                                     tgaagccagtaccTTTTGAAACCGCTACACTGCCGAGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGAAGGGCGAAGCCGGATTCGCCTGGTTAGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTAAGTTCTTGCGTACCCGACCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGGACAGCGGCTCGACTGGCGGGGCCATCCGCAGGTTAcc                                                                                                                                                       >  2:30036/13‑251 (MQ=255)                                                                                    ttGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTGGCCTGGTTTGCCCGCTGTGCCGTACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTTCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAatg                                                                                                                                        >  1:134452/1‑251 (MQ=255)                                                                                       aaaCCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAatgatg                                                                                                                                     <  1:23751/251‑1 (MQ=255)                                                                                               aNCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTAc                                                                                                                               <  2:357208/249‑1 (MQ=255)                                                                                               aCCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTAcc                                                                                                                              <  1:323688/250‑1 (MQ=255)                                                                                                               tGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCggg                                                                                                             >  1:171841/1‑251 (MQ=255)                                                                                                                  cgacCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAg                                                                                                          <  1:33807/251‑1 (MQ=255)                                                                                                                     ccGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAAt                                                                                                       <  2:217194/251‑1 (MQ=255)                                                                                                                      cGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCATGGTTTGCCCGCTGTGCCGAACAGGAAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATg                                                                                                      >  1:161165/1‑251 (MQ=255)                                                                                                                       gTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGGTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATTATGGACTACCCTGGCGTAATTAGCGGGCAGAATgc                                                                                                     >  2:73869/1‑251 (MQ=255)                                                                                                                           tttGCGACCCTCATGAAAGCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCTCCTGGCGTGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCTGGCTGAATGCgct                                                                                                  <  2:284808/250‑1 (MQ=255)                                                                                                                              gCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATGCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCgctgct                                                                                               <  2:203183/250‑1 (MQ=255)                                                                                                                                                 cAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTc                                                                           >  1:133637/1‑251 (MQ=255)                                                                                                                                                             cGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGc                                                               >  2:302183/1‑251 (MQ=255)                                                                                                                                                              gAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCt                                                              <  1:271944/251‑1 (MQ=255)                                                                                                                                                                                 ttGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCggg                                            <  2:136406/250‑1 (MQ=255)                                                                                                                                                                                         tgtgCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTTCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCggg                                            <  2:178762/242‑1 (MQ=255)                                                                                                                                                                                         tgtgCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCggg                                            >  1:178762/1‑242 (MQ=255)                                                                                                                                                                                           tgCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTAACGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGTCAGAATGCGCTGCTCGATGAACTGAATGCATTTCGCCACCTGACGCTGGACGGTCACTGCACGGGTTTGgggggt                                  >  1:329312/1‑250 (MQ=255)                                                                                                                                                                                            gCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGAATTTCGCAACCTGACGCTGGACGGTCACTGCCCGGGTTTGggtggca                                 >  1:197847/1‑248 (MQ=255)                                                                                                                                                                                                           caaACCCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCCCTGCCCGGGTTTGGGGGGTAAAGAACTTAACGCCt                  >  1:292165/1‑250 (MQ=255)                                                                                                                                                                                                            aaaaCCAGTACTTACAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCtat                >  1:295840/1‑251 (MQ=255)                                                                                                                                                                                                                     acttaacggTCAGTTCTTGCGTACCCGCGCTGGAAGACTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCAGCAGATTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGGAATTAGCGGGGAGAATGCGATGCTCGATAACATGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGGGGGAAAGAACATAAAGCCTATATTACTGCg       >  1:113869/8‑251 (MQ=255)                                                                                                                                                                                                                          cAGGTCAGTTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAt  >  1:243759/1‑251 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   TGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAT  >  NC_000913/3844188‑3844655

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