Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,216	A→G	noncoding (71/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,216	0	A	G	86.5%	103.7 / 10.6	37	noncoding (71/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/21);  minor base . (2/3);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

CGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTA  >  NC_000913/4606084‑4606351                                                                                                                                     |                                                                                                                                         cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTTATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAAc                                                                                                                                   >  1:1204492/1‑140 (MQ=255)      ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCa                                                                                                                              <  1:12419/140‑1 (MQ=255)        ggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACg                                                                                                                            <  1:107551/140‑1 (MQ=255)            aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGt                                                                                                                       <  1:5065/140‑1 (MQ=255)             cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTTATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGt                                                                                                                       <  2:1204492/140‑1 (MQ=255)             cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGt                                                                                                                       >  2:965342/1‑140 (MQ=255)              ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTa                                                                                                                      <  2:442894/140‑1 (MQ=255)                      ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACt                                                                                                              <  2:355298/140‑1 (MQ=255)                              gTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                                                                      <  1:965342/140‑1 (MQ=255)                               tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                                                                    >  1:751102/1‑140 (MQ=255)                                       cTAACACCTGAAGCTGGCGCGTCTTCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                             <  2:728146/140‑1 (MQ=255)                                        tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                                                                            >  2:696645/1‑140 (MQ=255)                                             ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCCTTAACCCCCACGTCCGTAAGAACCCTAACACCTGAAGCTAGCGCGTCTa                                                                                       >  2:965644/1‑140 (MQ=255)                                                gAAGCTAGTGCGTCTACCAATTCCGGCACCTTCGCTCAATCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                                    <  2:979364/140‑1 (MQ=25)                                                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                             >  1:20314/1‑129 (MQ=255)                                                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                             <  2:20314/129‑1 (MQ=255)                                                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                                                  >  2:343797/1‑140 (MQ=255)                                                       gcgtGTGTAGCAATTCCTCCGCCTTCGCGCCGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                             <  2:980017/140‑1 (MQ=18)                                                       gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                             <  1:930233/140‑1 (MQ=255)                                                                 cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                                   <  2:955183/140‑1 (MQ=255)                                                                 cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                                   <  1:343797/140‑1 (MQ=255)                                                                     tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                               >  2:127829/1‑140 (MQ=255)                                                                        ggcACCTTAGCTCAGTCATCTAACTTTTTTTGATATCGCCACGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                            <  2:1052458/138‑1 (MQ=18)                                                                        gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                            <  2:941360/140‑1 (MQ=255)                                                                        gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGa                                                           <  2:751102/140‑1 (MQ=255)                                                                          cACCTTCGCTCAGTCATCTTACTTTTTTTGATATCGCCACGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                          <  2:574662/140‑1 (MQ=255)                                                                              ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                      >  1:41470/1‑140 (MQ=255)                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                             >  2:1064178/1‑140 (MQ=255)                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                             >  1:770695/1‑140 (MQ=255)                                                                                         aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                                >  2:17144/1‑134 (MQ=255)                                                                                         aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                                <  1:17144/134‑1 (MQ=255)                                                                                            ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                        >  1:697833/1‑140 (MQ=255)                                                                                                   ttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                 <  1:965644/140‑1 (MQ=255)                                                                                                        atatCGCCTCGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑gggg                           <  1:290045/140‑1 (MQ=255)                                                                                                                 cGTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTg                  <  2:770695/140‑1 (MQ=255)                                                                                                                  gTTTGGTGCGAGGGGGGGGGCTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGa                 <  2:1091230/140‑1 (MQ=255)                                                                                                                      ggTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAAccc             >  2:722460/1‑140 (MQ=255)                                                                                                                        tGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAAccccc           <  1:696645/140‑1 (MQ=255)                                                                                                                         gCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCa          <  1:127829/140‑1 (MQ=255)                                                                                                                                  ggggCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTa  <  1:1064178/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                         CGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTA  >  NC_000913/4606084‑4606351

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