Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,985,726	A→G	intergenic (‑134/‑121)	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,726	0	A	G	97.6%	127.2 / ‑6.6	42	intergenic (‑134/‑121)	yqeF/yqeG	putative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/19);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCTTTCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGT  >  NC_000913/2985486‑2985959                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          ccTTTCACGACCATACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAAc                                                                                                                                                                                                                                 <  2:403308/250‑1 (MQ=255)         gACCAGACTACCAAGTGCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATATGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCTCGCCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTGTTATCACAAAGGAAATATGCCTGAGAAGCAGTCAGAGACATAACTGGCTCTTAGGGTTTGCAACCACTAAcc                                                                                                                                                                                                                         <  2:118786/251‑1 (MQ=255)                  aCCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACACCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACACAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCCCCAATa                                                                                                                                                                                                                 <  2:286092/250‑1 (MQ=255)                     aaGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAgg                                                                                                                                                                                                             >  1:185839/1‑251 (MQ=255)                            aCGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGAt                                                                                                                                                                                                      <  1:367778/251‑1 (MQ=255)                             cggcggAAGGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACAGAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATa                                                                                                                                                                                                     <  2:72399/251‑1 (MQ=255)                                   aaTGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGt                                                                                                                                                                                               <  1:396561/251‑1 (MQ=255)                                    aTGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGtt                                                                                                                                                                                              <  1:253427/251‑1 (MQ=255)                                                cGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTATTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGACTGAGAAGCAGTCAGAGACATAACTGGCACGAAGGGATGGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGTGTTAATCTTCCATAc                                                                                                                                                                                  >  2:42989/1‑251 (MQ=255)                                                       cGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCt                                                                                                                                                                           >  1:147890/1‑251 (MQ=255)                                                                 cGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCGCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTTGTTTATGGt                                                                                                                                                                 <  2:362927/251‑1 (MQ=255)                                                                         gtCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCGCCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTTGTTTATGGTAAATTGcc                                                                                                                                                         <  1:217658/251‑1 (MQ=255)                                                                                        aCAATCACAACGTCTTTCATCACAACCGCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTaa                                                                                                                                          <  1:344149/251‑1 (MQ=255)                                                                                                         catcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCTTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGTCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTCAATTTGTAGATGATACGtt                                                                                                                         >  1:47212/1‑251 (MQ=255)                                                                                                         catcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTGACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGtt                                                                                                                         >  1:6366/1‑251 (MQ=255)                                                                                                            cacaACCTCGCACCAATTAACAGTGCAATAATAAATGATTGTTAATTACAGTTAGCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATTCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATGGCCCTCCAGTCTGGTTAATTTGTAGATGATACGTTCta                                                                                                                      <  2:155907‑M1/251‑3 (MQ=255)                                                                                                                    cGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAATGTGATTTTTATCTCAAAGGAAATATGCCTGAGCAGCAGTCAGATACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaacc                                                                                                              <  1:417573‑M1/251‑11 (MQ=255)                                                                                                                    cGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaac                                                                                                               <  1:215909‑M1/250‑10 (MQ=255)                                                                                                                           aTTAACAATGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtacaaaaaccccgacat                                                                                                       >  2:207444‑M1/1‑234 (MQ=255)                                                                                                                              aaCAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttg                                                                                                    >  1:381759‑M1/1‑231 (MQ=255)                                                                                                                                             aaTGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTAAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTACTTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagc                                                                                      >  2:263411‑M1/1‑216 (MQ=255)                                                                                                                                                      taataaCAGTTATCTTAAGTGTTTAAAAAAAGTGATTTTTAGCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATCACTGTCACGTAGCGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTGCtaaaaaaaccccgacatttgccggggttgtgagcagaacgtaa                                                                             <  2:144405‑M1/250‑44 (MQ=255)                                                                                                                                                                tATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatttt                                                                   >  2:421078‑M1/1‑197 (MQ=255)                                                                                                                                                                 aTCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaacccccgccatttcccggggttgtgagcaaaaggaaatgcttattttcc                                                                 >  2:40251‑M1/1‑196 (MQ=255)                                                                                                                                                                 aTCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTAAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttac                                                                 >  2:449806‑M1/1‑196 (MQ=255)                                                                                                                                                                                  aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaatcccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggga                                                <  1:4759‑M1/251‑73 (MQ=255)                                                                                                                                                                                  aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggg                                                 >  2:10633‑M1/1‑179 (MQ=255)                                                                                                                                                                                   aaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTATGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGGAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattgtaccggacgcatcgccgggaa                                               <  2:325838‑M1/251‑74 (MQ=255)                                                                                                                                                                                   aaGTGATTGTTATCACAAAGGAAATTGGCCTGAGCAGAAGTCAGAGACATAATGGGCACGTAGGGGTTGCAACCACTAACCCACCAATAGATGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTAGGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaagaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaa                                               <  2:426071‑M1/251‑74 (MQ=255)                                                                                                                                                                                             gctcaCAAAGGCAATATGCCTGAGGAGCAGTCAGAGAGATAACTGGCACGTAGGGGTTGCAACCACTAACCCACCAATAAAGGGGGAGATAGAGCGTTAATCTCCCATACTTAACATGAGTTATGGGAAAGTGCACTCCATTCTGTTTAATTTGTAGCTGATACGTTCtaaaaaaaccccgacatttcccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacg                                     >  2:227‑M1/3‑168 (MQ=255)                                                                                                                                                                                              aTCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgt                                    >  1:336754‑M1/1‑167 (MQ=255)                                                                                                                                                                                                  caAAGGAAATATGCCTGAGCAGCAGTCAGAGATATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGACCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgacgcg                                >  2:236985‑M1/1‑163 (MQ=255)                                                                                                                                                                                                             tGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAAGTGCCCTACATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgcggttgtgagcataacgtaatgcttattttaccggacgcatcgcctggcacagaataacgtcgcggatggtatggc                     >  2:190541‑M1/1‑152 (MQ=255)                                                                                                                                                                                                             tGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGATTGCCACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggc                     >  2:131314‑M1/1‑152 (MQ=255)                                                                                                                                                                                                               ccTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTAAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGACACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctg                   >  1:333966‑M1/1‑150 (MQ=255)                                                                                                                                                                                                                  gAGCAGCAGTCAGAGACATAACTGGCCCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatttgaccggacgcatcgccgggaacagaataacgtcgctgatggtatggctgtt                 <  2:124152‑M1/250‑104 (MQ=255)                                                                                                                                                                                                                      agcagTCAGAGACATAACTGGCACGTAGGGATTGCAACCACTAATCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaacaaccccgacatttgcctgggttgtgagcataacgtaatgcttattttaccgcacgcatcgccgggaacagataaacgtagaggatggtatggctgtttgtgg            >  2:227565‑M1/1‑143 (MQ=255)                                                                                                                                                                                                                           tCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGc                                                                                                                                                          >  2:220944/1‑104 (MQ=255)                                                                                                                                                                                                                           tCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGc                                                                                                                                                          <  1:220944/104‑1 (MQ=255)                                                                                                                                                                                                                                gaCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgccgcggatggtatggctgttggtgaacagcattac  >  1:99865‑M1/1‑133 (MQ=255)                                                                                                                                                                                                                                            aCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcg                                  <  2:310550‑M1/207‑87 (MQ=255)                                                                                                                                                                                                                                            aCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcg                                  >  1:310550‑M1/1‑121 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          GCTTTCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGT  >  NC_000913/2985486‑2985959

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