Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,729,440	2 bp→TT	noncoding (105‑106/160 nt)	xylZ ←	small RNA XylZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,729,440	0	C	T	98.4%	224.3 / NA	64	noncoding (106/160 nt)	xylZ	small RNA XylZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/35);  minor base G (0/1);  total (28/36)
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
*	NC_000913	3,729,441	0	C	T	98.4%	217.0 / NA	64	noncoding (105/160 nt)	xylZ	small RNA XylZ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/35);  total (28/36)
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

GAGCTTTCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGCCAT  >  NC_000913/3729197‑3729679                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                               gagCTTTCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTAtt                                                                                                                                                                                                                                          <  1:410562/251‑1 (MQ=255)  agCTTTCATTGCGCGATCCGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTAttt                                                                                                                                                                                                                                         >  1:380329/1‑251 (MQ=255)      ttCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTc                                                                                                                                                                                                                                      <  2:416305/250‑1 (MQ=255)              gcgATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATa                                                                                                                                                                                                                             <  1:431517/251‑1 (MQ=255)                 aTCAGTTGCCTGCAACCACTGTTCCGGGTCTTGTTCCGACCAGAGTAGATGCGGGAGCGACACGGTCAGCCTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAGCTTTTACGCCCGAGGTGCCAAGGTCTCTCCCGATATACATATAGATCGTTCCTTAAAAAAATGCCCGGAATCGCTACCGATAACCGGGCCAACGGACTGACCAGTTAGTTGTTATTTGTCGAACAGATAATg                                                                                                                                                                                                                          >  2:265957/1‑251 (MQ=255)                     gTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTTGATGCGGGCGCGAAACGGTCAGCTTTTCCGGTTGCGCAGCAACC‑TCCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGTGGTGTCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTAGCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGttt                                                                                                                                                                                                                      <  1:531643/251‑1 (MQ=255)                     gTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGGTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGttt                                                                                                                                                                                                                      <  1:315928/251‑1 (MQ=255)                        gcctcaactcacTGTTCCGGGTCTTGTTCAGACCAGAGGGGATGCGGGCGCGCAACGGCCAGTTTTGCCGTTTGCGCAGTAACCCACGT‑GCCCTGCTCGTTGAGCAAAATAACTGTTACGGCCGAGTTGCCAAGATCTATCCGGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACTGGGCCAACGTACTGCACAGTTTGTTGTTATTTGTCGAACAGATAATGGCTTAcc                                                                                                                                                                                                                   <  1:129215/242‑1 (MQ=255)                            gccACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAg                                                                                                                                                                                                                 <  1:12587/249‑1 (MQ=255)                              caccacTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGAttt                                                                                                                                                                                                             <  2:322716/251‑1 (MQ=255)                                ccacTGGTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTc                                                                                                                                                                                                           <  2:464321/251‑1 (MQ=255)                                        ccGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑CCCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAGCAGATAATGGTTTTCCAGATTTTCCAgttgtt                                                                                                                                                                                                   <  1:136702/251‑1 (MQ=255)                                        ccGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATACCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttgtt                                                                                                                                                                                                   >  1:328994/1‑251 (MQ=255)                                                   tCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAcc                                                                                                                                                                                         >  1:207873/1‑250 (MQ=255)                                                         cAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTg                                                                                                                                                                                   <  2:553451/250‑1 (MQ=255)                                                         cAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTACTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGCCCACTCTg                                                                                                                                                                                   >  1:553451/1‑250 (MQ=255)                                                             gTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCa                                                                                                                                                                              <  1:369787/251‑1 (MQ=255)                                                                   gCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGgag                                                                                                                                                                        <  2:30554/251‑1 (MQ=255)                                                                             aaCGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAatgatg                                                                                                                                                              >  1:302584/1‑251 (MQ=255)                                                                              aCGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGt                                                                                                                                                             <  1:211335/251‑1 (MQ=255)                                                                               cGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTAACAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGtt                                                                                                                                                            <  1:424472/251‑1 (MQ=255)                                                                                   cAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTg                                                                                                                                                        >  2:157997/1‑251 (MQ=255)                                                                                     gCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAg                                                                                                                                                      <  1:526522/251‑1 (MQ=255)                                                                                        tttCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCat                                                                                                                                                   >  1:51174/1‑251 (MQ=255)                                                                                             gTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAACATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTg                                                                                                                                              >  1:510584/1‑251 (MQ=255)                                                                                                tGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGc                                                                                                                                            <  1:83489/250‑1 (MQ=255)                                                                                                tGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGc                                                                                                                                            <  1:52159/250‑1 (MQ=255)                                                                                                tGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGAGCCAAGATCTATCCCGATATACATATCGATCGCTCCTTAAAAAAATGCCCGGTATCGCTACCGATAAACGGGCCAAAGGACTGAACAGCGAGTTGTTATTTGTCGAATAGAAAATGGATTACCAGATTTTCCAGTTGTTCCTGGCGACAACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCCTAGTTGGCt                                                                                                                                           >  1:45432/1‑251 (MQ=255)                                                                                                 gcgcAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATACCCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCTCTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCt                                                                                                                                           <  1:333408/250‑1 (MQ=255)                                                                                                     agcaACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATc                                                                                                                                      >  1:246955/1‑251 (MQ=255)                                                                                                        aaCC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGc                                                                                                                                   >  2:127737/1‑251 (MQ=255)                                                                                                                 cACCCGGCTCGTTGAGCAAAATAACTTGTACGGCCGAGGTGCCAAGAGCTAACCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTTCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCGGGCGACCACTCTGATGCACCGGAGACAAAGGATGTTCCTGAGCATATTTGGCTAAAGGTGCCAGTGACa                                                                                                                           <  2:452305/251‑1 (MQ=255)                                                                                                                              gagaaaaatatctgtgaCGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGTTATCGCTACCGATAGCCGGGCCAACGGACTGCACAGTGAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGGTTTTCCAGCTGTTCCTGGCGACCGCTCTGATGCAACTGAGGCAAATGATGGTCCTGAGCATACTGGGCGAAATCTGCCAGCGACATTTGGCATTTCa                                                                                                               <  2:265908/234‑1 (MQ=255)                                                                                                                               aGCAAAATAACTTTTACGCCCGAGGGGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGTCATTTGGCCTTTCAgg                                                                                                             <  1:454955/251‑1 (MQ=255)                                                                                                                                           tttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGccc                                                                                                 >  2:75970/1‑251 (MQ=255)                                                                                                                                               cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAAtt                                                                                             >  1:392377/1‑251 (MQ=255)                                                                                                                                                    gAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAACCACATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGCGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTAGCt                                                                                         >  1:258642/1‑250 (MQ=255)                                                                                                                                                    gAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCt                                                                                         >  1:156516/1‑250 (MQ=255)                                                                                                                                                            aaGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCTCAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGcc                                                                                <  1:127737/251‑1 (MQ=255)                                                                                                                                                                   aTCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGa                                                                             <  1:98546/247‑1 (MQ=255)                                                                                                                                                                      ccGATATACATATCGACCGTTCCTTACAAAAACGCCCGGTATCGCTACCGATAACCGGGCCAACGGGCTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGTCTACAACTCTGCTGCCCTGGAGACAAATGATGTTCATGAGCATATTTGTCTACATCTGCCAGTGACATTTGGCCTCTCAGAATTTGCTGGACCAATTCGCTATTCCAGACGGAGTAACGc                                                                      >  2:257962/1‑251 (MQ=255)                                                                                                                                                                       cGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTAGTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCATCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCt                                                                     <  1:325740/251‑1 (MQ=255)                                                                                                                                                                        gATATGCATATCGATCATTCCATAAAAAAATGCCCGATATCGCTACCGATAACCGAGCCAACGGACTGCCCAGTTCATTGTTATTCGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGCGCATATTTGGCTACATCTGCCAGCGACATTTGGCCTTTCAGGATTTGCTGGACCAATCCGCTATTCCAGCCTGAATAAAGCtt                                                                    >  2:184533/1‑250 (MQ=255)                                                                                                                                                                              cagatcgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGAt                                                              <  2:11159/248‑1 (MQ=255)                                                                                                                                                                                tatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGc                                                            >  1:397896/1‑251 (MQ=255)                                                                                                                                                                                 atcgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGc                                                            <  2:397896/250‑1 (MQ=255)                                                                                                                                                                                 aacgatcgTCCCGTAAAAAAAGGCCAGGAAGCGCTACCGATAACCGGGCCAACGGGATGCACAGGTAGGCGATATTTGTCGAACTGATAATGGTTTTCCAGATTCTCAAGTTGTTCCTGGGGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGGATATTTGGATAAATCTGCCAGAGACATCTGGACGTTCAGGATTTGATGGCCTAATACGCTAATCCAGCCGCAACTACGCTGCGCGATGCg                                                           <  2:467727/249‑1 (MQ=255)                                                                                                                                                                                             ttAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCa                                                                                                               <  1:78032/187‑1 (MQ=255)                                                                                                                                                                                             ttAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCa                                                                                                               >  2:78032/1‑187 (MQ=255)                                                                                                                                                                                                    aaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCt                                        >  2:8454/1‑251 (MQ=255)                                                                                                                                                                                                    aaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCAAGCTCGCCATCt                                        >  1:112660/1‑251 (MQ=255)                                                                                                                                                                                                       gCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGGGGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGTGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGTCCCAATTCGCTATTCCAGCCTGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTc                                      <  2:520818/250‑1 (MQ=255)                                                                                                                                                                                                           ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCt                                                                                                      >  2:341403/1‑182 (MQ=255)                                                                                                                                                                                                           ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCt                                                                                                      <  1:341403/182‑1 (MQ=255)                                                                                                                                                                                                           ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTAGTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTGCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCAGATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTTCTGGCCCAATTCTCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTGCAATCa                                 <  2:493658/251‑1 (MQ=255)                                                                                                                                                                                                             tATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATg                               >  1:267444/1‑251 (MQ=255)                                                                                                                                                                                                              aTCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGGCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCGTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTGATCCAGCTCGCCATCTCCAATCATg                               <  2:510584/250‑1 (MQ=255)                                                                                                                                                                                                                       gATAACCGGGCCAACGGACTGCAAAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCa                      >  2:115857/1‑250 (MQ=255)                                                                                                                                                                                                                         tAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAAtt                   >  1:138082/1‑251 (MQ=255)                                                                                                                                                                                                                              gggCCCACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTc                >  2:263044/1‑249 (MQ=255)                                                                                                                                                                                                                                       gACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGc     <  2:380329/251‑1 (MQ=255)                                                                                                                                                                                                                                       gACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGc     >  1:312735/1‑251 (MQ=255)                                                                                                                                                                                                                                       gACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGc     <  1:35430/251‑1 (MQ=255)                                                                                                                                                                                                                                           gCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGCCAt  <  1:152317/250‑1 (MQ=255)                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                               GAGCTTTCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACC‑ACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGCCAT  >  NC_000913/3729197‑3729679

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