Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,049	T→A	16.0%	intergenic (‑189/+10)	argQ ← / ← argZ	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,049	0	T	A	16.0%	153.8 / 15.6	55	intergenic (‑189/+10)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (30/16);  new base A (3/6);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGA  >  NC_000913/2817902‑2818197                                                                                                                                                    |                                                                                                                                                     taAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg                                                                                                                                                    >  1:515203/1‑150 (MQ=255)   aaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCAAGCGCCTCGCCCTACGGGCCGATGCTAAAGCAACGATATCATCCCTAGAGCCTGCGGTTCATTCCGCTACATCTACGATACCTTCATTGCTGAAGATTTCGga                                                                                                                                                  >  2:782799/1‑150 (MQ=18)        atatCAATTAAGAGATGGTGCATCCGGGTGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAata                                                                                                                                             <  1:984320/150‑1 (MQ=255)           ccAAATAACAGATGCTGCACCCGCGAGGAGTACTCGCCTTGCCGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAACTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATgg                                                                                                                                          <  1:1130772/149‑1 (MQ=18)              aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg                                                                                                                                                    >  1:741941/1‑137 (MQ=255)              aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg                                                                                                                                                    <  2:741941/137‑1 (MQ=255)                tAATAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCAt                                                                                                                                     <  2:540717/150‑1 (MQ=255)                    agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                                 <  1:47379/150‑1 (MQ=255)                         gTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggagga                                                                                                                            >  2:1107559/1‑150 (MQ=255)                          tGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCATGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGAt                                                                                                                           >  1:1298348/1‑150 (MQ=255)                           gcctCCGGGCGGATGCCGCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCAGCCCTCGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGAtt                                                                                                                          <  2:398891/147‑1 (MQ=17)                             aTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCg                                                                                                                        >  2:1255475/1‑150 (MQ=255)                                 gggAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAAcc                                                                                                                    >  1:1257275/1‑150 (MQ=255)                                 gggAGGATGACTCAC‑‑TGCCGCTCGCCCGTCGGGCCGTTGCTCCCGCACCGTTATCATCCCTAGTGCCTGCGGCTAAGTCCGCTACCTCTTCGATACCTCCATTGCTGAAGTTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTc                                                                                                                  <  2:114776/150‑1 (MQ=11)                                   gaggaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAAGTTATAATCCCTGGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGGAGATTTCGGAGGTTATGGTGCTTCCGGGAGGGTTCGAACCTc                                                                                                                  >  2:1387121/1‑150 (MQ=17)                                      gaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                               >  2:71542/1‑150 (MQ=255)                                                  tgacgttCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCg                                                                                                   >  2:534398/7‑150 (MQ=255)                                                     ggCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATAATCCCTAGTGCCTGCGGTTAATTCCGCTACCACTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                                                >  1:281673/1‑150 (MQ=255)                                                       cTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                                              <  1:1342295/150‑1 (MQ=255)                                                       cTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                                              <  1:594263/150‑1 (MQ=255)                                                                   gCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                  >  2:775288/1‑150 (MQ=255)                                                                   gCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                  >  2:907047/1‑150 (MQ=255)                                                                     cGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                                                >  1:1183525/1‑150 (MQ=255)                                                                        tGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCGACATCTTAGATACCTGCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATACCAACCTCCGACCGAGCGGGTCGTGGCCGAGTCGTCTATCCAGCt                                                                             >  2:1349817/1‑150 (MQ=17)                                                                                  aCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                                   >  1:650251/1‑150 (MQ=255)                                                                                          atcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                                                           <  1:534398/150‑1 (MQ=255)                                                                                             ccTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaa                                                        >  2:994670/1‑150 (MQ=255)                                                                                              cTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaa                                                       >  1:906823/1‑150 (MQ=255)                                                                                                    ccTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttatt                                                 >  1:356323/1‑150 (MQ=255)                                                                                                    ccTGCGGTTAATTCCGCTACCTCGTCGAGACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttatt                                                 <  1:1140446/150‑1 (MQ=255)                                                                                                     cTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTCGCCGAGTACACTATCCAGCTGAGCTACGGATGCATCGGGAAatttattt                                                >  1:361296/1‑150 (MQ=255)                                                                                                      tGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTAgttt                                               >  2:35577/1‑150 (MQ=255)                                                                                                       gCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGGACTCTATCCCGCTGAGCTCCGGAGGCAGCGGGAAAGTTATTTTa                                              >  1:209764/1‑150 (MQ=14)                                                                                                           ttAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCAGACCGCTCGCGTCGTGGCCGGGTACTCTGTACCGCTGAGCTACGGATGCATCGGGAAAGGTAGGTTActga                                          >  2:1309437/1‑149 (MQ=14)                                                                                                              aTTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGATCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGAATCGGGAAATTTATTTTACTGCTGa                                       >  1:774348/1‑150 (MQ=255)                                                                                                                tCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGata                                     >  2:409241/1‑150 (MQ=255)                                                                                                                       tctctTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAc                              <  1:260618/149‑1 (MQ=17)                                                                                                                         tctTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                            >  2:109058/1‑150 (MQ=18)                                                                                                                         tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                            >  2:1263833/1‑150 (MQ=255)                                                                                                                           ttCGATACCTTCATTGCTGAAGGTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGct                          <  1:503896/150‑1 (MQ=38)                                                                                                                           ttCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGct                          <  1:503870/150‑1 (MQ=18)                                                                                                                            tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                         <  2:281673/150‑1 (MQ=255)                                                                                                                                 aCCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTATTTTACTGCTGATATTGATACCGCTCCTaaa                    >  1:1199347/1‑150 (MQ=255)                                                                                                                                  ccTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg                   <  2:628561/150‑1 (MQ=18)                                                                                                                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                    <  1:79177/148‑1 (MQ=18)                                                                                                                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                    >  2:79177/1‑148 (MQ=18)                                                                                                                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc                  >  2:543344/1‑150 (MQ=18)                                                                                                                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc                  <  1:950466/150‑1 (MQ=18)                                                                                                                                        ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                                            >  2:445882/1‑71 (MQ=39)                                                                                                                                        ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                                            <  1:445882/71‑1 (MQ=39)                                                                                                                                           cTGAAGGTTACGGAGAATATGGTGCATCCGGGTGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa          <  2:5008/150‑1 (MQ=16)                                                                                                                                           cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa          <  2:732862/150‑1 (MQ=255)                                                                                                                                              aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                                                 >  2:1402230/1‑60 (MQ=255)                                                                                                                                              aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                                                 <  1:1402230/60‑1 (MQ=255)                                                                                                                                                gATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATGCAGCTACTATAGACATATCAAGTAag     >  2:1047376/1‑150 (MQ=255)                                                                                                                                                   ttCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAgaga  <  2:650251/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGA  >  NC_000913/2817902‑2818197

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