Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	587,830	C→T	G51D (GGC→GAC)	ybcH ←	putative exopolysaccharide secretion system auxiliary component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	587,830	0	C	T	100.0%	132.4 / NA	39	G51D (GGC→GAC)	ybcH	putative exopolysaccharide secretion system auxiliary component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/16);  total (23/16)

CCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCG  >  NC_000913/587584‑588066                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             ccAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTg                                                                                                                                                                                                                                          >  2:38342/1‑251 (MQ=255)   aGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGcc                                                                                                                                                                                                                                        >  2:322698/1‑250 (MQ=255)           cgcTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAAc                                                                                                                                                                                                                                >  1:18113/1‑251 (MQ=255)               ccATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGa                                                                                                                                                                                                                             <  2:77431/250‑1 (MQ=255)                cATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACt                                                                                                                                                                                                                           >  1:315873/1‑251 (MQ=255)                             gCTGGAGATCGGCAGCCAGCAGGAGATTAAGATAGCTTTCCAGCTCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTg                                                                                                                                                                                                               <  2:381880/250‑1 (MQ=255)                                                  ggCGATTAATATAGCTTGCCAGCTCTGCGAATGACGCTGTCTGGTGCATAAAAAATGCCGGACCGGCGTTCAGCCAGACAATAAGCTTCAGACCAGCCTGTTGCGCTGCTGCGGCCCGCTTAAACAATAACGTGCGATGGTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGGATCGAAGTCTTGCAAACTTAACTTACTCATCAACCCCTGCCACTGTGGATCGGTAATCTGa                                                                                                                                                                                         <  2:405937/251‑1 (MQ=255)                                                      cttAAAATAGCTTTCCAGCGCTGCGGCAGACTGTTTCTGGTGCATAAAAAATTCAGGACCGGCGTTCAGCCCGCCAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGGGCGCTGTTCTGGCTGGGTAAATGCATCGCCGGCACGGGTCCATTGCAAAACAAGGGTATCGAGGTCTTGCAAACGTAACTGACACATCAGCCACTGCCACTGGGAATCGGAAATCTGAcgat                                                                                                                                                                                     >  2:51600/2‑248 (MQ=255)                                                         aaGATAGCTTTCCAGCGCTGCGGACGACTGACTCTGGCGCATAAAAAATTCCGGACCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGCGAGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTATTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTAtctc                                                                                                                                                                                  >  1:175232/1‑251 (MQ=255)                                                           gATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACACGGGAATCGAAGTCTTGCAAACGAAACTGACTCATCAGCCCCTGCCACTGGGAATCGGTAATCTGACTATCTCgg                                                                                                                                                                                >  2:315324/1‑251 (MQ=255)                                                            aTAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGt                                                                                                                                                                               <  1:344909/251‑1 (MQ=255)                                                               gCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAACAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCACAACAAGGGTATCGACGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATAGGTAATCTGACTATCTCGGtttt                                                                                                                                                                            >  2:249167/1‑251 (MQ=255)                                                                            gCGGACGCCTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGcc                                                                                                                                                                   >  2:276578/1‑247 (MQ=255)                                                                                             gTGCATAAAAAATTCCGGATCGGCGTTCAGCCAGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCa                                                                                                                                               >  1:68460/1‑250 (MQ=255)                                                                                              tGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCaa                                                                                                                                              >  1:92208/1‑250 (MQ=255)                                                                                                cATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAgg                                                                                                                                           <  1:180892/251‑1 (MQ=255)                                                                                                                             cGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGc                                                                                                              >  1:284785/1‑251 (MQ=255)                                                                                                                              gACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTCGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGc                                                                                                              >  2:287415/1‑250 (MQ=255)                                                                                                                                        ttCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                                                            >  2:22355/1‑210 (MQ=255)                                                                                                                                        ttCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                                                            <  1:22355/210‑1 (MQ=255)                                                                                                                                              ccAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACtt                                                                                             <  1:144906/251‑1 (MQ=255)                                                                                                                                               cAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACtt                                                                                             >  2:65079/1‑250 (MQ=255)                                                                                                                                                     gTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAt                                                                                       <  1:31174/250‑1 (MQ=255)                                                                                                                                                                   gCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTcc                                                                        <  1:234428/251‑1 (MQ=255)                                                                                                                                                                             aCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTgaga                                                              <  1:243762/251‑1 (MQ=255)                                                                                                                                                                              cAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAAt                                                                                                   <  1:48793/213‑1 (MQ=255)                                                                                                                                                                              cAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAAt                                                                                                   >  2:48793/1‑213 (MQ=255)                                                                                                                                                                                  aaCGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcg                                                          <  1:152478/250‑1 (MQ=255)                                                                                                                                                                                        cgcTGTTCTGGCTGGGTAAATGCAGCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttgtt                                                   <  1:410831/251‑1 (MQ=255)                                                                                                                                                                                                      ggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGc                                                                                                              <  2:272938/178‑1 (MQ=255)                                                                                                                                                                                                      ggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGCATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGc                                                                                                              >  1:272938/1‑178 (MQ=255)                                                                                                                                                                                                          aaTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGACAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTg                                 >  2:277132/1‑251 (MQ=255)                                                                                                                                                                                                              cATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGAtt                              >  1:180201/1‑250 (MQ=255)                                                                                                                                                                                                                        gcGGGTCCATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAGCTTACTCATCAGCCCCTGCCACTGGGTATCGGTAGTCTGACTGTCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCTTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTaa                   <  2:175232/250‑1 (MQ=255)                                                                                                                                                                                                                              ccATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTAt             >  2:63655/1‑251 (MQ=255)                                                                                                                                                                                                                               cATTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATg            >  1:386592/1‑251 (MQ=255)                                                                                                                                                                                                                                aTTGCAAAACAAGGGTATCGAAGTCTTGCAAACGTAACTGCCTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGt           >  2:45170/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaCAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACt    >  2:306110/1‑251 (MQ=255)                                                                                                                                                                                                                                         cAAGGGTATCGAAGTCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAATCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCg  <  2:6845/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             CCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTCAGACCAGCCTGTTGCGCAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCG  >  NC_000913/587584‑588066

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