Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,044,642	T→C	intergenic (+18/‑296)	asnT → / → yeeJ	tRNA‑Asn/inverse autotransporter adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,044,642	0	T	C	92.7%	110.1 / 1.1	41	intergenic (+18/‑296)	asnT/yeeJ	tRNA‑Asn/inverse autotransporter adhesin
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (17/21);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-02

CATCACGCTAATGATACAGACTCGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATTTACATTTTT  >  NC_000913/2044401‑2044885                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    cATCACGCTAAATAAATAGACTCGTTTTCGGCGAAGAAAGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACGTATGAAATTTAGTGTTGACATACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTTTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAgcgc                                                                                                                                                                                                                                            <  2:131628/251‑1 (MQ=25)          aattaTACAGACTCGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCa                                                                                                                                                                                                                                     <  2:99848/245‑1 (MQ=255)          aaTGATACAGACTCGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCa                                                                                                                                                                                                                                     <  2:509189/249‑1 (MQ=255)              aTACAGACTCGTTTTCGGCGACGAATTTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAATGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGACCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATGCGCTCTTTTAGCGCAATGTCACTGAcc                                                                                                                                                                                                                               <  2:485235/251‑1 (MQ=255)                acagacctgTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATGAACCAATTGAAATGACTTACGCAATCTATTTTTGACAGACAAGGTACCGCTAAGTAATAGTCATCCCGTTCACACTATTCAGCTGTAGTTCATTCAGTAGACCGGTGGCCTCTTAATTCGTATTTCACTTGTTCGAGCCCATTGATATGAGCCAAATTCCAAAAAGCCGCTCTTCTATCGCAATGTCACTgtcttt                                                                                                                                                                                                                             <  2:509692/243‑6 (MQ=255)                 cAGACTCGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTa                                                                                                                                                                                                                            <  1:183439/251‑1 (MQ=255)                       cGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGa                                                                                                                                                                                                                       <  2:481308/250‑1 (MQ=255)                            tCGGCGAAGAATGTTCATTAATTAACATCTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGAAAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCGGTAGTTCAGTCGGTAGAACGGCGGACTGTGAATCCGTATGTCACTGGTTCGAGTCCAGTCAGCGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTATTTGAACAtt                                                                                                                                                                                                                  <  2:319328/250‑1 (MQ=37)                               gCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGTTAGAACGGAGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGttt                                                                                                                                                                                                              <  1:178437/251‑1 (MQ=255)                                   cgaATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTa                                                                                                                                                                                                          <  1:26915/251‑1 (MQ=255)                                         ttCATTAATTAACAACTTTGCAGATTAGTTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCCGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCACTGTCACTGCCCTTATTTGAACCTTGTTTTTTAACGGAt                                                                                                                                                                                                    >  2:100958/1‑251 (MQ=37)                                                 ttaaCAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGtt                                                                                                                                                                                            >  1:118611/1‑251 (MQ=255)                                                   aacaacTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCTCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGtttt                                                                                                                                                                                          <  2:466409/251‑1 (MQ=255)                                                    ccaacTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGAACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTAGTTCGCGTCCAGTCAGGGGAGCCAAATTCTAAAAATTCGCCCTTTTAGCGCAATGTCACTGCCCTTAGTTGAACATTGTTTTTTAACGGATAACGGGttttt                                                                                                                                                                                         >  2:407150/2‑251 (MQ=37)                                                           tGCAGATTAATTCACCAATTGAAATCACTTATGAAATTTAGTGTTGACAGCCAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGAATGTTAATCCGTATGTCACTGGTTCGAGTCCAGACAGAGAAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCAGGTTTTTAACATCt                                                                                                                                                                                  >  1:444505/1‑251 (MQ=37)                                                            gCAGATTAATTAACCAATTTAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTATAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGTAGCCAAATTCTAAAAATGCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCt                                                                                                                                                                                  <  1:323713/250‑1 (MQ=255)                                                                  taattaaCCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAgcgc                                                                                                                                                                           >  1:278332/1‑251 (MQ=255)                                                                                      cTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAAc                                                                                                                                                                                                                     <  1:335500/189‑1 (MQ=255)                                                                                      cTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAAc                                                                                                                                                                                                                     >  2:335500/1‑189 (MQ=255)                                                                                              aTTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGtttt                                                                                                                                               >  1:138881/1‑251 (MQ=255)                                                                                                                    aCCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAtt                                                                                                                          <  1:455016/250‑1 (MQ=255)                                                                                                                                  ttCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATtat                                                                                                            >  1:27613/1‑250 (MQ=255)                                                                                                                                       cccGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACCTCTTAAGCGCCCTCGACCTTTATGGTCGAGGGCGTTTTGCTATGAACCCCATCACCCTTTTCCCCTCGATTATAAAACt                                                                                                      >  2:432659/1‑251 (MQ=255)                                                                                                                                         cGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTg                                                                                                    >  1:511590/1‑251 (MQ=255)                                                                                                                                                               ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCCCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGAGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGATATGAACGCCCTCACCATTTTCCCCTCGCTTATAAAACTTGAGTTATTCAGTAGTCTCCCctc                                                                              >  2:55865/1‑251 (MQ=255)                                                                                                                                                                     cGGTAGAACGGCGGACTGATAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGTGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTACCCTCTTGCa                                                                         >  1:199744/1‑250 (MQ=255)                                                                                                                                                                                cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCa                                                              >  2:388462/1‑250 (MQ=255)                                                                                                                                                                                 ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGCTAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGtctc                                                                                   <  2:100199/228‑1 (MQ=255)                                                                                                                                                                                 ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGtctc                                                                                   >  1:100199/1‑228 (MQ=255)                                                                                                                                                                                 ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCtt                                                                            >  2:445592/1‑235 (MQ=255)                                                                                                                                                                                 ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCtt                                                                            <  1:445592/235‑1 (MQ=255)                                                                                                                                                                                      gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGCCCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACt                                                         <  2:134061/249‑1 (MQ=255)                                                                                                                                                                                                    cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGtt                                          <  1:55865/250‑1 (MQ=255)                                                                                                                                                                                                            cGAGTCCAGTCAGAGTTGACAAATTCTAAAAAGTCGCTCGTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTCCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCGTGCAACTCGCACCCAAAACTGCCTAACGAAAAGTTATTAAttt                                  <  2:74086/250‑1 (MQ=255)                                                                                                                                                                                                              aGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCTCAATGTCACTGACCTGAGTTGAATATTGTTTTTTAACGGACTGCGGGTTTTTAACATCTAAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAAtttt                                 <  2:180632/249‑1 (MQ=255)                                                                                                                                                                                                                   aGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCCACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAAtaat                          >  2:400156/1‑248 (MQ=255)                                                                                                                                                                                                                         aggagCCAAATTCTAAAAAGTCGCTCTTTTAGCGCAATGGCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCACCACCATTTTCCCCGCGATTATAAAACTTGAGTTATTCAATAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCt                    <  1:211154/251‑1 (MQ=255)                                                                                                                                                                                                                          ggagCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGGTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTa                   <  2:251440/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaaaTTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAGCTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATTTACAtt     >  2:520763/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaaaTTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAAAGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATAtt                       >  1:545064/1‑233 (MQ=255)                                                                                                                                                                                                                                           aTTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATTTACAttttt  >  1:92107/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    CATCACGCTAATGATACAGACTCGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATTTACATTTTT  >  NC_000913/2044401‑2044885

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