Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,269,459	G→A	A16T (GCT→ACT)	aphA →	acid phosphatase/phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,269,459	0	G	A	96.9%	96.9 / ‑6.3	33	A16T (GCT→ACT)	aphA	acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/15);  total (16/17)
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

CTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTA  >  NC_000913/4269215‑4269697                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               cTTACATCGAGTAATAAAAATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTNTa                                                                                                                                                                                                                                           >  2:592592/1‑250 (MQ=255)    aCATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAAc                                                                                                                                                                                                                                        >  1:38317/1‑250 (MQ=255)                        tttatagcaaaaaaGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTCTGGTCAGATGGTTTTTGCAACAAATCTCACAATAACAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGcc                                                                                                                                                                                                                   >  2:349931/9‑251 (MQ=255)                                                         ggTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAAc                                                                                                                                                                                  <  1:365403/251‑1 (MQ=255)                                                                                      acaAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCcat                                                                                                                                                     >  2:10586/1‑251 (MQ=255)                                                                                        aaaaaCCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAAGCTCACAATAAAAAGTGTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTGCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCcattc                                                                                                                                                   <  1:106867/251‑1 (MQ=255)                                                                                          aaaCCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCcattcat                                                                                                                                                 <  2:38317/251‑1 (MQ=255)                                                                                                      aTCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAATCAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCCCCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTGTCGGTc                                                                                                                                     <  2:415830/251‑1 (MQ=255)                                                                                                           tccttTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGcac                                                                                                                                 <  1:566103/250‑1 (MQ=255)                                                                                                                        aGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGc                                                                                                                   <  1:10586/251‑1 (MQ=255)                                                                                                                          cTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAACTAGcat                                                                                                                 >  2:575194/1‑249 (MQ=255)                                                                                                                                  aTGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCGCCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAggg                                                                                                          <  2:510382/250‑1 (MQ=255)                                                                                                                                                  tGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCgg                                                                                          <  2:17649/250‑1 (MQ=255)                                                                                                                                                              cacaATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGat                                                                               <  2:10460/249‑1 (MQ=255)                                                                                                                                                                       aaGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGCTGacac                                                                    >  2:112099/1‑251 (MQ=255)                                                                                                                                                                       aaGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGaaac                                                                    >  1:518289/1‑248 (MQ=255)                                                                                                                                                                         gTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACgg                                                                  >  1:248543/1‑251 (MQ=255)                                                                                                                                                                           ttCAACATACTTACTATGTAGGGAAAAATATGCGCATGATCACACAGGTAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTGAACCCTGGGACTAACGTTTCCAGGCTTGCTTAACAGGCCCCCATTCATTGGGTTTCGGTCTCACAAATTGAAAATAGCCTCTCAGGGCGTCCGCCAATGGCGTTGTGGTTTGATATCGATGACACGGt                                                                 <  2:298252/250‑1 (MQ=255)                                                                                                                                                                           ttCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACggg                                                                 >  2:447532/1‑249 (MQ=255)                                                                                                                                                                               aCATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGGCACGGTACttt                                                            >  1:562582/1‑251 (MQ=255)                                                                                                                                                                                          aTTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCg                                                 <  2:518289/251‑1 (MQ=255)                                                                                                                                                                                           tttAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCgg                                                >  1:552714/1‑251 (MQ=255)                                                                                                                                                                                                ggAAAAATATGCGCCAGATCCCCCAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAGCAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCCCCGCTTAACCCTGGGACTACCGTTGCCAGGCTTGCTGAACATGCCCCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATCGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTGGATATCGATGACACGGTACTTTGTTCCAGTCCGGGCTTc                                           <  1:519877/251‑1 (MQ=255)                                                                                                                                                                                                   aaaaTATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCCCCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTgg                                        <  2:447615/251‑1 (MQ=255)                                                                                                                                                                                                   aaaaTATGCCCAAGATCACACAGGCAATCATTTCCGTTTTCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTgg                                        <  2:140111/251‑1 (MQ=255)                                                                                                                                                                                                        atGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACt                                                              <  2:381600/224‑1 (MQ=255)                                                                                                                                                                                                             cAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCaaa                               >  1:301896/1‑250 (MQ=255)                                                                                                                                                                                                                 aTCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCCCCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAc                          <  2:58857/251‑1 (MQ=255)                                                                                                                                                                                                                       caGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATTTCGATGACACGGTACTTTTTTTCAGTCCCGGCTTCTGGCGCGGCAAAAAAACCTtctc                    >  2:614016/1‑251 (MQ=255)                                                                                                                                                                                                                        aGGCAATCAGTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGAGGGGGTTGATATCGGTGACACGGGACTTTTTTCCAGTCCGGGCTTTTGGCGCGGCAAAAAAACCTTTTCg                   >  2:142793/1‑251 (MQ=255)                                                                                                                                                                                                                                 gTGCCGTTTGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTTATATCGATTAAACGGGACATTTTTTCAGGCCCGGCTTTTGGCGCAGGACAAAAACCtttctcacacaactc          >  2:100946/1‑238 (MQ=255)                                                                                                                                                                                                                                        ttGCTTATTGTTCACTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGGTGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGAtt   >  1:424335/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCTTATTGTTCTCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTa  <  1:349931/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               CTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTA  >  NC_000913/4269215‑4269697

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