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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,322	1	.	G	84.1%	165.1 / 24.5	63	noncoding (80/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (5/5);  new base G (22/31);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑A‑GG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCT  >  NC_000913/4606174‑4606461                                                                                                                                                        |                                                                                                                                                cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑A‑gg                                                                                                                                               >  2:505711/1‑151 (MQ=255)    aCTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑A‑gg‑gg                                                                                                                                            <  1:738575/150‑1 (MQ=255)     cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTT‑GGTGCG‑A‑gg‑ggg‑g                                                                                                                                         >  1:285281/1‑151 (MQ=255)            tGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑A‑gg‑ggg‑gg‑ggcctt                                                                                                                                 >  1:642041/1‑146 (MQ=255)              atatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGa                                                                                                                               <  2:747949/151‑1 (MQ=255)               tatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGa                                                                                                                               >  2:678310/1‑150 (MQ=255)                      tCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCGCGAGCAATATCAACGTGGTTTTT‑GGAGCG‑A‑GGCGGG‑GG‑G‑ACTTGAACCCCCAc                                                                                                                       >  1:777110/1‑151 (MQ=255)                       cGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCa                                                                                                                        <  2:210726/149‑1 (MQ=255)                       cGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTAGTTTTT‑GGTGCG‑A‑GG‑GGGCGG‑G‑ACTTGAACCCCCAc                                                                                                                       >  1:448436/1‑150 (MQ=255)                         tttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCAc                                                                                                                       <  1:971725/148‑1 (MQ=255)                          tgggTGCGTGGGGGGGGACTTGAACCGCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATAAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCa                                                                                                                        <  1:938857/144‑1 (MQ=21)                           tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCACGAGCAATATCAACGAGGTTCTT‑GCACCG‑A‑GGCGGG‑GG‑C‑ACTTGAACCCCCACGt                                                                                                                     >  1:81401/1‑148 (MQ=17)                            ggTGCGAGGGACGCGGCTTGCCCCCCCCCGTGCGCAACAACACTACCACCTGAAGCTAGCCCGTCTACCAATTCCGCCACCTACGCCCGGCGCGAGCACTATCACCGTGGTTTTT‑GGTTCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTc                                                                                                                    <  1:156590/148‑1 (MQ=16)                              tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTa                                                                                                                >  2:854308/1‑150 (MQ=255)                               gCGGGGGGGGGGCATTCAACCCCACAGTCCGTAAGAACACTAAGACCTGAAGCTAGCGCTTCTTCCAGTTCCGCCACCTGCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAg                                                                                                              <  1:644924/151‑1 (MQ=25)                               gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTaa                                                                                                               <  1:968701/150‑1 (MQ=255)                                  aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGac                                                                                                           >  2:142366/1‑151 (MQ=255)                                  aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGac                                                                                                           >  2:118824/1‑151 (MQ=255)                                   gggggtggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTCCGCCCAGTACGAGCAATAGCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCa                                                                                                                        <  1:685690/137‑1 (MQ=21)                                     ggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTCCGCCCAGTGCGAGCTATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                        <  1:497128/151‑1 (MQ=255)                                      gggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTa                                                                                                       <  2:395117/151‑1 (MQ=255)                                       ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                        >  2:1147384/1‑149 (MQ=255)                                        gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAac                                                                                                     <  2:730694/151‑1 (MQ=255)                                           aCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCGCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                                                                  <  1:560488/151‑1 (MQ=255)                                              tGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                                                                 <  2:910110/149‑1 (MQ=255)                                               gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTGGCGCGTCTACCAATTCCGCCACCTTCGCCCGGTGCGTGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGaa                                                                                              <  2:607286/151‑1 (MQ=255)                                                  cccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCACTATCAACGTGGTTTTT‑GCTCCG‑A‑GGCGCG‑GG‑C‑TCTTGAACCCCCACTTCCGTACGGTCACTAactctcga                                                                                               >  1:13417/1‑141 (MQ=18)                                                   ccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTGCGGCCAGTGCGAGCAACATCAACGTGGTTTTT‑GCGGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                          <  1:402281/151‑1 (MQ=37)                                                       aCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATCAACGTGGTTTTT‑GGGGCG‑A‑GGCGGG‑GG‑G‑ACGTGAACCCCCACGTCCGTAAGGATACTAACACCTGAAGCTAgcg                                                                                       >  1:156804/1‑150 (MQ=37)                                                        cGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTACGGACACTAACACCTGAAGCTAgcg                                                                                       >  1:155864/1‑149 (MQ=255)                                                        cGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                                     <  2:423387/151‑1 (MQ=255)                                                             gTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTGTGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACTCCCACGTCCGGAAGGTCACTAATACCTGAAGGTAGCGCGTCTAc                                                                                >  1:714615/1‑151 (MQ=25)                                                                      cTAACACCTGAAGCTAGCGCGTCTCCCAATTCCGCCACCTCCGCCCAGTGCGAGCAATAGCAACGTGGTTATT‑GGTGGG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                             <  1:1189507/145‑1 (MQ=25)                                                                        aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                                     <  2:139273/151‑1 (MQ=255)                                                                         acacCAGAAGATAGCGCTTCTACCAATTACGCAAACTCAGCCAAGTGCGAGCAATAGCAACGTGGATTTT‑TGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                         <  1:681478/146‑1 (MQ=18)                                                                           acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                                                                  >  2:767011/1‑151 (MQ=255)                                                                                         cgtcctccACTTCCGCCACCTTCGCCCGGTGCGAGCAATATCAACGTGGTTTTT‑TGTGCC‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                                    <  1:887015/145‑1 (MQ=18)                                                                                            cTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTc                                                 >  2:942765/1‑151 (MQ=255)                                                                                               ccAAGTCCGCCACCTGCGCCCACTGCCTGCAATAGCAACGTGGTTATT‑GGTGCC‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCtt                                               <  1:1144945/150‑1 (MQ=18)                                                                                                cAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGTTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAtt                                                  <  2:1095245/146‑1 (MQ=255)                                                                                                         cACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                    <  2:156838/151‑1 (MQ=255)                                                                                                          aCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGGGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc                                   >  1:1650/1‑151 (MQ=255)                                                                                                            cTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                    <  2:448543/148‑1 (MQ=255)                                                                                                             ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAc                                >  2:967886/1‑151 (MQ=255)                                                                                                             ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAc                                >  2:896655/1‑151 (MQ=255)                                                                                                              tCGCCCCGTGCGATCATTATCAACGTGGGTTTT‑GGTTCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTa                                 <  1:330904/149‑1 (MQ=21)                                                                                                               cGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGCCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAcc                               >  2:361693/1‑150 (MQ=255)                                                                                                               cGCCAAGTTCGAGCAATAGCAACCTCGTTTTT‑CGTACG‑A‑GG‑GGG‑GGTG‑CCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCa                              <  1:645968/151‑1 (MQ=18)                                                                                                                 cccAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg                            >  1:559448/1‑151 (MQ=255)                                                                                                                  ccAGTGCGGACACTCTCACCCTAGTTGTT‑GGTGCGCG‑GG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt                                  <  1:920502/144‑1 (MQ=17)                                                                                                                   cAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg                           >  2:1169144/1‑150 (MQ=255)                                                                                                                    aGTGCGAGCAATAGCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg                           <  2:1053831/149‑1 (MQ=255)                                                                                                                     gttcGAACAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCAc                             <  1:808520/143‑1 (MQ=255)                                                                                                                        cgggCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgc                       <  2:972875/146‑1 (MQ=255)                                                                                                                        ccaGCAATAGAAACATGGTTTTT‑GGGGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg                      <  2:243807/148‑1 (MQ=21)                                                                                                                         gAGCAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcgc                     >  2:533974/1‑150 (MQ=255)                                                                                                                            cAATATCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTc                 >  1:1152422/1‑151 (MQ=255)                                                                                                                              atatCAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCt                >  2:105843/1‑150 (MQ=255)                                                                                                                                atCAACCGGGTTTTT‑GGTGCG‑A‑GG‑GGG‑GG‑GCACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAtt                  <  1:1137544/146‑1 (MQ=255)                                                                                                                                 tCAACGTCGTTTTTCGGGTCA‑A‑GG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTc                 <  1:704612/146‑1 (MQ=21)                                                                                                                                  cAACGTGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTgg           >  1:376532/1‑151 (MQ=255)                                                                                                                                    cctggttctttgtgTGCG‑A‑GG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtt         <  1:327559/138‑1 (MQ=255)                                                                                                                                       tGGTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtttt       <  1:1147124/150‑1 (MQ=255)                                                                                                                                         gTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTc      >  2:836693/1‑149 (MQ=255)                                                                                                                                         gTTTTT‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGACGCGCATTCTAGTGGTTTTc      >  2:90891/1‑149 (MQ=255)                                                                                                                                          ttgtt‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCa     <  1:518904/146‑1 (MQ=255)                                                                                                                                           tttt‑GGTGCG‑AGGG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCt  >  2:61150/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                CTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCG‑A‑GG‑GGG‑GG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCT  >  NC_000913/4606174‑4606461

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