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	100.0%	159.9 / NA	44	intergenic (‑224/‑31)	yqeF/yqeG	putative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/18);  total (26/18)

ACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATC  >  NC_000913/2985573‑2986049                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          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)                                                                                                                                         gaCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgccgcggatggtatggctgttggtgaacagcattac                                                                                            >  1:99865‑M1/1‑133 (MQ=255)                                                                                                                                                     aCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcg                                                                                                                            >  1:310550‑M1/1‑121 (MQ=255)                                                                                                                                                     aCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcg                                                                                                                            <  2:310550‑M1/207‑87 (MQ=255)                                                                                                                                                        tagggtttGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgatacc                                                                             >  1:369818‑M1/4‑118 (MQ=255)                                                                                                                                                         agggttTGCAACCACTAACCCACCAATAGAGGGGCAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTACATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgtttgggaacagcattaccatacggtcgaaaacc                                                                            >  2:319520‑M1/3‑117 (MQ=255)                                                                                                                                                                gCAACCTCTAACCCACCAAGAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGTTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtactgcttattttaccgggcgcatcgccgggaatagaataacgtctcggatggtatggctgttggtgaacatcattaccatacggtcgataccaattccc                                                                      <  2:425793‑M1/250‑141 (MQ=255)                                                                                                                                                                     cACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacct                                                                <  1:310627‑M1/251‑147 (MQ=255)                                                                                                                                                                              accaATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggc                                                       <  2:62557‑M1/251‑156 (MQ=255)                                                                                                                                                                                     gagGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgatac                                                                              <  1:367105‑M1/221‑133 (MQ=255)                                                                                                                                                                                     gagGGGTAGATAGGGAGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgatac                                                                              >  2:367105‑M1/1‑89 (MQ=255)                                                                                                                                                                                              aTAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgtttgtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaaaccagtttcc                                        >  1:418716‑M1/1‑80 (MQ=255)                                                                                                                                                                                                  ggCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacagttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattccaagacctgctgtcggcggtaagccatgttccagtgc                                   <  1:298775‑M1/251‑176 (MQ=255)                                                                                                                                                                                                        aaTCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgac                             <  1:449806‑M1/251‑182 (MQ=255)                                                                                                                                                                                                                    ttAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgaagtaatcttcatc                 >  1:128456‑M1/1‑58 (MQ=255)                                                                                                                                                                                                                      aaCCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcattgg               >  2:25587‑M1/1‑56 (MQ=255)                                                                                                                                                                                                                         cTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgaaatcttcatcgt               >  2:173282‑M1/1‑53 (MQ=255)                                                                                                                                                                                                                                aTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcg      <  1:326921‑M1/250‑205 (MQ=255)                                                                                                                                                                                                                                   gTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgtaccagtgcggtgacgtaatcttcatcgtagaacatcgcttc  >  1:127064‑M1/1‑43 (MQ=255)                                                                                                                                                                                                                                   gTAAAATGCCCTCCAATCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttc  <  1:40251‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          ACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATC  >  NC_000913/2985573‑2986049

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