Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,729,441	C→T	noncoding (105/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,441	0	C	T	80.0%	92.8 / 21.6	45	noncoding (105/160 nt)	xylZ	small RNA XylZ
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (17/19);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CGCCCAGAGCTTTCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGC  >  NC_000913/3729191‑3729676                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            cGCCCAGAGCTTTCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCCCAGTTAGcc                                                                                                                                                                                                                                             <  1:402423/251‑1 (MQ=255)  gCCCAGAGCTTTCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCg                                                                                                                                                                                                                                            >  1:52632/1‑251 (MQ=255)            ttCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTc                                                                                                                                                                                                                                   <  2:107584/250‑1 (MQ=255)               aTTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTGCCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAAc                                                                                                                                                                                                                               <  1:227446/251‑1 (MQ=255)                    gcgATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGAt                                                                                                                                                                                                                           >  1:108200/1‑250 (MQ=255)                      gATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAAt                                                                                                                                                                                                                        >  1:291353/1‑251 (MQ=255)                                ctgccACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTACAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAcc                                                                                                                                                                                                                >  2:163791/1‑249 (MQ=255)                                       cacTGTTCCGGGTCTTTTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTcc                                                                                                                                                                                                       <  1:469185/251‑1 (MQ=255)                                        acTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCa                                                                                                                                                                                                      <  1:1366/251‑1 (MQ=255)                                           gTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgtt                                                                                                                                                                                                   >  2:49279/1‑251 (MQ=255)                                                                agagTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGAt                                                                                                                                                                              <  1:182407/251‑1 (MQ=255)                                                                   gTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCa                                                                                                                                                                           <  1:112611/251‑1 (MQ=255)                                                                                            cTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTACTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGAACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGCAGACAAATGATGTTCCTGAGc                                                                                                                                                  >  2:17329/1‑251 (MQ=255)                                                                                              tttCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCat                                                                                                                                                >  2:338265/1‑251 (MQ=255)                                                                                               ttCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCat                                                                                                                                                <  2:299064/250‑1 (MQ=255)                                                                                                  cGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATAttt                                                                                                                                            <  1:471497/251‑1 (MQ=255)                                                                                                             caACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTg                                                                                                                                 <  2:378650/251‑1 (MQ=255)                                                                                                              aaCCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGATGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGc                                                                                                                                <  1:309756/251‑1 (MQ=255)                                                                                                                  accTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGt                                                                                                                            >  1:118461/1‑251 (MQ=255)                                                                                                                            gCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAAGGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTgg                                                                                                                   <  1:204001/250‑1 (MQ=255)                                                                                                                                  tGAGCACAATAACTTTTACGCCCGAGGAGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCGCGGTATCGCTACGGATAAACGGGCCAACGGACTGCACAGTTAGTTGTTATGTGACGAACGGATAATGCGTTGCCAGATTTTCCAGTTGTTCTTAGCGACCACTCTGATGCATCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTTCAAGTGACATTAGCCCtttg                                                                                                             >  1:175419/1‑249 (MQ=255)                                                                                                                                   gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATCCATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGGCATTTGGCCTTTCAg                                                                                                           <  1:30052/251‑1 (MQ=255)                                                                                                                                   gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTc                                                                                                             >  2:413242/1‑249 (MQ=255)                                                                                                                                       aaaaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACAACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGAtt                                                                                                       >  1:268180/1‑251 (MQ=255)                                                                                                                                        aaaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGAttt                                                                                                      <  1:85435/251‑1 (MQ=255)                                                                                                                                            aaCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCt                                                                                                   <  1:191089/250‑1 (MQ=255)                                                                                                                                                  tACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAAGGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCa                                                                                             <  2:29782/250‑1 (MQ=255)                                                                                                                                                     gCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTc                                                                                         >  1:90459/1‑251 (MQ=255)                                                                                                                                                        cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCt                                                                                      >  1:136338/1‑251 (MQ=255)                                                                                                                                                        cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCt                                                                                      >  1:428067/1‑251 (MQ=255)                                                                                                                                                               ccAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAACTGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCATTTCAGGATTTGCTGGCCCAATTCGCTATTcc                                                                                 >  2:338442/1‑249 (MQ=255)                                                                                                                                                                 aaGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGcc                                                                             <  2:291353/251‑1 (MQ=255)                                                                                                                                                                            cGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCt                                                                  <  2:136338/251‑1 (MQ=255)                                                                                                                                                                              atatACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTgc                                                                <  1:194207/251‑1 (MQ=255)                                                                                                                                                                                   cATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGAt                                                           >  1:173652/1‑251 (MQ=255)                                                                                                                                                                                          atcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATATCCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTGGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTa                                                    <  1:257048/251‑1 (MQ=255)                                                                                                                                                                                                  ttAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTc                                            >  1:66022/1‑251 (MQ=255)                                                                                                                                                                                                     aaaaaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGcc                                         >  1:203345/1‑251 (MQ=255)                                                                                                                                                                                                        aaaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATc                                      <  1:371493/251‑1 (MQ=255)                                                                                                                                                                                                         aaTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCt                                     <  1:99195/251‑1 (MQ=255)                                                                                                                                                                                                          aTGCCCGGTATCGCTACCTATAACCGGTCCAACTGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGACTTTCCAGTTGTTCCTGGCTACCACTCTGATGCACCTGAGACAAATGATGTTCCTGTGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGTATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCtt                                    <  1:174448/251‑1 (MQ=255)                                                                                                                                                                                                            gCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTc                                   <  2:345891/250‑1 (MQ=255)                                                                                                                                                                                                                      gCTACCGATAACCGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATgcgcg                        >  1:287428/1‑251 (MQ=255)                                                                                                                                                                                                                         aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATCTTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTg                     >  2:175065/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGGGCCAACGGACTGCACAGTTAGTTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCa            <  2:108200/251‑1 (MQ=255)                                                                                                                                                                                                                                            gACTGCACAGTAAGGTGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGTTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCTCAATTCACTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAACCATGCGCTCTGCAATTTTCAGCGCCAGTGc  >  2:237126/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            CGCCCAGAGCTTTCATTGCGCGATCAGTTGCCTGCCACCACTGTTCCGGGTCTTGTTCCGACCAGAGTGGATGCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGGCTAAATCTGCCAGTGACATTTGGCCTTTCAGGATTTGCTGGCCCAATTCGCTATTCCAGCCGGAATAACGCTGCGCGATGCGTTTATCCAGCTCGCCATCTTCAATCATGCGCGCTGCAATTTTCAGCGCCAGTGC  >  NC_000913/3729191‑3729676

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