Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,026,424	C→T	10.2%	A220V (GCG→GTG)	guaD →	guanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,424	0	C	T	10.2%	139.1 / 6.0	49	A220V (GCG→GTG)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base T (2/3);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGT  >  NC_000913/3026280‑3026572                                                                                                                                                 |                                                                                                                                                     gCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAgatggcgatgg                                                                                                                                                 >  1:572722/1‑150 (MQ=255)    aCGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGgcgc                                                                                                                                              <  1:814232/150‑1 (MQ=255)          cGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACAAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGTGATGGCGCAACGcc                                                                                                                                        <  2:384502/150‑1 (MQ=255)           ggATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGccg                                                                                                                                       <  2:615486/150‑2 (MQ=255)            gATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTg                                                                                                                                      <  1:842499/150‑1 (MQ=255)              ttATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGTGATGGCGCAACGCCTGaa                                                                                                                                    >  1:1323521/1‑150 (MQ=255)              ttATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGTGATGGCGCAACGCCTGaa                                                                                                                                    >  1:1328325/1‑150 (MQ=255)                   tGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAaagaag                                                                                                                               >  1:380241/1‑150 (MQ=255)                    gCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAagaaga                                                                                                                              >  1:1362357/1‑150 (MQ=255)                     cTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGaaa                                                                                                                                   <  1:1056027/144‑1 (MQ=255)                     cTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGaaa                                                                                                                                   >  2:1056027/1‑144 (MQ=255)                         acacTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATc                                                                                                                         >  1:1251821/1‑150 (MQ=255)                         acacTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATc                                                                                                                         <  2:1362357/150‑1 (MQ=255)                           acTGCCGAAAGCAGCTATCACCATAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCg                                                                                                                       <  2:1045847/150‑1 (MQ=255)                           acTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCg                                                                                                                       >  1:528616/1‑150 (MQ=255)                                 gAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACg                                                                                                                 >  2:1176206/1‑150 (MQ=255)                                 gAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACg                                                                                                                 <  2:510139/150‑1 (MQ=255)                                    agcagcTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTgg                                                                                                              <  1:802689/150‑1 (MQ=255)                                           aTCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATAGGTGGGTACATa                                                                                                       <  2:1045835/150‑1 (MQ=255)                                             cACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGTGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATAcc                                                                                                     <  2:1328325/150‑1 (MQ=255)                                             cACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGTGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATAcc                                                                                                     <  2:1323521/150‑1 (MQ=255)                                                     caaagAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCtg                                                                                             >  2:1102305/1‑150 (MQ=255)                                                      aaagAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATACGGATACGTGGGTACATACCCATCTCtgt                                                                                            >  2:776685/1‑150 (MQ=255)                                                      aaagAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGATATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCtgt                                                                                            >  1:334348/1‑150 (MQ=255)                                                        agAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGa                                                                                          <  1:733041/150‑1 (MQ=255)                                                                          gCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGc                                                                        <  1:776685/150‑1 (MQ=255)                                                                           cacaAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGcc                                                                       >  2:257827/1‑150 (MQ=255)                                                                           cacaAAAATGGTAGACAGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCATGAACAGATGGCGATGGCGAAACGCCTGAAAGAAGAATATCCGGATAAGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGcc                                                                       >  2:397686/1‑150 (MQ=255)                                                                             caAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTg                                                                     >  1:520637/1‑150 (MQ=255)                                                                                   tGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGaa                                                               >  2:649854/1‑150 (MQ=255)                                                                                     gtcgtcTGCTATATTCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAAt                                                             >  2:1420580/1‑150 (MQ=255)                                                                                               atatGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATcc                                                   <  1:1338306/150‑1 (MQ=255)                                                                                               atatGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATcc                                                   <  2:528616/150‑1 (MQ=255)                                                                                                tatGCGATTACGCCACGCTTCGCCCCGACCTCGTCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCt                                                  >  1:985359/1‑150 (MQ=255)                                                                                                   gCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGAc                                               <  1:1176206/150‑1 (MQ=255)                                                                                                    cGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGatc                                              <  1:620097/150‑3 (MQ=255)                                                                                                    cGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGAcc                                              >  2:1310841/1‑150 (MQ=255)                                                                                                     gATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCa                                             <  2:572722/150‑1 (MQ=255)                                                                                                       ttACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCatg                                           <  2:520637/150‑1 (MQ=255)                                                                                                             cACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGGCCAACGCCTGAAAGAAGAAGATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGtt                                     >  2:1194716/1‑150 (MQ=255)                                                                                                               cGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTAt                                   >  1:168154/1‑150 (MQ=255)                                                                                                                   tCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGATTTATCCTGTCCATGATTGTTATCTgg                               <  2:1348000/150‑1 (MQ=255)                                                                                                                            ccTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTAcc                      >  1:842235/1‑150 (MQ=255)                                                                                                                              tCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGt                           <  2:1174074/143‑1 (MQ=255)                                                                                                                              tCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGt                           >  1:1174074/1‑143 (MQ=255)                                                                                                                                  ctcCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGt                >  2:95272/1‑150 (MQ=255)                                                                                                                                    ccTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTAc              >  2:1295552/1‑150 (MQ=255)                                                                                                                                    ccTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTAc              >  1:56781/1‑150 (MQ=255)                                                                                                                                     cTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACg             >  2:1003837/1‑150 (MQ=255)                                                                                                                                     cTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACg             <  1:1310841/150‑1 (MQ=255)                                                                                                                                               ggcgatggcgCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCgg   <  1:1420580/150‑1 (MQ=255)                                                                                                                                                gcgatggcgCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGt  >  1:1326184/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACAGATGGCGATGGCGCAACGCCTGAAAGAAGAATATCCGGATACGTGGGTACATACCCATCTCTGTGAAAACAAAGATGAAATTGCCTGGGTGAAATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGT  >  NC_000913/3026280‑3026572

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