Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,985,816	T→C	intergenic (‑224/‑31)	yqeF ← / → yqeG	putative acyltransferase/putative transporter YqeG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,985,816	0	T	C	98.1%	194.3 / NA	53	intergenic (‑224/‑31)	yqeF/yqeG	putative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (34/18);  minor base G (0/1);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGT  >  NC_000913/2985570‑2986056                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 ccGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGtt                                                                                                                                                                                                                                              <  1:237084/251‑1 (MQ=255)         cacaACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCCATAGTACATGATTGTTAATTACAGTTATCTTCATTTTTTAAAAAAAGTGATTTTTATCACAACGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGAAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGt                                                                                                                                                                                                                                      >  2:406064/1‑251 (MQ=255)            aaCGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGt                                                                                                                                                                                                                                      >  2:112004/1‑248 (MQ=255)            aaCGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGGTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAga                                                                                                                                                                                                                                   <  1:328265/251‑1 (MQ=255)                   ttCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATAc                                                                                                                                                                                                                             <  2:161730/250‑1 (MQ=255)                    tcatcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACg                                                                                                                                                                                                                            <  1:315171/250‑1 (MQ=255)                           aaCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaa                                                                                                                                                                                                                    <  2:363801‑M1/251‑6 (MQ=255)                           aaCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaa                                                                                                                                                                                                                     >  1:132803‑M1/1‑246 (MQ=255)                            aCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTGAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTATATAGGGCGTTAATCTCCCATACTTAACCTGGTTGATGGTAAATTGCCCTCCATTGTGTTTAATGTTTGGATGAGACGTTCtaaaaa                                                                                                                                                                                                                   <  2:159114‑M1/251‑7 (MQ=255)                              cTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaa                                                                                                                                                                                                                  >  1:116749‑M1/1‑243 (MQ=255)                                cGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaacc                                                                                                                                                                                                               <  1:300182‑M1/251‑11 (MQ=255)                                cGCACCAATTAACAGTGCAATAATAAATTATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCCCTAACCCACCAATAGAGGTGTAGATAGGGCTTTAATCTCCCATACTTAACCTGGTTTATGGTAAAGTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaac                                                                                                                                                                                                                <  2:292230‑M1/250‑10 (MQ=255)                                 gCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGGTATATTAATTGTTTAAAAAAAGAGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGACAGAGACATAACTGGCACGTAGGGTTTGCACCCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaacaaccc                                                                                                                                                                                                              >  1:229120‑M1/1‑240 (MQ=255)                                    ccAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAACGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgc                                                                                                                                                                                                           >  2:409865‑M1/1‑237 (MQ=255)                                          aaCAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttg                                                                                                                                                                                                     >  1:32263‑M1/1‑231 (MQ=255)                                                     aGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACACAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAAACCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtg                                                                                                                                                                                          >  1:369975‑M1/1‑220 (MQ=255)                                                        aaaTGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgag                                                                                                                                                                                        <  2:117412‑M1/250‑34 (MQ=255)                                                                       aCAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCATAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGGTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgctta                                                                                                                                                                        <  2:224302‑M1/251‑50 (MQ=255)                                                                       aCAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGAAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCGGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatgttccggggttttgagcataacgtaaggctta                                                                                                                                                                        <  2:337691‑M1/251‑50 (MQ=255)                                                                               cTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccg                                                                                                                                                                >  1:308049‑M1/1‑194 (MQ=255)                                                                                  aaTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggac                                                                                                                                                             >  2:182486‑M1/1‑191 (MQ=255)                                                                                           aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccg                                                                                                                                                    >  2:85000‑M1/1‑182 (MQ=255)                                                                                               aaGTGATTTTTATCACAACGGACATATGCCTGAGCAGCAGACAGAGACATAACTGGCACGTAGGGTTTGCAACCACTACCCCCCCCATAGAGGGATAGATAGGACGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTACTTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggtgttgtgagcataacgtaatgcttattttaccggacgcatcgccgggg                                                                                                                                                 >  2:104470‑M1/1‑178 (MQ=255)                                                                                                  tGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACAAATAGCGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTAGGGTAAATTGCCCTCCATACTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtactgcttattttaccggacgcatcgccgggaaa                                                                                                                                               >  2:254806‑M1/1‑175 (MQ=255)                                                                                                           tCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGCTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtc                                                                                                                                    >  1:114081‑M1/1‑166 (MQ=255)                                                                                                                aaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcgg                                                                                                                                >  2:351631‑M1/1‑161 (MQ=255)                                                                                                                  ggAAATATGACTGAGCAGCAGTCCGCGACATAACTGGCACGTAGGGTTTGCAAACACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatg                                                                                                                             >  2:325586‑M1/1‑159 (MQ=255)                                                                                                                          gCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggct                                                                                                                     >  1:337819‑M1/1‑151 (MQ=255)                                                                                                                           ccTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctg                                                                                                                    >  2:29510‑M1/1‑150 (MQ=255)                                                                                                                                  agcagTCAGAGACATAACTGGCACGTAGTGGTTGCAACCACTAACCCACCAATAGAGGGGAAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggctatgagcataacgtaatgcttattttaccggacgcatcgccgggaagagaataacgtcgcggatggtatggctgttgggg                                                                                                              >  2:248349‑M1/1‑143 (MQ=255)                                                                                                                                      gTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacag                                                                                                         <  1:4245‑M1/251‑113 (MQ=255)                                                                                                                                                  aCTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacg                                                                                             >  1:249910‑M1/1‑127 (MQ=255)                                                                                                                                                     ggCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggt                                                                                           <  1:270492‑M1/250‑127 (MQ=255)                                                                                                                                                     ggCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacgg                                                                                            <  1:351631‑M1/249‑126 (MQ=255)                                                                                                                                                        aCGTAGGGTTTGCAACCACTAACCCACCCATAGCGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgat                                                                                       >  1:350389‑M1/1‑121 (MQ=255)                                                                                                                                                          gtagggtttgCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgatac                                                                                     >  1:122977‑M1/5‑119 (MQ=255)                                                                                                                                                                       ccACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataccgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgtcggtgcacagcattaccatacggtcgataccaattcccagacc                                                                        >  1:324941‑M1/1‑106 (MQ=255)                                                                                                                                                                            aaCCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgc                                                                     <  1:79623‑M1/249‑149 (MQ=255)                                                                                                                                                                            aaCCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACATGGTTTATGGAAAATTGCCCTCCATCCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattataccgaacgcatcgccgggaacagaataacgtcgcggatggtaaggctgttggtgaacagccttcccatacggtcgataccaattcccccacctgctg                                                                   >  2:290351‑M1/1‑101 (MQ=255)                                                                                                                                                                                   caATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcgg                                                            >  1:357785‑M1/1‑94 (MQ=255)                                                                                                                                                                                     aTAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggt                                                           >  2:366263‑M1/1‑92 (MQ=255)                                                                                                                                                                                        aaggGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGGTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtggataccaattcccagacctgctgtcggcggtaagc                                                       <  2:393553‑M1/250‑163 (MQ=255)                                                                                                                                                                                                                 ccATACTTAACCTGGGTTATGGTAAATTGCCCTCCATTCTGTTTACTTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagactaacgtcgcggatggtatggctgttggtgaacagcattaccatacgatcgataccaaatcccagacctgctgtcggcggcaagccatggtccagtgcggtgccgtaatc                              >  2:149506‑M1/1‑64 (MQ=255)                                                                                                                                                                                                                  cATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatct                             <  2:142174‑M1/251‑189 (MQ=255)                                                                                                                                                                                                                     cctTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtcagccatgttccaatgcggagacggactcttc                           >  2:166046‑M1/2‑60 (MQ=255)                                                                                                                                                                                                                          aCCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaagcttcatcgta                     <  2:101662‑M1/251‑197 (MQ=255)                                                                                                                                                                                                                            cTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgt                      >  1:294640‑M1/1‑53 (MQ=255)                                                                                                                                                                                                                                   aTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaaaatcgc            >  1:256965‑M1/1‑46 (MQ=255)                                                                                                                                                                                                                                         aaTTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcgtc      >  1:211847‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                                          aTTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcggagaacatcgcttcgtcg     >  2:192842‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                                           ttGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttcca                                              <  2:363950‑M1/209‑172 (MQ=255)                                                                                                                                                                                                                                           ttGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttcca                                              >  1:363950‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                                             gCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcgtcgtcc  >  1:68071‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 CCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGT  >  NC_000913/2985570‑2986056

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: ‑