Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,633,497	G→T	S83R (AGC→AGA)	yjjX ←	ITPase/XTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,633,497	0	G	T	97.9%	167.8 / NA	47	S83R (AGC→AGA)	yjjX	ITPase/XTPase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (22/24);  total (22/25)
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

ACCGCAGGCAACGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACC  >  NC_000913/4633420‑4633572                                                                              |                                                                            aCCGCAGGCAACGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCt                                                                           <  1:939222/80‑1 (MQ=255)  ccGCAGGCAACGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGt                                                                         <  2:939223/81‑1 (MQ=255)   cGCAGGCAACGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGt                                                                         >  1:939288/1‑80 (MQ=255)      aGGCAACGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCAc                                                                      >  2:939286/1‑80 (MQ=255)         cAACGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATc                                                                  >  2:939285/1‑81 (MQ=255)          aaCGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATc                                                                  >  2:939287/1‑80 (MQ=255)            cGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGAt                                                               <  1:939224/81‑1 (MQ=255)             ggCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATg                                                              <  1:939225/81‑1 (MQ=255)             ggCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATg                                                              >  1:939280/1‑81 (MQ=255)                aaaGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGc                                                             >  1:939299/1‑79 (MQ=255)                aaaGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCa                                                           <  1:939226/81‑1 (MQ=255)                 aaGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAg                                                          <  2:939227/81‑1 (MQ=255)                 aaGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAg                                                          >  1:939320/1‑81 (MQ=255)                  aGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGc                                                         >  1:939312/1‑81 (MQ=255)                      gCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTc                                                      <  1:939228/80‑1 (MQ=255)                      gCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCa                                                     <  2:939229/81‑1 (MQ=255)                        cGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAAt                                                   >  2:939324/1‑81 (MQ=255)                        cGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAAt                                                   <  1:939230/81‑1 (MQ=255)                         gAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATc                                                  <  2:939231/81‑1 (MQ=255)                            cgcgCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGc                                                >  1:939283/1‑80 (MQ=255)                             gcgcCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCa                                              >  1:939305/1‑81 (MQ=255)                             gcgcCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCa                                              <  2:939232/81‑1 (MQ=255)                             gcgcCTCGCCGCGCTGGCTGGCGGTTTCAATCACCACCCAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCa                                              <  2:939233/81‑1 (MQ=255)                              cgcCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCAc                                             >  2:939303/1‑81 (MQ=255)                              cgacTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCAc                                             <  2:939234/78‑1 (MQ=255)                               gcCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCAcc                                            >  1:939308/1‑81 (MQ=255)                                   cgccgcGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCaa                                         >  1:939297/1‑80 (MQ=255)                                    gccgcGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCaaaa                                       <  2:939235/81‑1 (MQ=255)                                    gccgcGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCaaa                                        <  2:939236/80‑1 (MQ=255)                                    gccgcGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCaaa                                        >  1:939315/1‑80 (MQ=255)                                    gccgcGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCaa                                         <  1:939237/79‑1 (MQ=255)                                    gccgcGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCCGTCACCATCGATGCCAGCTTCAATCGCCACCCaaa                                        >  2:939292/1‑80 (MQ=255)                                      cgcgCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAAt                                     >  2:939313/1‑81 (MQ=255)                                       gcgcTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAAt                                     >  2:939311/1‑80 (MQ=255)                                       gcgcTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAAt                                     >  1:939319/1‑80 (MQ=255)                                            ggctggcGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCct                               >  2:939309/1‑81 (MQ=255)                                              ctggcGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCctct                             <  2:939238/81‑1 (MQ=255)                                                 gtgTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTg                            <  2:939239/77‑1 (MQ=255)                                                 gcGTTTTCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTgg                           <  2:939240/80‑1 (MQ=255)                                                      ttCAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTa                      >  2:939298/1‑80 (MQ=255)                                                        cAATCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTaaa                    <  2:939241/80‑1 (MQ=255)                                                         aaTCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAAcg                  >  1:939296/1‑81 (MQ=255)                                                         aaTCACCACCCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAAcg                  <  2:939242/81‑1 (MQ=255)                                                             accaccCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAAcggcgg              <  2:939243/81‑1 (MQ=255)                                                             accaccCAACTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAAcggcgg              <  1:939244/81‑1 (MQ=255)                                                                      cTGAAAGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCt     <  1:939245/81‑1 (MQ=255)                                                                         aaaGTTCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTAcc  <  1:939246/81‑1 (MQ=255)                                                                              |                                                                            ACCGCAGGCAACGGCAAAGTAGCCGAACGCGCCTCGCCGCGCTGGCTGGCGTTTTCAATCACCACCCAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACC  >  NC_000913/4633420‑4633572

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