Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,591,647	G→A	R14H (CGC→CAC)	dapE →	succinyl‑diaminopimelate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,591,647	0	G	A	97.4%	118.5 / NA	39	R14H (CGC→CAC)	dapE	succinyl‑diaminopimelate desuccinylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/12);  minor base T (0/1);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCG  >  NC_000913/2591418‑2591892                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      cGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCtg                                                                                                                                                                                                                                   >  2:210477/1‑250 (MQ=255)  gTGGTACAACCTGGAGTAAACTGGTCGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTACTGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGATTCAGTGATTCCAGTTCTCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTACTTCAACGCCCTTCCCTGCGTCCtgat                                                                                                                                                                                                                                 >  2:122325/1‑251 (MQ=255)     gTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAACATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTgatgat                                                                                                                                                                                                                              >  2:73818/1‑251 (MQ=255)                 gTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGcc                                                                                                                                                                                                                    >  1:51277/1‑249 (MQ=255)                          aCGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTttgttg                                                                                                                                                                                                         <  2:142270/251‑1 (MQ=255)                                    ccGCAATAAAAGCTCCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGGCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGAGTCCAGTTATCAGCAAGTTTTCCATGAGGTGTAGTCTATGTCGTTCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGTCCTTCCGTGAGTCCTGTTGATGCAGTATGCCAGGCTTTGTTGATTGAACGtt                                                                                                                                                                                               <  2:377331/251‑1 (MQ=255)                                             aaTCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCg                                                                                                                                                                                       >  2:296240/1‑250 (MQ=255)                                                    gATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGttt                                                                                                                                                                               >  1:284702/1‑251 (MQ=255)                                                         ggCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCACACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCAAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGCGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCg                                                                                                                                                                          >  1:50697/1‑251 (MQ=255)                                                             tctGCGGCGGCATTAATGACTGAAATGCCTGCACTTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGCTGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTg                                                                                                                                                                       >  2:191525/1‑250 (MQ=255)                                                               tGCGGCGGCATTAACGTCTGACATGCCTGCCATTATCAAACGTCCATTCCTCTGAGTGCACGGTAAGCCTATTCTTCTGGGATTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGCCACCACCGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGGCCGGCTTTGTTGGTTGGACCTTTGCAGGCGATCGGTTTTACCGTTGAAc                                                                                                                                                                    >  2:261786/1‑251 (MQ=255)                                                                      gcATTATTGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGAt                                                                                                                                                                                                                       >  2:376984/1‑193 (MQ=255)                                                                      gcATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGAt                                                                                                                                                                                                                       <  1:376984/193‑1 (MQ=255)                                                                        aTTAATGACTGAAATACCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGATTATTGAGCTGACCCAAAAGCTTATTCACCGCCCTTCCCTGCGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCCGGCGTTCGTTTTTACCGTTGAACACATGGACt                                                                                                                                                           >  2:82850/1‑251 (MQ=255)                                                                                          tGCAATTATCAAACGTTCATTGCTCTGAGCGCCCGGTAAGACTATGCTGCTGGGATTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGCGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTCTTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCAAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGAATTTGCGGATAGGCAGAAtt                                                                                                                                         >  2:347778/1‑251 (MQ=255)                                                                                                  tCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGc                                                                                                                                  >  1:138348/1‑250 (MQ=255)                                                                                                   caaccggccATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGAGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCAt                                                                                                                                >  2:316637/8‑251 (MQ=255)                                                                                                         tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTg                                                                                                                          >  1:259176/1‑251 (MQ=255)                                                                                                                        gCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGt                                                                                                           <  1:16112/251‑1 (MQ=255)                                                                                                                           cGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGtt                                                                                                          >  1:74977/1‑249 (MQ=255)                                                                                                                           cGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGtt                                                                                                          >  1:325696/1‑249 (MQ=255)                                                                                                                            ggTAAGCCTATGCTGCTGGGCTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGGTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGc                                                                                                       >  1:90333/1‑251 (MQ=255)                                                                                                                                    tATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAACCGTTAGCATTTGCCg                                                                                               >  2:246223/1‑251 (MQ=255)                                                                                                                                     aTGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGGTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGcc                                                                                                <  2:288193/249‑1 (MQ=255)                                                                                                                                             gggTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATAc                                                                                        <  1:40759/249‑1 (MQ=255)                                                                                                                                                       gATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTgcc                                                                            <  1:210477/251‑1 (MQ=255)                                                                                                                                                             atttATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTgg                                                                      <  2:116443/249‑1 (MQ=255)                                                                                                                                                                     gCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGCTGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCg                                                              >  2:298677/1‑251 (MQ=255)                                                                                                                                                                          ttttCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCg                                                          <  2:11606/250‑1 (MQ=255)                                                                                                                                                                                 tGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATc                                                   >  1:287728/1‑250 (MQ=255)                                                                                                                                                                                         aGTCTATGGCGTGCCCGGTTATTGAGCTGACACAACGGCTTATTCACCGCCCTTCCCTGAGTCCTTATGATGCAGGATGCCAGGCTTTGTTGATTGAACGGTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCAGACCGACGTGGTGCCTACTGGCTACGCCGAGCGTTGGAGCAATccccc                                           <  2:147434/250‑1 (MQ=255)                                                                                                                                                                                          gTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATCGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGt                                         >  2:90484/1‑251 (MQ=255)                                                                                                                                                                                              aTGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGTTTGa                                     <  1:72435/251‑1 (MQ=255)                                                                                                                                                                                                         ggTTATTGAGCTGACACCACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTc                          >  1:95293/1‑251 (MQ=255)                                                                                                                                                                                                               tGAGCTGACACAACAGCTTATTCACCGCCCTTCCCTGAGTCCTGAGGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGGTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTTGCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCCCCATTCGTGAc                     <  1:261786/250‑1 (MQ=255)                                                                                                                                                                                                                       cacaACAGCTTATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAGTTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGtt             >  1:350165/1‑250 (MQ=255)                                                                                                                                                                                                                                 tATTCTCCGCCCTTCCCTGAGGCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGGTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCACGTTATTCGGgcg    <  2:24434/249‑1 (MQ=255)                                                                                                                                                                                                                                 tATTCACCGCCCTTCCCTGAGTCCTGATGATGCCGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGcgcg  >  1:160423/1‑251 (MQ=255)                                                                                                                                                                                                                                 tATTCACCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGcgcg  >  2:299427/1‑251 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGCGTGGGCAGGGTGAAACGTTAGCCTTTGCCGGGCATACCGACGTGGTGCCGCCTGGCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCG  >  NC_000913/2591418‑2591892

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