Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,756,976	G→A	R273Q (CGA→CAA)	intA →	integrase IntA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,756,976	0	G	A	96.1%	176.1 / NA	51	R273Q (CGA→CAA)	intA	integrase IntA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/20);  minor base C (2/0);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGCTGTTCCTAGTGTAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGATTTCCACCT  >  NC_000913/2756731‑2757220                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     tgatGCTGTTCCTAGTGTAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATg                                                                                                                                                                                                                                                  >  1:536960/1‑250 (MQ=255)               gtgtAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGAc                                                                                                                                                                                                                                   >  1:254801/1‑251 (MQ=255)                 gtAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGAcc                                                                                                                                                                                                                                  >  1:378819/1‑250 (MQ=255)                         cGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTc                                                                                                                                                                                                                         >  1:23177/1‑251 (MQ=255)                         cGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTc                                                                                                                                                                                                                         >  2:259279/1‑251 (MQ=255)                            aaaaGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTACAATGAAGATGAACCGGGACCATACGGTTCCAc                                                                                                                                                                                                                      >  2:259115/1‑251 (MQ=255)                                      gagaAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATgg                                                                                                                                                                                                            >  1:549574/1‑251 (MQ=255)                                       agaAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATggg                                                                                                                                                                                                           >  2:7773/1‑251 (MQ=255)                                                  aaaaaaaaCATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTAt                                                                                                                                                                                                <  2:349617/248‑1 (MQ=255)                                                     aaaaaCATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGaa                                                                                                                                                                                                                                                >  2:403443/1‑200 (MQ=255)                                                     aaaaaCATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGaa                                                                                                                                                                                                                                                <  1:403443/200‑1 (MQ=255)                                                      aaaaCATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTg                                                                                                                                                                                            <  1:74536/251‑1 (MQ=255)                                                      aaaaCATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTAGCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTg                                                                                                                                                                                            <  2:18342/251‑1 (MQ=255)                                                           aTGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAat                                                                                                                                                                                       <  1:77449/251‑1 (MQ=255)                                                             gCCAAGCATCCGGCCGGATCAACTGCCGCAGCTAATGCACACCAGGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAatg                                                                                                                                                                                      <  1:135812/250‑1 (MQ=255)                                                                aaGCACCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCGGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAAtgatga                                                                                                                                                                                  <  2:306948/251‑1 (MQ=255)                                                                             gATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAgtgg                                                                                                                                                                     >  1:547445/1‑251 (MQ=255)                                                                                     tCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCg                                                                                                                                                               <  1:244909/249‑1 (MQ=255)                                                                                       cGCAGCTAAGGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCAAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGaa                                                                                                                                                           <  2:500442/251‑1 (MQ=255)                                                                                                 gCACACCATGCGTACGGCAAGTATCAGCATGGCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCACCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTc                                                                                                                                                 <  2:36685/251‑1 (MQ=255)                                                                                                  caccccATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTcc                                                                                                                                                <  2:23177/247‑1 (MQ=255)                                                                                                  cacaCCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCACATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAACTGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTcc                                                                                                                                                >  2:159355/1‑251 (MQ=255)                                                                                                                     gtaacAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTGCGAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaac                                                                                                                             <  1:178369/247‑1 (MQ=255)                                                                                                                           gCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGCAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCTAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaacc                                                                                                                        >  1:441347/1‑250 (MQ=255)                                                                                                                            catgcccccaAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCGGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGAAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaaccaa                                                                                                                      <  2:230515/243‑1 (MQ=255)                                                                                                                                      agtggccgtttCATGTGGAAACTTCAACCCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATCCGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCa                                                                                                            <  2:464531/242‑1 (MQ=255)                                                                                                                                        aTGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCaa                                                                                                           >  2:144807/1‑250 (MQ=255)                                                                                                                                           cctgttactgtggaaactgcaacccATCACCCGCACTGCCGAGGCTGCTGAGGCCCGATGGGATGAAGCCGATTGCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGATCATACGGTTCCACTAGCTGATGGGGCTCTTGCTATTCGGGAAATGATGAAGCCTCTCAGTGGTGGCCGAAAAGTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAg                                                                                                       <  2:90901/228‑1 (MQ=255)                                                                                                                                              gttaaTGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTg                                                                                                     >  2:184586/5‑250 (MQ=255)                                                                                                                                               tacaaggggaaaCTTCTAACCACCCCCCGCCCTGCCGAAGCTGCTGAGGTGAGAGTGGATGAGATCGATTTCAATGCTAGCCAATGTAAAATTCATGCATCTCAACTGAAAATGAACCGGGACCATACGGTTCCCCTATTTGAGGGGGCGCTTGCTATTCTGGAAATGATGAAGCCTCGCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGa                                                                                                    <  2:229111/240‑1 (MQ=255)                                                                                                                                                 cATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAAt                                                                                                  <  1:385141/250‑1 (MQ=255)                                                                                                                                                  aTGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTCTTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGTATgc                                                                                                >  2:158664/1‑251 (MQ=255)                                                                                                                                                         aaCTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGCGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGCACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGCTGAAGCCTCTCATTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTc                                                                                         >  2:329748/1‑251 (MQ=255)                                                                                                                                                                                    aaGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGTCCGAGAATTTATCTTGCCTAGCCGTATCAAGCCCAACCAACCAATGAAGAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACtt                                                              <  1:330167/251‑1 (MQ=255)                                                                                                                                                                                    aaGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACtt                                                              >  2:28917/1‑251 (MQ=255)                                                                                                                                                                                            gAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTcaca                                                      <  1:401014/251‑1 (MQ=255)                                                                                                                                                                                                cTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGACTTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATCGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGc                                                  >  2:48319/1‑251 (MQ=255)                                                                                                                                                                                                    aTGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGc                                              >  1:411398/1‑251 (MQ=255)                                                                                                                                                                                                    aTGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGc                                              >  2:474645/1‑251 (MQ=255)                                                                                                                                                                                                                 gATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGTACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGACTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCCATGACTAGCCAACCAGTGAATGCAGCACTCAAGCATGCTGGCTTAGGAGGTGTCCTTGTTTCACACGGCTTGCGTTCTATCGcccg                                 >  2:499185/1‑249 (MQ=255)                                                                                                                                                                                                                 gATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCACCCAATGAATAGCAAAACAGTGAATGCAGCACTCACGCGTGCTGGCTTGGGAGGTGTACTTGTTTCACACGGTTTGCGTTCTATCGcc                                   >  2:297937/1‑249 (MQ=255)                                                                                                                                                                                                                  aTTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGt                                >  1:313822/1‑251 (MQ=255)                                                                                                                                                                                                                      cAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACgg                            >  1:234791/1‑251 (MQ=255)                                                                                                                                                                                                                        aTGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGc                           >  1:489280/1‑250 (MQ=255)                                                                                                                                                                                                                               cGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATg                   >  1:51747/1‑251 (MQ=255)                                                                                                                                                                                                                                     gAAAATTCATGCAGCTACAATGAAGATGAACCGGGACCATACGGTTCCCCTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTAAGTGGTGACCGAGCATTTATCTTTCCTAGACGTATCAAGCCCAACAAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGagga               >  2:295952/1‑249 (MQ=255)                                                                                                                                                                                                                                       aaaTTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGaggaagg            >  2:299136/1‑250 (MQ=255)                                                                                                                                                                                                                                        aaTTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGAt          >  1:427727/1‑251 (MQ=255)                                                                                                                                                                                                                                         aTTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGAtt         >  2:38014/1‑251 (MQ=255)                                                                                                                                                                                                                                           tCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGCGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCACTGAGGAAGGTTTTc       >  2:563851/1‑251 (MQ=255)                                                                                                                                                                                                                                                cAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGATTTCCACCt  <  1:299136/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     TGATGCTGTTCCTAGTGTAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGATTTCCACCT  >  NC_000913/2756731‑2757220

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