Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,215	(G)8→9	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,207	1	.	G	96.0%	162.8 / 0.1	50	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (27/21);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606068‑4606341                                                                                                                                              |                                                                                                                                            gTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCCTCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑Ga                                                                                                                                              >  1:1180086/1‑140 (MQ=255)  tAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑g                                                                                                                                           <  2:1158424/140‑1 (MQ=255)      gATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑ggg‑‑g‑g                                                                                                                                    <  1:226799/140‑1 (MQ=255)       aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑ggg‑‑g‑gg                                                                                                                                   <  2:1486199/140‑1 (MQ=255)             ggTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑ggg‑‑g‑ggggggac                                                                                                                             <  1:571982/140‑5 (MQ=255)              gTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑ggg‑‑g‑ggggggctt                                                                                                                            >  1:954251/1‑135 (MQ=255)               ggCGGGGGGGGGGACTTCAACCCCCCCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTg                                                                                                                           <  2:896279/139‑1 (MQ=21)                    ggggggggGCTTGAACCCCCCCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGGCTTCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAAccc                                                                                                                      <  2:5962/140‑1 (MQ=25)                       gggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCa                                                                                                                   <  1:592591/140‑1 (MQ=255)                        ggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCAc                                                                                                                  <  1:1246241/140‑1 (MQ=255)                         gggACTTGAACCCCCGCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACg                                                                                                                 <  2:775393/140‑1 (MQ=255)                         gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCGGAAGCTAGCGCGACCACCAACTCCGCCACCTTCGCCCAGTCAACTTACTGTTGGTGATGGCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTTAACACCCACg                                                                                                                 >  2:141601/1‑140 (MQ=17)                              ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTCATCTTACTTTTTTTGATATCGCCTCGTCTGGTGC‑GA‑‑GGGT‑G‑GGGG‑ACTTGAACCCCCACGTCCga                                                                                                            >  2:343546/1‑139 (MQ=255)                                  aaccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGa                                                                                                        >  1:1510390/3‑140 (MQ=255)                                      ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGG‑‑G‑GGGGCTCTTGAACCCCCCCGTCCGCAAGTacac                                                                                                    >  2:778151/1‑140 (MQ=25)                                         cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                 <  1:1328958/140‑1 (MQ=255)                                                    gcaggaaCACCTGAAGCTAGCGAGTCTACCAATTCCGCCACCTTCGCACAGGCATCTTACTATTTTGGATGTCGCCTCGTTTGGTGC‑GA‑‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                     <  2:517297/135‑1 (MQ=25)                                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATCTTACTTTTTAAGAAATCCCCCCGTTGGGGGC‑GA‑‑GGG‑‑GTGGGC‑ACTTCAACCCCCACCGACGAACGACCCCTTACACCTGAAGCTAg                                                                                    >  2:1486729/1‑140 (MQ=16)                                                          acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGTTATCGCCTCGTTTGGTGC‑GA‑‑GGGCTG‑GGGG‑ACTTGCATCCCCACGTCCGCTAGTATACTAATCCCTCAAGCtggatc                                                                                 >  2:37581/1‑135 (MQ=21)                                                                    cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                      >  1:327914/1‑140 (MQ=255)                                                                      aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                    >  2:1314500/1‑140 (MQ=255)                                                                        cgcgTGTACCAAATCCGGCGCCGACGAACAGTCATCTTACTTTTTTTGATATCGCCGCGATTGGGGCGGG‑‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:971358/140‑1 (MQ=18)                                                                        cgcgTCTACCAATTCCGCCGCCTACGCACAGTCAACTTACCTTTTTTGATAGCGCCGCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:634058/140‑1 (MQ=21)                                                                        cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  >  1:281868/1‑140 (MQ=255)                                                                            tcaaCCAATTCCGCCGCCTTCGCACAGTCATCTAACTATTATTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                              <  2:1473360/137‑1 (MQ=21)                                                                            tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                              <  1:404413/140‑1 (MQ=255)                                                                               aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAAACCCAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                           >  2:267459/1‑140 (MQ=25)                                                                                 cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                         <  1:1314500/140‑1 (MQ=255)                                                                                  aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  1:950897/140‑1 (MQ=255)                                                                                  aaTTCCGCCACCTTCGAACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  1:1422683/140‑1 (MQ=255)                                                                                  aaTTCCGCCACCTTCAGTCAATGAACTTACTTTGGTTGAAATCGCCGTGATTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  2:985553/140‑1 (MQ=14)                                                                                   aTTCGGCCACCTTCGCACGGTCATGTTGCTTTTTTGGATATCGCCTCGTTTGGGGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                       <  2:954251/140‑1 (MQ=18)                                                                                    ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCTCCTCGGTTGGTGA‑GA‑‑GGG‑‑G‑GGGGCACTTGAACCCCCAACGCCGCAGGCACACCAACACCTGAAGCTAGCGCGTCTACCACTTCCGCCACCTTCGCCCa                                                      >  2:324007/1‑140 (MQ=21)                                                                                     tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGGG‑‑G‑GGGG‑ACTTGAACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                                                      >  2:144407/1‑140 (MQ=255)                                                                                     tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑CCTTGAACCCCCAAGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                                     >  1:255159/1‑140 (MQ=255)                                                                                       cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   >  1:153263/1‑140 (MQ=255)                                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGG‑‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                 <  1:184266/140‑1 (MQ=255)                                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                 >  2:123591/1‑140 (MQ=255)                                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                 >  1:622643/1‑140 (MQ=255)                                                                                          cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGa                                                >  2:640127/1‑140 (MQ=255)                                                                                               tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                           >  1:1033953/1‑140 (MQ=255)                                                                                                 gCACAGTCAACTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat                                         <  2:881050/140‑1 (MQ=255)                                                                                                     aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                     <  1:694145/140‑1 (MQ=255)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAACAACg                                    >  2:1247396/1‑140 (MQ=255)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACgg                                   >  1:362490/1‑139 (MQ=255)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACga                                   >  1:835105/1‑139 (MQ=255)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                   >  1:924486/1‑140 (MQ=255)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                   >  1:611294/1‑140 (MQ=255)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGG‑‑G‑GGGGAACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTATCGCGTCAACCAATTCCGCCACCCTCGCCCAGTACCAACAATATGAACga                                   >  2:1008075/1‑139 (MQ=18)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                   <  1:778151/140‑1 (MQ=255)                                                                                                          tCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                >  2:501476/1‑140 (MQ=255)                                                                                                           cTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCCAGCTATAACAACGAGGtt                               >  2:147897/1‑140 (MQ=18)                                                                                                            ttACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttt                              <  1:343546/140‑1 (MQ=255)                                                                                                                 ttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                         >  1:430284/1‑140 (MQ=255)                                                                                                                 ttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                         >  1:1052640/1‑140 (MQ=255)                                                                                                                 ttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                         >  1:707714/1‑140 (MQ=255)                                                                                                                 ttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                         >  1:462725/1‑140 (MQ=255)                                                                                                                 ttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTACGAGCAATATCAACGTGGTTTTTGGt                         >  2:564981/1‑140 (MQ=255)                                                                                                                    tttGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCCCTGCCAACAATATCCACGTGGTTTTTGGTGcc                      >  2:78333/1‑139 (MQ=18)                                                                                                                     ttGATATCGCCTCGTTTGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                     <  1:267459/140‑1 (MQ=255)                                                                                                                     ttGATATCGCCTCGTTTGGGGC‑GAGGGGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                      <  1:207182/140‑1 (MQ=255)                                                                                                                            cGCCTCGTTTGGTGC‑GAGGGGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg               <  1:144407/140‑1 (MQ=255)                                                                                                                            cGCCTCGTTTGGTGC‑GAGGGGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCGCGAGCACTATCAACGAGGTTTTTCGTGCGAgggggg               >  2:683798/1‑140 (MQ=18)                                                                                                                                   tttGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACTCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGa       <  1:1456858/140‑1 (MQ=255)                                                                                                                                    ttGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGACCCCCCACGTCCGTAAGAACACTATCACCTGAGGCTAGCGCGTCTTCCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGaa      <  1:324007/140‑1 (MQ=21)                                                                                                                                     tGGTGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAc     <  2:8711/140‑1 (MQ=255)                                                                                                                                        tGC‑GAG‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATCAACGAGGTTTTTGGTGCGAGGGGGGGGCCTGGGAcccc  >  2:1105017/1‑140 (MQ=18)                                                                                                                                              |                                                                                                                                            GTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGG‑‑G‑GGGG‑ACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606068‑4606341

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