Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,785,316	C→T	R19R (CGC→CGT)	gpmM →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,785,316	0	C	T	100.0%	121.2 / NA	39	R19R (CGC→CGT)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

CGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGAAAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGC  >  NC_000913/3785079‑3785562                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       cGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCACGACGAAAAAAAAGCGCAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAgg                                                                                                                                                                                                                                           >  1:121398/1‑251 (MQ=255)         acatccagcGCGGTAACTGGAATGCAATTTTGCGGGGCACGACGAAAAAAAAGCGCAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGTATAACGc                                                                                                                                                                                                                                    <  2:270617/246‑1 (MQ=255)              cagcGCGTTAACTGGAATGCAATTTTGCGGGGCACGACGAAAAAAAAGCGCAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCAtttt                                                                                                                                                                                                                              <  2:329177/251‑1 (MQ=255)                gcgcGGTAACTTGAATGCAATTGTGCGGGGCACGACGAAAAAAAAGCGCAAATCTGGCATCGGCTTGCCGCGAGCTGTCGTATATTTTGATCTGCAGAGGCTATTTTATCTATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCGTGAGTATTAGGTTGGCGCAATGCCGGGTTCTAAAAAACTTGTGTCACTGGGGTTTCTGTATGGCTATGGCTAGCTTGAAGTACAGCAGTATAGCGCCATTTTTa                                                                                                                                                                                                                            <  2:428853/251‑1 (MQ=255)                                gCAATTTTGCGGGGCACGACGAAAAAAAAGCGCAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCggg                                                                                                                                                                                                            >  1:213765/1‑250 (MQ=255)                                              cacgacgaAAAAAAAGCGCAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACtgt                                                                                                                                                                                              <  1:154843/249‑1 (MQ=255)                                                              cgcaaaTCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGc                                                                                                                                                                               >  1:338789/4‑250 (MQ=255)                                                               gcaaaTCAGGCATCGGCCTGCCGCGAGCGGTCGTATATTTTGATCCTCAGAGGCTATTTTAACCATTAAGCTGTAGTAAAATTACGCCAATTTTGACTCTTGAGTATGAGGTTGTAGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCAt                                                                                                                                                                             >  2:222820/3‑251 (MQ=255)                                                                                  gCCGCGAGCGGTCGTATAGTTTGATCGTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTcc                                                                                                                                                           <  2:10991/250‑1 (MQ=255)                                                                                       gAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGCTTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTg                                                                                                                                                     >  2:306634/1‑251 (MQ=255)                                                                                        aGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCCACTTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGGGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGACATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTgg                                                                                                                                                    >  2:107959/1‑251 (MQ=255)                                                                                         gCGGTCGTATATTTTGAGCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGa                                                                                                                                                   <  2:344064/251‑1 (MQ=255)                                                                                              cGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTc                                                                                                                                               <  1:50681/250‑1 (MQ=255)                                                                                               gaatatTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCgg                                                                                                                                             >  2:109075/3‑251 (MQ=255)                                                                                                 atatTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGAACTGGTGATTCTGGATGGCTATGGCTATCGTGAGGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTACTGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGCCGGTc                                                                                                                                           >  1:156849/1‑251 (MQ=255)                                                                                                              cAGATGCTATTTTAGCGATTCAGCTGTAGTAAAATTACGCAAATGTTGACTCTTGAGTATGATGGTGTCGCAATGTCGGTTGCTAAAAAACCGATGGTTCTGGTGATTCTGGTTGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTGTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGTCCAATCCGCCGCATACCCTGATCGACCCGTCCTGTCTGGAAGTCTGGCTGGCGGACCGTCa                                                                                                                              <  2:411176/251‑1 (MQ=255)                                                                                                              cAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGAAATGGATGCACTGTGGGCCAATCGTCCGCATACACTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTc                                                                                                                               >  1:203893/1‑250 (MQ=255)                                                                                                                   gCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTATAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCCATCGTCCGCATACCCTGATCGACGATTCAGGTCTGGAAGTCGGTCTGCCTGACCGGCAGATgg                                                                                                                         >  1:438265/1‑251 (MQ=255)                                                                                                                     tATTTTATCGCTTCAGCTGTAGTAAAATTACGCAACTTTTGACTCTTGCGTCTGAGGTTGTCGCAATGTCGGTTGCTAAAAAACCTATGGTACTGGGGATTCTGGTTGGCTATGGCTATCGTGAAGCACAGAAGGATCACGCCATTTTTCGAGCTAAAACCCCGGTAATGGATGCACTGTGGCCCAATCGTCCGCATACCCTGATCGACGCTTCCGGACAGGAAGACGGTCTGCCTGTCCGTCAGATGGGt                                                                                                                       >  1:423847/1‑251 (MQ=255)                                                                                                                        tttATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAAc                                                                                                                    >  1:77194/1‑251 (MQ=255)                                                                                                                                                           ttGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGcc                                                                                                       >  1:389734/1‑229 (MQ=255)                                                                                                                                                           ttGACTCTTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCGATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAAGCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCGGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGcc                                                                                                       <  2:389734/229‑1 (MQ=255)                                                                                                                                                                tctTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCg                                                                            >  1:255202/1‑251 (MQ=255)                                                                                                                                                                tctTGAGTATGAGGTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCg                                                                            >  2:337192/1‑251 (MQ=255)                                                                                                                                                                             gTTGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGAcc                                                                >  1:346510/1‑250 (MQ=255)                                                                                                                                                                             gTTGTCGCAATGTCGGTGTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGAcc                                                                <  2:346510/250‑1 (MQ=255)                                                                                                                                                                               tGTCGCAATGTCGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGACGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGGGCCGGCCGCCTCGTGGATCAGGACCTg                                                              >  2:268832/1‑250 (MQ=255)                                                                                                                                                                                                       aCCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGACGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAg                                     <  2:57990/251‑1 (MQ=255)                                                                                                                                                                                                         cTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCt                                                                                                                                                                        <  2:310682/118‑1 (MQ=255)                                                                                                                                                                                                         cTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCt                                                                                                                                                                        >  1:310682/1‑118 (MQ=255)                                                                                                                                                                                                          tagggTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTGCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCAGCGTGTa                                                                        <  2:306673/210‑1 (MQ=255)                                                                                                                                                                                                          tATGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTATGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTa                                                                        >  1:306673/1‑213 (MQ=255)                                                                                                                                                                                                           aTGGTACTGGTGATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATc                                  <  2:213765/250‑1 (MQ=255)                                                                                                                                                                                                                actggttgttcgtgATGGCTATGGCTATCGTTAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATAGTCCGCATGCCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTATGCCACGTTAACCTGGGTGCCTGCCGCATCTTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGctt                            <  2:421521/238‑1 (MQ=255)                                                                                                                                                                                                                      gATTCTGGATGGCTATGGCTATCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCTCTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGgtttcttt                       <  2:355337/250‑1 (MQ=255)                                                                                                                                                                                                                               tggctatggctaTCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATcc               >  1:34558/1‑249 (MQ=255)                                                                                                                                                                                                                                ggctatggctatCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGGCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGGTAACCTGGGTGCCGGCCGCATCGTGTATCAGTACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCGAATccag             <  2:295648/250‑3 (MQ=255)                                                                                                                                                                                                                                    atggctatCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGc          <  1:394352/249‑1 (MQ=255)                                                                                                                                                                                                                                          atCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGgc  >  1:208729/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       CGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGAAAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGC  >  NC_000913/3785079‑3785562

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