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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,094	1	.	G	95.9%	157.7 / ‑0.8	49	noncoding (72/87 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (19/28);  total (20/29)
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.58e-01

GCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGG‑‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAC  >  NC_000913/4605854‑4606344                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                             gCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGAGGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAAc                                                                                                                                                                                                                                                    <  1:365242/250‑1 (MQ=255)  cTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAAccc                                                                                                                                                                                                                                                  >  2:292915/1‑251 (MQ=255)   tGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAAcccc                                                                                                                                                                                                                                                 >  2:87242/1‑251 (MQ=255)            ggcggATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGt                                                                                                                                                                                                                                        <  2:576270/247‑1 (MQ=255)                gATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGGCCGTAAgg                                                                                                                                                                                                                                    >  2:173474/1‑251 (MQ=255)                       gCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCCAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGAAAGGACCCTa                                                                                                                                                                                                                              >  2:392240/1‑250 (MQ=255)                       gCGGATGTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCAGGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATGGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                                                                                                                                             <  1:697218/251‑1 (MQ=255)                         ggATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAac                                                                                                                                                                                                                            >  1:557518/1‑250 (MQ=255)                              tATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGGTGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGAAACTAACACCTg                                                                                                                                                                                                                       >  2:229232/1‑250 (MQ=255)                                 cGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                                                                                                                                                                                                   >  1:137672/1‑251 (MQ=255)                                 agTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                                                                                                                                                                                                                      >  2:400418/2‑249 (MQ=255)                                                    tGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGGGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                                                                                                                                                 <  1:590442/250‑1 (MQ=255)                                                       cagtCTGATGGGATTAGTAGTGACGACGATGAGGCGGTAGCACAAACCGTTGCGGCATTACATTTTTCTACCAACGGGCCTGGCACTTATTGCTGGTATCGTGCTGATGATAAGGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGGGAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCGGAAGCGAGCGCGTCGACAGATTCCGc                                                                                                                                                                                              <  2:324545/247‑1 (MQ=255)                                                             gATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                                                                                                                                                          <  1:132091/248‑1 (MQ=255)                                                                gggATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCTGCAGTACATTTGGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                                                                                                                                                    <  1:120826/251‑1 (MQ=255)                                                                  gATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAAGTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAAGTGGTGCGAGGTGGGGGG‑‑A‑CTTGAGCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                                                                                                                                                                                  <  2:722437/251‑1 (MQ=255)                                                                  gATTAGGATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGATCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                                                                                                                                                                                  <  1:292915/251‑1 (MQ=255)                                                                         aTTGACGCCGATGATGCTGTATCACAAACCGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACTCATCt                                                                                                                                                                           >  2:249110/1‑251 (MQ=255)                                                                              cGCCGATGTTGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTGCCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGAGGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtt                                                                                                                                                                      <  2:182388/251‑1 (MQ=255)                                                                               gCCGATGATGCTGTATCACAAACAGTTGCGGCATTACCTTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCACGGAAGAAACAAGAAAGGAAGTATAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAAGACACTAAAACCTGAAGCAAGCGCGTCTACAAATTCCGCCACCTTCGCACAGTCATCTTACttt                                                                                                                                                                     >  1:392929/1‑251 (MQ=38)                                                                                  gatgatGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTCTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGG‑‑GGCTTGAACCCCCACGTCCGTAAGGACACTAACCCCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTTCGCACAGACATCTTACtttttt                                                                                                                                                                  >  2:654264/1‑251 (MQ=38)                                                                                     gatgCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTg                                                                                                                                                                <  2:215176/250‑1 (MQ=255)                                                                                     gatgCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAAGCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTg                                                                                                                                                                <  2:179300/250‑1 (MQ=255)                                                                                      atgCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGat                                                                                                                                                              >  1:47886/1‑251 (MQ=255)                                                                                               cacaAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                     <  2:557518/251‑1 (MQ=255)                                                                                                     cAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                                                                                                                                                                                               >  1:114460/1‑187 (MQ=255)                                                                                                     cAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                                                                                                                                                                                               <  2:114460/187‑1 (MQ=255)                                                                                                       gTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCCCTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGc                                                                                                                                             <  1:10142/251‑1 (MQ=255)                                                                                                          gCGGCATTTCATTTTGCTACCATCGTGCCTGGCAATTTGTGCTGGCATAATGCTGATGGTAATGAATTTGAATCAGGGATGAAAAATCAAGTAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGGGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATGCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                                                                                                          <  2:730029/251‑1 (MQ=37)                                                                                                             gCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                       >  1:302551/1‑251 (MQ=255)                                                                                                               aTTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                                     >  1:470041/1‑251 (MQ=255)                                                                                                                         gCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAAc                                                                                                                           <  1:300162/251‑1 (MQ=255)                                                                                                                           tACCATCGTGCATGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCAAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAAcc                                                                                                                          <  2:363952/250‑1 (MQ=255)                                                                                                                            cccATCTTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAAACAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAAcc                                                                                                                          <  2:214071/248‑1 (MQ=255)                                                                                                                               agcGTGCCTGGCCCTTATTGCTGGCATCATGCTTATTATAATGAATTTTAATCAGGGATGAAAAATCAAGGAAGAAACAAAAAAGGAAGTAAAGATAATTGGGGGGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc                                                                                                                     <  2:302551/249‑1 (MQ=37)                                                                                                                                 cGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACg                                                                                                                    <  2:47886/250‑1 (MQ=255)                                                                                                                                    gCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                                                                                                                <  1:147520/251‑1 (MQ=255)                                                                                                                                    gCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                                                                                                                >  1:389878/1‑251 (MQ=255)                                                                                                                                    gCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTGGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                                                                                                                <  1:111169/251‑1 (MQ=255)                                                                                                                                                 tGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                    <  1:311742/250‑1 (MQ=255)                                                                                                                                                     ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                               <  1:373596/251‑1 (MQ=255)                                                                                                                                                                     tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGTCTCGTTTGGTGCGAGGGTGGGGACTTGAACCCCCACGTCCGTAATAACACTAACACCTGGAGCTAGCCCGTCt                                                                                >  1:449247/1‑250 (MQ=38)                                                                                                                                                                           aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                         <  1:411370/251‑1 (MQ=255)                                                                                                                                                                               gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                      <  1:400418/251‑1 (MQ=255)                                                                                                                                                                                        aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                            <  1:464902/251‑1 (MQ=255)                                                                                                                                                                                                         aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATCTTACTTATTTTGATATCGCTTCGTTGGGCGCGAGGGGGGGGACTTGAACCCCCCCGTCCGTACGAACACTAACACCTGAAGCTAACGCGTCTACCAATTCCGCCACCTTCCCCCAGTGCGAACAATATc                                           >  2:227059/1‑251 (MQ=25)                                                                                                                                                                                                                 aaGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGAAAGGAAACTAACACATGAGGCCAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTAATTTTTTCGATATCGCCTAGTTTGGTGCGAGGGGGGGGACTTGAACACCCACGTCCGTAAGGAAACTAACACCCGAAGCTAGCGCGTCGACCAATTCCGCCAGCTTCGCCAAGTGAGCGCAATATCAACGCGGt                                   >  1:56593/1‑251 (MQ=21)                                                                                                                                                                                                                                aTTGGTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACTCTAACACCTGACGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGCGTTTGGAGCGAgggg                    >  1:183592/1‑251 (MQ=37)                                                                                                                                                                                                                                  ggggcgaggg‑gggggg‑‑A‑CTTGACCCCCCGCGTCCGTAAGGACACTAACACCTGATGCTAGCGCGTCTACCGATGCCGCCACCTTCGGACATTCATCTTACTGTTTGTGATGTCGCCTCGTTTGGTGCGAGGGGGGGGACTTGATCCCCCACGTCCGTAAGAACACTAACGCCTGAAGCTAGCGCGTCTACCAATTCCGGCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGGGCGAggggggg                 <  1:364020/243‑1 (MQ=255)                                                                                                                                                                                                                                  gggTGCGAGG‑GGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg                  <  1:666899/250‑1 (MQ=255)                                                                                                                                                                                                                                     gcgaggg‑gggggg‑‑A‑CTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAACAACGTGGTTTTTGGTGCGAGGGGGGGGa               <  1:229232/245‑1 (MQ=255)                                                                                                                                                                                                                                     gcgaggg‑gggggg‑‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTACCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATGTCAACGTGGTTTTTGGCGCGAgggggtgggc              >  1:27055/6‑246 (MQ=255)                                                                                                                                                                                                                                       gaggg‑gggggg‑‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCCTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGAGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGTGCAATATCAACGTGGTTTTTGGACCGAGGGGGGGGCCtt            >  1:139171/4‑251 (MQ=37)                                                                                                                                                                                                                                          gg‑gggggg‑‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTTCGAGCAATATAAAGGTGGCTTTTGTTGCGAGGAGGGGGACTTGa          >  1:427162/1‑250 (MQ=37)                                                                                                                                                                                                                                                  g‑‑A‑CTTGACCCCCCACGTCCGTAAGGACAGTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc  <  1:268812/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                             GCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGG‑GGGGGG‑‑A‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAC  >  NC_000913/4605854‑4606344

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