Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,624	C→T	noncoding (71/76 nt)	lysT →	tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,624	0	C	T	98.0%	142.7 / NA	53	noncoding (71/76 nt)	lysT	tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/22);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAG  >  NC_000913/780382‑780864                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  cgATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATc                                                                                                                                                                                                                                           <  1:515799/251‑1 (MQ=255)  gATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAGTCg                                                                                                                                                                                                                                          >  1:18683/1‑251 (MQ=255)          gCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAACGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTacggcg                                                                                                                                                                                                                                  >  2:36108/1‑246 (MQ=255)           cATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTgg                                                                                                                                                                                                                                 <  1:384391/251‑1 (MQ=255)              aaCACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGaa                                                                                                                                                                                                                               <  1:505086/250‑1 (MQ=255)              aaCACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGGGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAg                                                                                                                                                                                                                              <  2:326932/251‑1 (MQ=255)                   acGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGgt                                                                                                                                                                                                                          >  1:464279/1‑250 (MQ=255)                    cGACCAGAAGTCGCATTATTTCTGGTCGTGGCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGgtag                                                                                                                                                                                                                        <  2:549877/251‑1 (MQ=255)                          gAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTCTCTACGGTGGAAGCGGTAGTAAAAc                                                                                                                                                                                                                  >  1:230186/1‑251 (MQ=255)                               cGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAACCGTGa                                                                                                                                                                                                              >  2:211492/1‑250 (MQ=255)                                      tttCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAAc                                                                                                                                                                                                      >  1:587099/1‑251 (MQ=255)                                           ggTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAGTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTg                                                                                                                                                                                                  <  1:410573/250‑1 (MQ=255)                                                     tgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCTCGACTCACCAATCGCTAAGGTGGAAGAGGTGGTAAAACGTGATGGATAACGTTGCATGAGCAACg                                                                                                                                                                                       >  1:143504/1‑251 (MQ=255)                                                      gCGAATCATAAGCAGTTGAGTGATTTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTTACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACgg                                                                                                                                                                                      <  2:464279/251‑1 (MQ=255)                                                      gCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAACTCAATAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGAACGACTCACCATTCGCTAAGGCGGAAGCGGTAGAAAAACGTGAAGGATAACGTTGCATAAGCAACgg                                                                                                                                                                                      >  2:382978/1‑251 (MQ=38)                                                        gAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGcc                                                                                                                                                                                    <  1:452930/251‑1 (MQ=255)                                                               aaGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAggg                                                                                                                                                                             <  1:185911/251‑1 (MQ=255)                                                               aaGCAGTTGAGTGAGCTACATCGAAATTGTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAGTGGTCGCAGGTGCGAATCCTGCACGACTCACCAATCGCTAAGGGGGAAGCGGGAGTAAAACGTGAAGGATAACGTTGCATGGGCAACGGGCCGAAggg                                                                                                                                                                             <  2:280826/251‑1 (MQ=38)                                                                        aGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGa                                                                                                                                                                     >  1:433099/1‑250 (MQ=255)                                                                              cTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAGCGTGAAGGATAACGTTGCAACAGCAACGGCCCGAAGGGCGAGCCGAACTCggg                                                                                                                                                              >  2:492060/1‑249 (MQ=38)                                                                                aCATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGccc                                                                                                                                                                                   >  1:155251/1‑228 (MQ=255)                                                                                aCATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGccc                                                                                                                                                                                   <  2:155251/228‑1 (MQ=255)                                                                                 cATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCa                                                                                                                                                          <  1:613877/251‑1 (MQ=255)                                                                                     gAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCt                                                                                                                                                      >  2:512995/1‑251 (MQ=255)                                                                                      aaaTTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAGCGTGAAGGATAACGTTGCATGAGCAACGACCCGAAGGGCGAGACGAAGTCGCG‑TCATCCTg                                                                                                                                                     >  1:91998/1‑251 (MQ=255)                                                                                              gttgCGCTCAATTCTGAAATCAGTAATATATTCCTCCGTAGCCACGGTATATCAAACAAAACGAAAGCAACGAAAAAGTGGGGCGTTAGCTCAGTTGGTAGAGCAGTTGACTTACAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGATCAACGGCCCGAAAGGCGAGACGAAGTCGAG‑TCATCCTGCACGAcc                                                                                                                                              <  2:587099/250‑1 (MQ=25)                                                                                                    cTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAGTCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑ACATTCTGCACGACccacctc                                                                                                                                        >  1:542359/1‑248 (MQ=255)                                                                                                        aGTCTGAACTCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTCAGGTGGACGCGGAAGTAAATCGTGAAGGAGAACGTTGCATGCGCAACGGCCCGAAGGGCGAGACGCAGTCGCG‑GCAACCTGAACGACCCCCCAATAAc                                                                                                                                    >  1:328317/1‑250 (MQ=25)                                                                                                               aaTCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACACACCGAAAGCAACGAAAACGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATAGCTAAGGTGGAAGCGGACGTAAAACGTGAAGGATAACGTTGCCTGAGCAACGGCCCGACGGGAGAGGAGAAGTCGAG‑GAATCCTGCACGAACCACTCCTAACCtagtt                                                                                                                              >  2:612859/1‑249 (MQ=25)                                                                                                                     aaaTATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTgtag                                                                                                                      <  2:460035/250‑1 (MQ=255)                                                                                                                       atatatGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGACGTCGAG‑TCATCCTGCACGACCCACCACTAACATAATTAGTTgtagt                                                                                                                     >  2:87472/1‑250 (MQ=255)                                                                                                                               cgccgTTGCCAGGGGATATCAAACAAACCGAAAGCAGCGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTGCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCCCCACTAACATAGTTAGTTGTAGTATCCAGCgt                                                                                                            <  1:11752/251‑1 (MQ=255)                                                                                                                               cgccgTTGCCACGGGATATCAACCAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATAGCTAAAGCGGAAGCGGAGGTAAAGCGAGAAGGATAACGTTGCATAAGCAACGGCCCGAAGGGCGAGGCGAAGTCGGG‑TCATCCTGCACGACCCAACACTAAAATAGTTAGTAGTAGTATCCAGCg                                                                                                             >  1:157680/1‑250 (MQ=25)                                                                                                                               cgccgTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCgt                                                                                                            >  1:595187/1‑251 (MQ=255)                                                                                                                                     tGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTTGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAGGTAGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTAt                                                                                                       >  1:134691/1‑250 (MQ=255)                                                                                                                                           gggATATCAAACAAACCTAAAGCAACGAAAAAGTGGGGCGTTAGCTCAGTGGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAAACGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCGGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGa                                                                                                <  2:438131/251‑1 (MQ=255)                                                                                                                                              atatCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCAAGACTCACCAATCGCTAAGGTTGAAGCGGTAGTAAAACGGTAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTa                                                                                             <  1:343982/251‑1 (MQ=255)                                                                                                                                                  caaacaaaCCGAAAGCAACGAAAAAGTGGTTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCTCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGTGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTGGTAGTATCCAGCGTAGTATCGGGTGATTAGCTc                                                                                         <  1:319388/251‑1 (MQ=255)                                                                                                                                                   aaacaaacCGAAAGCAACGAAAAAGTGGGTAGTTAGCTCAGTTAGAACAGACGTTGACTTTTAATCAATTGGTAGCAGGTTCGAATCCTGCACGACTCAACACTCGCTAAAGTGGAAGCGGTAGTAAAAAGTGACGGATAACGTTGCATGAGCAACGACCCGAAGGGCGAGACGAAGTCGCG‑TAATCCTGAACGACCACCCCCTAACAAAGTTAGTTGTCGTAACCAGCGTAGCATCGGGTGATTAGCTc                                                                                         >  1:577916/1‑250 (MQ=255)                                                                                                                                                    aacaaacCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTc                                                                                         >  1:446270/1‑249 (MQ=255)                                                                                                                                                                        aaaGTGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCCGGTTCGGGTACTGCACGACACACAGGTTACTTCTGGGGAGTCGGCAGTAAAACGTGTAGGATAAAGTTGCATAAGCAACGGCCCGAAGGGCGAGACGAAGGAGAG‑ACATCCTGCACGACCCACCACTAGCATAGTTGGTTGTAGTATGCAGCGTATTATCTGGCGAATAGCttaggtggaatagcacatccgta                                                                    >  2:73284/1‑228 (MQ=12)                                                                                                                                                                         aaggggggcGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACGCCCCAATCGTTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTAc                                                                  <  1:211492/243‑1 (MQ=255)                                                                                                                                                                          aGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTAc                                                                  >  1:74731/1‑250 (MQ=255)                                                                                                                                                                          aGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACa                                                                 <  1:620891/251‑1 (MQ=255)                                                                                                                                                                             gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGGGTAGAGCGCGGCTCAGGAATCGCTTGGGGGGCAGGGGGGTTGCAGCGTGATGGAAAACGTTGCATGAGCAACGGCCCGACGGGCGAGCCGAGGTCGAG‑GCATCCTGCACGACGCGACACTAACATAGTTAGTTGTAGTATCACGCGTAGTATCGGGGGATTAGCTCAGTTGGGAGAGCACCTCCCTTACAagg                                                              >  2:106729/1‑251 (MQ=255)                                                                                                                                                                                         cAGTTGGTAGCGAAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCACAAAGGCGTAAGCGGTAGTAAAACGTGCAGGATAACGTTGCATGAGCAACGGCCCGAAGGGAGAGACGACGTCGAGATC‑TCCTACAAGCCCCCCCCCTAACAGAGGTAGTTGGAGTATCTAGCGTGTTATCGGATGATTAGCTCAGCTGGAAGAGCTCCTCCCATACAaggcgcgcgtcatcg                                                  >  2:605975/1‑239 (MQ=16)                                                                                                                                                                                                       gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCTCTAAGGTGGAGGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGTCGAGCCGAAGTCGAG‑TCATCCTACACGACCCCCCTCTAACATAGTTCGTTGTAGTATCCAGCGTAGTATTAGGACATTAGCTCAGCGGGGGGAGCACTTCCATTACAAGGAGGGGGTAGGCGGTTCGagcacggct                                    >  2:170791/1‑243 (MQ=21)                                                                                                                                                                                                                 aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGCCGACGTCGAG‑TCATCCTGCACGACACACCACTCACATAGTTAGCTGTAGCATCAAGCGTAGTATCGGGCGATTAGCTCAGATGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGGTCCCGCAATcaaccaccc                          >  1:268288/1‑250 (MQ=37)                                                                                                                                                                                                                    cAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCg                       >  1:221044/1‑251 (MQ=255)                                                                                                                                                                                                                                    cGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGtt       <  1:512995/251‑1 (MQ=255)                                                                                                                                                                                                                                    cGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGtt       >  1:77660/1‑251 (MQ=255)                                                                                                                                                                                                                                         ccTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTag  <  1:159812/251‑1 (MQ=255)                                                                                                                                                                                                                                         ccTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGGCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTAGCGGGTGAGTAGCTCAGCTGGGAGAGCACCTCCCTGACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTa   <  2:595187/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  CGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAG‑TCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAG  >  NC_000913/780382‑780864

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