Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,215	(G)8→7	noncoding (72/87 nt)	leuP ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,208	0	G	.	97.5%	167.8 / ‑4.3	40	noncoding (79/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (16/23);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAG  >  NC_000913/4606062‑4606353                                                                                                                                                   |                                                                                                                                                  aaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgtggg                                                                                                                                               >  2:258933/1‑151 (MQ=255)        tAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTg                                                                                                                                       >  2:409760/1‑151 (MQ=255)                aTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAAccccc                                                                                                                                >  2:131001/1‑150 (MQ=255)                        gAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCAATCTTACTTTTTTAGATATCGCCTCGTTTGGTGCG‑GCGGGGGGACTTTACAGCCCACGt                                                                                                                            >  2:236355/1‑146 (MQ=18)                                        aCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                                                                       >  2:432400/1‑151 (MQ=255)                                          ccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                                                                     <  2:656249/151‑1 (MQ=255)                                                       aGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                                                                        <  1:628917/151‑1 (MQ=255)                                                               aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                                  <  2:298353/149‑1 (MQ=255)                                                                      gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCATACTTTTTTTGATATCGTCTCGTTTCGCGCCA‑TCGGGGGCCTTGTCACCCCACCTCCGTAAGAACACTAACACCTGACGCGAGCGCGTCGTCCAATTCCCCCAc                                                                           >  1:98214/1‑149 (MQ=21)                                                                           tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCTCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                                    >  1:287882/1‑151 (MQ=255)                                                                                cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCTAg                                                                 >  1:109619/1‑149 (MQ=255)                                                                                 gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                              <  2:352682/151‑1 (MQ=255)                                                                                     aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                             <  2:554180/148‑1 (MQ=255)                                                                                      ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCa                                                         <  2:209728/151‑1 (MQ=255)                                                                                       cAATTCCGCCTCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                           <  1:633229/148‑1 (MQ=255)                                                                                        aaTTCGGCCTCCTTCGCACAGTCATCTTACTTTAATTGATATCGCATCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                       <  1:591406/151‑1 (MQ=21)                                                                                        aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                       >  2:561584/1‑151 (MQ=255)                                                                                                   cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTGCGAGCAATATCAACGAGGt                                            >  1:521564/1‑151 (MQ=255)                                                                                                         acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                      >  2:166779/1‑151 (MQ=255)                                                                                                         acaGGCATCTTACATTTTTTGATATCACATAGTTTGTTGGGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                       <  1:391672/150‑1 (MQ=21)                                                                                                          caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                                     <  2:125603/151‑1 (MQ=255)                                                                                                             tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                  >  2:171421/1‑151 (MQ=255)                                                                                                                 cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCTAGCGCGAGCAATATCAACGTGGTTTTTGGAGCGAgg                               >  1:72624/1‑150 (MQ=255)                                                                                                                   tACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggg                            >  2:53669/1‑151 (MQ=255)                                                                                                                   tACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg                             <  1:591841/150‑1 (MQ=255)                                                                                                                   tACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg                             <  1:287105/150‑1 (MQ=255)                                                                                                                    aCTGTTTTTGATATAGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggg                              <  1:660782/148‑1 (MQ=255)                                                                                                                     cTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg                          >  2:662142/1‑151 (MQ=255)                                                                                                                       ttttttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggg                         <  2:10485/150‑1 (MQ=255)                                                                                                                       ttttatGATAACGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg                          <  1:477929/149‑1 (MQ=255)                                                                                                                        tttttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGa                        <  2:623115/150‑1 (MQ=255)                                                                                                                         ttttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggg                         <  1:80261/148‑1 (MQ=255)                                                                                                                             gATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACt                      <  2:397413/147‑1 (MQ=255)                                                                                                                              atatCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCCTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGaa                  <  1:409310/150‑1 (MQ=255)                                                                                                                              agcaggcgcagTTTGGGACGA‑GGGGGGGACGTGAACCCCCACGTCCGTAAGACCACTAACACCTGAAGCTAGCTCGGCTACCAAGTCCGCCACCATCGCCCAGTGCGAGCAATATCGACGTGGTGTTTGGTGCGAGGGGGGGGACTTGaa                  <  1:52164/140‑1 (MQ=21)                                                                                                                               tatCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCACGAGCAATATAAACGAGGTTCTTCGAGCGAggggggggcca                      >  1:328523/1‑142 (MQ=18)                                                                                                                                atCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGTGAGCAATATCAACGAGGTTTTTGCTCCGAggcgcgtggtcgtg                   >  1:151793/1‑135 (MQ=18)                                                                                                                                    ccTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc           <  2:661390/151‑1 (MQ=255)                                                                                                                                      tCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCa            <  2:255635/148‑1 (MQ=255)                                                                                                                                        gTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTCTTGGAGCGAGGCGGGGGCCTTGACCCCCCACCTcc       >  1:636492/1‑151 (MQ=18)                                                                                                                                             gTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAg  <  2:113434/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  AAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAG  >  NC_000913/4606062‑4606353

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