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.8%	300.8 / NA	81	intergenic (‑224/‑31)	yqeF/yqeG	putative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (45/35);  minor base A (1/0);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAG  >  NC_000913/2985571‑2986054                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               cGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGGGCAATAGTAAATGATTGTTAATTACAGTTAGCTTAAGTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGATCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATATAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAAGTGCCCTCCATTCTGttt                                                                                                                                                                                                                                           <  2:78184/251‑1 (MQ=255)  gACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTa                                                                                                                                                                                                                                          <  1:84993/251‑1 (MQ=255)       tCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTg                                                                                                                                                                                                                                     >  2:95758/1‑251 (MQ=255)              gTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAgatga                                                                                                                                                                                                                              <  1:27135/251‑1 (MQ=255)                  ttCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAgatgat                                                                                                                                                                                                                             <  1:36787/248‑1 (MQ=255)                  ttCATCAAAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAACAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATATAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACg                                                                                                                                                                                                                          >  2:130925/1‑251 (MQ=255)                   tcatcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACg                                                                                                                                                                                                                          <  1:223629/250‑1 (MQ=255)                    catcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGtt                                                                                                                                                                                                                        <  2:57985/251‑1 (MQ=255)                     ataacaACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATATCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTc                                                                                                                                                                                                                       >  2:371729/4‑251 (MQ=255)                       cacaACCTCGCACCAATTAACAGTGAAATAGTAAATGATTGTTAAATACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTAGCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCGCCAATAGAGGGGTAGATATGGCGTTAATCTTCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGAGACGTTCt                                                                                                                                                                                                                      <  2:126942‑M1/250‑2 (MQ=255)                        acaACCTCGAACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTCATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGCCACGTAGGGTTTGCAACCAATAACCCACCAATAGAGGGGAAGATAGGGCGTTAATCTCCCATACTTAACCTAGTTTATGGTAAATTGCCCTCAATTCTGTTTAATTTATAGATGATACGTTCta                                                                                                                                                                                                                     >  1:477069‑M1/1‑248 (MQ=255)                          aaCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaa                                                                                                                                                                                                                   >  2:99457‑M1/1‑246 (MQ=255)                                 cACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGAAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTCATTTGTAGATGATACGTTCtaaaaaaacccc                                                                                                                                                                                                           >  2:12344‑M1/1‑239 (MQ=255)                                    cAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATCTGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGTTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCAGTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                                                                                                        <  1:156443‑M1/251‑16 (MQ=255)                                    cAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATATAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTTGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                                                                                                        <  1:96294‑M1/251‑16 (MQ=255)                                         aaCAGTGCGATAGTAAATGATTGGTCATTCCAGTTATCTTAATTGTTTAAAAAAAGGGATTTTTATCACAAAGCAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATCGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTGTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttt                                                                                                                                                                                                   >  2:34796‑M1/1‑231 (MQ=255)                                          aCAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGAGATTTTTATCACAAAAGAAATATGCCTGAGCAGCAGTTAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAACTTGCCCTACATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttac                                                                                                                                                                                                  >  2:34463‑M1/1‑230 (MQ=255)                                            aGTGCAATAGTAAATGATTGTTAATTACAGTAATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAACGAAATATTCCTGAGCAGCAGTCAGAGACATAACTTGCACGTAGGATTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTACCATACTTAACCTGGGTTATGGTAAATTGCCATCCCTTCTGTTTAATTTGTAGATGATACGTTAtaaaaaaaccccgacatttgccg                                                                                                                                                                                                >  2:166427‑M1/1‑228 (MQ=255)                                             gTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATGCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgg                                                                                                                                                                                               <  2:364487‑M1/251‑25 (MQ=255)                                                         aTGATTGTTAATTACAGTTCTCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcat                                                                                                                                                                                   >  1:364389‑M1/1‑215 (MQ=255)                                                         aTGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcat                                                                                                                                                                                   >  1:486294‑M1/1‑215 (MQ=255)                                                           gATTGTTAATTACAGTTATCTTAAGTGGTTAAAAAAAATGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaatccccgacatttgccggggttttgagcatat                                                                                                                                                                                 <  2:132715‑M1/251‑39 (MQ=255)                                                             ttGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacg                                                                                                                                                                               >  2:40564‑M1/1‑211 (MQ=255)                                                               gTTAATTACATTTATCTTAATTGATTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGAAGGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGAGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaacccgacatttgccggggttgtgagcataacgt                                                                                                                                                                              >  2:469474‑M1/1‑209 (MQ=255)                                                               gTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgta                                                                                                                                                                             <  2:289982‑M1/251‑43 (MQ=255)                                                                  aaTTACAGTTATCTGAATTGGTTGAAAAAATTGATGTTTATCACAAAGGAAAGATTCCTGAGCAGCAGTCAGTGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGGTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCAGTCTGTTTAATTTGTTGATGATACTGTCgaaaaaaaccccgacatttgccggggttgtgatcataacgtaaa                                                                                                                                                                           <  2:27701‑M1/250‑45 (MQ=255)                                                                   aTTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgc                                                                                                                                                                         >  2:39735‑M1/1‑205 (MQ=255)                                                                   aTTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgc                                                                                                                                                                         <  1:61745‑M1/251‑47 (MQ=255)                                                                                tAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttac                                                                                                                                                                <  1:371729‑M1/247‑56 (MQ=255)                                                                                  aTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacg                                                                                                                                                          >  1:9541‑M1/1‑190 (MQ=255)                                                                                   ttGTTTAAAAAATGTGATTTTTATCACAAAGGCAATATGCCTGAGCAGAAGTCAGAGCCATAACTGGCACGTAGGGTTTGCAACCACAAACCCACCAATAGCGGGGTAGATAGGGCGTTAATCTCTCATACTTAACCTGGTTTATGGTACATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgcaggggttgtgatcatcacgtaatgcttattttaccggacgc                                                                                                                                                         >  2:162452‑M1/1‑189 (MQ=255)                                                                                    tGTTTAAAAAAAGTTATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACTTAGGGTTTTCAACCACTACCCCACCAATAGAGGGGTAGATAGGGCGGTAATCTCCCATACTTAACCTGGTTTATGTTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggtttggagcataacgtaatgcttattttaccggacgc                                                                                                                                                         <  2:412365‑M1/250‑63 (MQ=255)                                                                                          aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccg                                                                                                                                                  >  1:443122‑M1/1‑182 (MQ=255)                                                                                          aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccg                                                                                                                                                  >  1:216703‑M1/1‑182 (MQ=255)                                                                                               aGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaac                                                                                                                                             >  1:133227‑M1/1‑177 (MQ=255)                                                                                                 tGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGAAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagaataacgtaatgcttattttacctgacgcatcgccgggaac                                                                                                                                             >  2:152686‑M1/1‑175 (MQ=255)                                                                                                  gATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacaga                                                                                                                                          <  2:175264‑M1/251‑78 (MQ=255)                                                                                                   aTTTTTATCCCAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgggggtgtgagcctaacagaatgcttatttgaccggacgcatcgcctggaacaga                                                                                                                                          >  2:293‑M1/1‑173 (MQ=255)                                                                                                   aTTTTTATCACAAAGGAACTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaaccga                                                                                                                                          >  2:414934‑M1/1‑173 (MQ=255)                                                                                                     ttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACACACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGATTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaacaaaccccgacatttgccggggttgtgagcataaagtaatgcttattttaccggacgcatcgccgggaacagaat                                                                                                                                        >  1:196554‑M1/1‑171 (MQ=255)                                                                                                        tATCACAAAGTAAATATGCCTGAGCAGCAGTAAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTGTATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggcgttgtgagcataacgtaatgcttattttaccggacgcatcgccggtaacagaataac                                                                                                                                     <  2:471725‑M1/250‑83 (MQ=255)                                                                                                         aTCACAAAGGAAATAGGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacg                                                                                                                                    <  1:218368‑M1/250‑84 (MQ=255)                                                                                                          tCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtc                                                                                                                                  >  1:410968‑M1/1‑166 (MQ=255)                                                                                                           cacaAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtc                                                                                                                                  <  1:118319‑M1/250‑86 (MQ=255)                                                                                                             caAAGTAAATATGCCTGAGCAGCAGGCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGTGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTTTAGATGATACGTTCtaaaaaaaccccgacatgtgccggggttgtgagcataacggaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcg                                                                                                                               <  2:369410‑M1/251‑89 (MQ=255)                                                                                                               aaGGAAATATGACTGAGCAGCAGTCAGAGCCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGAAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGACCTCCATTCTGTTTAATTTGTAGAGGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgagg                                                                                                                              >  2:118254‑M1/1‑161 (MQ=255)                                                                                                                   aaaTATGCCTGAGCAGCAGTCAGAGACATAANTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                                                                                                                         <  2:242286‑M1/251‑95 (MQ=255)                                                                                                                   aaaTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                                                                                                                         >  1:31707‑M1/1‑157 (MQ=255)                                                                                                                   aaaTATGCCTGAGCAGCAGTCAAAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacctttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgg                                                                                                                          >  1:331104‑M1/1‑157 (MQ=255)                                                                                                                   aaaTATGACTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                                                                                                                         >  2:131873‑M1/1‑157 (MQ=255)                                                                                                                     atatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtat                                                                                                                       <  1:88996‑M1/251‑97 (MQ=255)                                                                                                                             gAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtgtggctgttg                                                                                                               <  1:34796‑M1/251‑105 (MQ=255)                                                                                                                              agcagcAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggaatggctgttgg                                                                                                              >  1:74543‑M1/1‑146 (MQ=255)                                                                                                                                     gTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaaca                                                                                                        <  1:400319‑M1/250‑112 (MQ=255)                                                                                                                                     gTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggggaaca                                                                                                        >  2:400319‑M1/1‑139 (MQ=255)                                                                                                                                      tCAGAGACATAACTGGCACGTAGGGATTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagc                                                                                                      >  1:253609‑M1/1‑138 (MQ=255)                                                                                                                                        agagACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgtgatggctgatggtgaaccgca                                                                                                     >  2:458817‑M1/1‑136 (MQ=255)                                                                                                                                          agaCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcatt                                                                                                   >  2:19828‑M1/1‑134 (MQ=255)                                                                                                                                              aTAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccat                                                                                              >  1:161931‑M1/1‑130 (MQ=255)                                                                                                                                              aTAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccat                                                                                              >  1:181142‑M1/1‑130 (MQ=255)                                                                                                                                                    ggCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtc                                                                                        <  1:175762‑M1/251‑128 (MQ=255)                                                                                                                                                          tagggcttgcaaccACTAACCCACCAATAGAGGGGTAGATA‑GGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtggaacaccttaccaaacgggcgataacc                                                                                 >  2:94128‑M1/7‑117 (MQ=255)                                                                                                                                                            gggtTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaa                                                                                <  1:360265‑M1/250‑136 (MQ=255)                                                                                                                                                              gTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaatt                                                                              <  2:159179‑M1/251‑138 (MQ=255)                                                                                                                                                                              ccaccaATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcg                                                              <  1:414934‑M1/251‑154 (MQ=255)                                                                                                                                                                                         ggggTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagcc                                                    >  1:19572‑M1/1‑87 (MQ=255)                                                                                                                                                                                           ggTAGATAGGGCGTTAATCTCCCATTCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatg                                                 <  2:274612‑M1/251‑167 (MQ=255)                                                                                                                                                                                               gATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataaagtcgcggatggtatggctgttggtgaacagcattaccataaggtcgataccaattcacagacctgctgtcggcgggaagccatgttcc                                             >  1:232022‑M1/1‑81 (MQ=255)                                                                                                                                                                                               gATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCATCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatgtgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcgtcggtaagccatgttcc                                             <  2:19572‑M1/251‑171 (MQ=255)                                                                                                                                                                                                  aGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagt                                          >  1:207305‑M1/1‑78 (MQ=255)                                                                                                                                                                                                   gtgccttcatctcccatACTTAACCTGGTTTATTGTAAATTGCCCTCCAGTCTGTTTAATTTGTAGATGATACGTTCtaaacacaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcaccgccgggaacagaataacgtcgcgtatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcgtcggtaagccatgttccattg                                         <  2:354174‑M1/243‑175 (MQ=255)                                                                                                                                                                                                       gTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctnctgtcggcgggaagccattttcccgtgcggt                                     >  2:246662‑M1/1‑73 (MQ=255)                                                                                                                                                                                                                ccATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattccca                                                                          <  1:36596‑M1/205‑142 (MQ=255)                                                                                                                                                                                                                    aCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTAtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggcacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagaactgctgtcggcggtaagccatgtttcagggaggtgacgtaatcttct                        >  2:465857‑M1/1‑60 (MQ=255)                                                                                                                                                                                                                             ggTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaac               <  1:461492‑M1/251‑201 (MQ=255)                                                                                                                                                                                                                                 tATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggggaacagaattaccatacgggcgataccaattccccgacctgctgtcggcggtaagcaatgttccagtgcggggacgtaatcttcatcgcagaacaacg           >  2:269331‑M1/1‑47 (MQ=255)                                                                                                                                                                                                                                  aTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgtttcaatgaggtgacgtaaatttcatcgtagaacaacgc          >  2:124212‑M1/1‑46 (MQ=255)                                                                                                                                                                                                                                  aTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacagtcgataccaattcccagccctgctgtcggcggtaatccatgttccagtgcggtgacgtaaccttcatcgtagaacatcgc          >  2:47768‑M1/1‑46 (MQ=255)                                                                                                                                                                                                                                     gTAAATTGCCCTCCATTCTGTTTATTTTGTANATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatttgtccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatactgtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcggagaacatcgcttc       <  2:229769‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                                      tAAATTGCCCTCCATTATGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttcccgtgcggnggcgtaatcttcatcgtagaaccccgcgttc      >  2:239871‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                                          ttGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccctgttccagtgcggtgacgtaatcttcatcgtagaacaacccttcgtctc  >  2:114839‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               CGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAG  >  NC_000913/2985571‑2986054

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