Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,046	A→G	18.7%	intergenic (‑186/+13)	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,046	0	A	G	18.7%	128.6 / 21.4	55	intergenic (‑186/+13)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (30/14);  new base G (3/7);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 3.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAG  >  NC_000913/2817902‑2818194                                                                                                                                                 |                                                                                                                                                     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                 >  2:79177/1‑148 (MQ=18)                                                                                                                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                 <  1:79177/148‑1 (MQ=18)                                                                                                                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc               <  1:950466/150‑1 (MQ=18)                                                                                                                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc               >  2:543344/1‑150 (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)                                                                                                                                                 |                                                                                                                                                     TAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAG  >  NC_000913/2817902‑2818194

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: ‑