Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,728,118	(C)8→7	coding (1254/1455 nt)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,118	0	C	.	98.3%	242.1 / ‑2.2	58	coding (1254/1455 nt)	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (36/21);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGC  >  NC_000913/3727868‑3728368                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           gaTATATACATGCGAATTATCCCCCACCCGGTCCGGAAGGGGATAACGTTTACGACATTAATGGCAGAAGTTGCTGAAAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTAAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGCGCTGGCACCACATcccccccc                                                                                                                                                                                                                                                            >  2:342517/1‑251 (MQ=255)     tataGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCACACATCACCCCC‑GTACg                                                                                                                                                                                                                                                       >  1:177013/1‑251 (MQ=255)        aGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGt                                                                                                                                                                                                                                                     >  2:139566/1‑250 (MQ=255)              gAATTATCCCCCCCCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGGAGAAGTTGCTGGTAGAGGCGACGGAACGTGTCTCGTCGTGGCTGAGAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTATTTGCTGCAACAATTCAATGAGCGATTTCTCTGTATTCGCCGCGAGCTGCGCCAGCCTTGCGGCGCCCAGGGCTGGCCCCTCATCCCCCCC‑TTCCGGTAATCGa                                                                                                                                                                                                                                               <  2:445013/250‑1 (MQ=255)                         cACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTACTCGAGCTGCTGACCGc                                                                                                                                                                                                                                   >  1:480923/1‑251 (MQ=255)                          aCCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                                                                                                                  >  2:595238/1‑251 (MQ=255)                           cccGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                                                                                                                                 >  1:507440/1‑251 (MQ=255)                            ccGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTCCGGTAATCGAGCTGCTGACCGCTg                                                                                                                                                                                                                                 >  1:577160/1‑250 (MQ=255)                             cGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCTGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCTCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGat                                                                                                                                                                                                                               <  2:499187/251‑1 (MQ=255)                               gTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCACATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGatat                                                                                                                                                                                                                             <  1:436065/251‑1 (MQ=255)                                 caggcaggGGATAACTGTTACGCCATTAATGGCAGAAGTTGCTGATAGAGTCGACGGAACGTTTCTCGTCGTGGCTGATAATCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                                                                                                                           <  2:240595/251‑1 (MQ=255)                                                 tttACGCCATTAATGGCAGCAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGcc                                                                                                                                                                                                           >  2:475084/1‑251 (MQ=255)                                                  ttACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGcc                                                                                                                                                                                                           <  2:98431/250‑1 (MQ=255)                                                    aCGCCATTAATGGCAGAAGTTGCTGATAGCGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGt                                                                                                                                                                                                        >  1:247630/1‑251 (MQ=255)                                                        cATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGAGGCTGATAAGCGGCATAAAGCCGCGCATCTGGGAGATGCGACTGTTCTCACGGGAGGTGCGGCAACAATTCAATGGGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGAGCTGGCCCCACATCCCCCCC‑GAACGGGAATCGAGCTGCTGACCGCTGCTATCCGCCAGAGTCTGACGCCAGTACTc                                                                                                                                                                                                    >  1:73913/1‑251 (MQ=255)                                                           tAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGCCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTACTCACt                                                                                                                                                                                                 >  2:83517/1‑251 (MQ=255)                                                                   gAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                                                                                                                                                         <  1:398789/251‑1 (MQ=255)                                                                    aaGTTGCCGATAGAGGCGACCGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCC‑GTACGGTAATCGAGCTGCCGACCGCCGATTTCCGCCCGCATCTGACGCCCGTACTCACTACCcgcccc                                                                                                                                                                                        >  2:410177/1‑251 (MQ=255)                                                                        tGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGcc                                                                                                                                                                                                                        <  2:578104/215‑1 (MQ=255)                                                                        tGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGcc                                                                                                                                                                                                                        >  1:578104/1‑215 (MQ=255)                                                                                     gACGGACCGTTTCTAGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCCACAATTAAAAGAGCGATTTCTCTGGATTCGCCCCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGCCCGCTGATATCCGCCAGCATCTGACGCCAGAACTCCCTACGCCCCCCGCCCCCACTCAACgca                                                                                                                                                                        >  2:471387/1‑248 (MQ=255)                                                                                        ggAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAacac                                                                                                                                                                    >  2:496143/1‑251 (MQ=255)                                                                                         gAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACt                                                                                                                                                                   <  1:595238/251‑1 (MQ=255)                                                                                             gTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCCTCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCCCCCCCCCCCCCATCAACGTAACACTTTGc                                                                                                                                                               >  2:670221/1‑251 (MQ=255)                                                                                                      cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACCCTTTGCGGTTTAAt                                                                                                                                                       >  2:186527/1‑250 (MQ=255)                                                                                                         ggCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTCCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAAACCg                                                                                                                                                   >  2:152411/1‑251 (MQ=255)                                                                                                             gATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg                                                                                                                                                <  1:414826/250‑1 (MQ=255)                                                                                                              aTAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAgc                                                                                                                                                                                                                                             >  1:581862/1‑156 (MQ=255)                                                                                                              aTAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAgc                                                                                                                                                                                                                                             <  2:581862/156‑1 (MQ=255)                                                                                                                  gcgccaaaaCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTATCTGGATCCGCCGCGACCTGCGCCAGCTTTGCTGCGCCCAGTGCTGGCCTAACA‑CCCCCCCCGTCCGGTAATCGAGCTTCTGACCGCTGATATCCGCCAGCATCTGATGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                                                                                                                          <  2:177013/244‑1 (MQ=255)                                                                                                                  gaggtataacgaggcgcAGCTGGTAGTGGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                                                                                                                          <  2:577160/238‑1 (MQ=255)                                                                                                                  gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                                                                                                                          <  1:469762/251‑1 (MQ=255)                                                                                                                            gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCAGCACCAATTCAATGAGCGATTTCTCTGGACTCGCCGCGATATGCGCCAGCCTTGCTGCGCCCAGCGCTGTCCCCACATCCCCCCC‑GTACGGCAATCGAGCTGCTGACCGATGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCACGCATGCACGACATcc                                                                                                                                 >  1:298373/1‑250 (MQ=255)                                                                                                                              tgagcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCCGCCTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAGTCAACGTAACACTTTGCGGTTTCATACCCCAGGCATGcccaacaaccctg                                                                                                                              >  2:32118/4‑239 (MQ=255)                                                                                                                               gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                                                                                                              <  1:218373/250‑1 (MQ=255)                                                                                                                                           gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg                                                                                                                 >  2:336259/1‑251 (MQ=255)                                                                                                                                               cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCAt                                                                                                              >  2:597402/1‑250 (MQ=255)                                                                                                                                               cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATa                                                                                                             <  1:139566/251‑1 (MQ=255)                                                                                                                                                       cTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTCCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTTTAATACCGCAGGCCAGCACGACATCCATGccctccacccgctcctaacccacgc                                                                                                     >  2:157826/1‑228 (MQ=255)                                                                                                                                                                gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGc                                                                                             >  1:311201/1‑250 (MQ=255)                                                                                                                                                                                      cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTggg                                                                      >  2:534421/1‑251 (MQ=255)                                                                                                                                                                                       gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTgggg                                                                     <  1:199371/251‑1 (MQ=255)                                                                                                                                                                                              ctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCCACGCATAACCCACGCCTTCCAGCACTGCCCGCGCCAGTTCATTGGGGCCATGtt                                                              >  2:548237/1‑251 (MQ=255)                                                                                                                                                                                                     cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGAGATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCTCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAg                                                        <  1:362243/250‑1 (MQ=255)                                                                                                                                                                                                             tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCCCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGACCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCa                                                >  2:367230/1‑250 (MQ=255)                                                                                                                                                                                                                  gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGaaa                                          <  1:186527/251‑1 (MQ=255)                                                                                                                                                                                                                   ccAGCATTGCTACGCCCAGTGCTGTCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCAAGCCTCTGACGCCAGAACTCACTACGCGCCCCGCCCCCACTCACCGTAACACTTTGCGGTTTAAAACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATACCCCACGCCTTCCAGCACTGCTCGCGCCCGTTCATTGGGGCCATGTTGATGAGTCAAACCAaaaaaa                                          >  2:363255/1‑250 (MQ=255)                                                                                                                                                                                                                          tgctgcGCCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCtt                                  >  1:480103/1‑251 (MQ=255)                                                                                                                                                                                                                             tgcGCCCAGACCTGGACCCACATACCACCC‑GTAAGGTAATCGCGCTGCTGACCGCTGATATCCGCAAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCATGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCCGTTCATTGGGGCCCTGTTGATGCGTCAAAACCAAGAAAACCCCCTTCGc                               >  2:118661/1‑251 (MQ=255)                                                                                                                                                                                                                               cgcCCAGTGCTGGCCACACATCCCCACC‑GTACGGTAATCGAGCTGCTGACCGCTGCTACACGACAGCATCTGACGCCAGCACACACTCCGCGCCCCGCACCCAATCAACGAAACACTCTGCGGTTTAATACCGCAGGCATGCCCGACCTCCATGCCATCTGCCAGCGCATAACCCACGCCTTACCGCACTGCTCGCGCGAGATCATTGGAGCCATGTtgat                                                           >  2:640026/1‑221 (MQ=255)                                                                                                                                                                                                                               cgcCCAGGGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCAACTTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtgat                                                           <  1:640026/221‑1 (MQ=255)                                                                                                                                                                                                                                gcCCAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGCCGCCAGTACTCAATACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTg                            >  1:71482/1‑251 (MQ=255)                                                                                                                                                                                                                                  ccAGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGGACTCAAGCCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCGCGCGCCCGTTCATTGGGGTTATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGGCTggg                          <  2:480923/251‑1 (MQ=255)                                                                                                                                                                                                                                    aGTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTgggg                         >  1:41441/1‑250 (MQ=255)                                                                                                                                                                                                                                     gTGCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGTGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGat                       >  1:366515/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCTCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATACCCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGACCATGTTGATGAGTAAAACCCAAGAAAACCCCCTTCGCCTGGGGat                       >  2:404412/1‑249 (MQ=255)                                                                                                                                                                                                                                       gCTGGCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGat                       <  1:404412/249‑1 (MQ=255)                                                                                                                                                                                                                                          ggCCCCACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTAATACCGAAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGATCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtg                  >  1:88265/1‑251 (MQ=255)                                                                                                                                                                                                                                              ccACATCCCCCCC‑GTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtgtg              <  1:90300/251‑1 (MQ=255)                                                                                                                                                                                                                                                      ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACg       <  1:496143/251‑1 (MQ=255)                                                                                                                                                                                                                                                           cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGc  >  1:117701/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           GATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGC  >  NC_000913/3727868‑3728368

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