Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,168,106	A→G	7.4%	N170D (AAC→GAC)	ycfJ →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,106	0	A	G	7.4%	155.2 / 3.1	54	N170D (AAC→GAC)	ycfJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base G (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

ACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGT  >  NC_000913/1167959‑1168254                                                                                                                                                    |                                                                                                                                                     aCTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCaa                                                                                                                                                     >  7:347771/1‑149 (MQ=255)   ttACACGACTACGCAACAGAGTTCTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGAACTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACg                                                                                                                                                   >  4:359367/1‑149 (MQ=255)      cacGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATACGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACggg                                                                                                                                                 <  4:36917/148‑1 (MQ=255)         gACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCACAATCCGCATGGACCGCGATCCGGGTACGCTGATGCCGCTTGATAGCAACGGGCAAc                                                                                                                                             >  6:239906/1‑149 (MQ=255)         gACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAAc                                                                                                                                             <  4:157799/149‑1 (MQ=255)         gACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAAc                                                                                                                                             >  8:248002/1‑149 (MQ=255)                           tAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAg                                                                                                                           >  1:90512/1‑149 (MQ=255)                              aaCGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGtat                                                                                                                        <  8:347771/149‑1 (MQ=255)                                     tATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAgg                                                                                                                 <  7:86915/149‑1 (MQ=255)                                           aaGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTAc                                                                                                           >  1:391608/1‑149 (MQ=255)                                            aGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGCTCCAGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTAct                                                                                                          >  3:54103/1‑149 (MQ=255)                                             gTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACggg                                                                                                                                                 <  8:202490/109‑1 (MQ=255)                                             gTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTActc                                                                                                         >  8:244851/1‑149 (MQ=255)                                             gTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGATAGATAGCAACggg                                                                                                                                                 >  7:202490/1‑109 (MQ=255)                                                aGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGAATGGACCGCGATCCGGTTACGCAGATCCCGCTAGATAGCGACGGGCAACTGATTTTGAATAACAAACTATAACAAGGCTGTACTCTg                                                                                                       >  5:304002/1‑148 (MQ=255)                                                    aaaaTGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAAtt                                                                                                  <  8:219533/149‑1 (MQ=255)                                                      aaTGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACTGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTg                                                                                                >  1:68780/1‑149 (MQ=255)                                                        tGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGc                                                                                              <  1:3502/149‑1 (MQ=255)                                                           tCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGcccc                                                                                           >  6:243215/1‑149 (MQ=255)                                                              gTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCa                                                                                        >  2:32188/1‑149 (MQ=255)                                                               ttATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCAt                                                                                       <  1:32188/149‑1 (MQ=255)                                                                  tgatgTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTc                                                                                     <  4:301233/148‑1 (MQ=255)                                                                            tataAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCt                                                                            >  4:100909/1‑147 (MQ=255)                                                                             ataAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAg                                                                         >  4:279737/1‑149 (MQ=255)                                                                              taAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGa                                                                          <  2:68780/147‑1 (MQ=255)                                                                              taAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAgg                                                                        <  2:90512/149‑1 (MQ=255)                                                                                 gATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGc                                                                     <  2:391608/149‑1 (MQ=255)                                                                                     ggCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCttt                                                                  <  3:359367/148‑1 (MQ=255)                                                                                        gATCATCAGTGAAAAATCCGCATGGACCGCGATCCGTGTACGCAGATCCCGCTAGATAGCGACTGGCAACTGATTTTGAATAACAAAGTATAACTAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGTGGCttttttt                                                              <  6:304002/149‑1 (MQ=255)                                                                                                gggCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCGACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCttt                                                                  >  7:252427/1‑137 (MQ=255)                                                                                                gggCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCGACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCttt                                                                  <  8:252427/137‑1 (MQ=255)                                                                                                gggCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTa                                                      >  5:256780/1‑149 (MQ=255)                                                                                                   cAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTAttt                                                   >  8:137238/1‑149 (MQ=255)                                                                                                    aaaaTCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGACTAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTc                                                  <  8:35892/149‑1 (MQ=255)                                                                                                    aaaaTCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTc                                                  <  7:244851/149‑1 (MQ=255)                                                                                                          cGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGt                                            >  4:2467/1‑149 (MQ=255)                                                                                                           gCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGtt                                           <  3:100909/149‑1 (MQ=255)                                                                                                               ggACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTAttt                                                   >  8:319990/1‑137 (MQ=255)                                                                                                               ggACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTAttt                                                   <  7:319990/137‑1 (MQ=255)                                                                                                                  ccGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAACGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGCCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCa                                    >  2:151870/1‑149 (MQ=255)                                                                                                                   cgcgATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCa                                                                                <  4:309298/104‑1 (MQ=255)                                                                                                                     cgATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGcc                                                                                             <  5:57618/89‑1 (MQ=255)                                                                                                                     cgATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGcc                                                                                             >  6:57618/1‑89 (MQ=255)                                                                                                                     cgATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATaa                                 <  3:279737/149‑1 (MQ=255)                                                                                                                      gATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTAttt                                                   >  3:357159/1‑130 (MQ=255)                                                                                                                      gATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTAttt                                                   <  4:357159/130‑1 (MQ=255)                                                                                                                                  cAGATCCCGCTATATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGtt                    <  6:356742/149‑1 (MQ=255)                                                                                                                                  cAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGtt                    <  1:53784/149‑1 (MQ=255)                                                                                                                                   aGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTc                   >  4:375636/1‑149 (MQ=255)                                                                                                                                   aGATCCCGCTAGATAGCAACGGGCAACTGATATTGAATAACAAAGTATAACAAGGCTGAACTCTGCAATTTGGCCCCTCAATCGCTCAGGCTGAGGGGCTTATATTGCGACTAATTTCACCAGATCGGGCCATAAACGCAAAGTCGTTc                   >  3:277188/1‑149 (MQ=255)                                                                                                                                           cTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATc           <  5:299145/149‑1 (MQ=255)                                                                                                                                            tAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCt          >  6:138616/1‑149 (MQ=255)                                                                                                                                              gATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGc        <  5:138616/149‑1 (MQ=255)                                                                                                                                               aTAGCAACGGGCAACTGATTTTGAATAACAAAGTATAAGAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGCCTTATTTCACCAGTTCGGGCCATACACGCAAAGTGGTTCCGGCAATCTGCa       >  2:120308/1‑149 (MQ=255)                                                                                                                                                  gCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATc     <  7:137238/148‑1 (MQ=255)                                                                                                                                                    aaCGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGt  <  5:243215/149‑1 (MQ=255)                                                                                                                                                    aaCGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGt  >  5:276472/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGT  >  NC_000913/1167959‑1168254

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