Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	549,698	A→G	G24G (GGA→GGG)	ylbF →	putative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	549,698	0	A	G	100.0%	217.2 / NA	63	G24G (GGA→GGG)	ylbF	putative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/37);  total (26/37)

GGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGA  >  NC_000913/549553‑549841                                                                                                                                                  |                                                                                                                                                ggCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGtgtg                                                                                                                                            <  2:830046/151‑1 (MQ=255)   cACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTgg                                                                                                                                           >  1:499784/1‑150 (MQ=255)       gCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGGGTGgcggcg                                                                                                                                      >  2:260228/1‑151 (MQ=255)             gTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGAt                                                                                                                                >  2:834409/1‑151 (MQ=255)                  ttGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGGTTAAcc                                                                                                                           >  2:321672/1‑151 (MQ=255)                  ttGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGGGTGGCGGCGGGCGATTaaac                                                                                                                           >  2:419563/1‑149 (MQ=255)                    gAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCtg                                                                                                                         <  1:868245/151‑1 (MQ=255)                       cAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCtgat                                                                                                                       >  1:594187/1‑150 (MQ=255)                        aGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGGTTAACCtgatg                                                                                                                      >  1:223044/1‑150 (MQ=255)                        aGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCtgatga                                                                                                                     <  2:817579/151‑1 (MQ=255)                        aGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCtgatg                                                                                                                      >  1:619061/1‑150 (MQ=255)                          gcgcTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACg                                                                                                                   <  1:260254/151‑1 (MQ=255)                            gcTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGa                                                                                                                 <  1:888628/151‑1 (MQ=255)                              tGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACggga                                                                                                                >  2:336155/1‑148 (MQ=255)                                 aaGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAgcg                                                                                                            <  2:145884/151‑1 (MQ=255)                                 aaGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAgcg                                                                                                            <  2:524443/151‑1 (MQ=255)                                    cACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAgcggcg                                                                                                         <  2:346341/151‑1 (MQ=255)                                    cACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAgcgg                                                                                                           >  1:517660/1‑149 (MQ=255)                                        ggCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAAc                                                                                                      <  1:452701/150‑1 (MQ=255)                                         gCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTg                                                                                                    <  2:109615/151‑1 (MQ=255)                                             aaaGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTaa                                                                                                <  1:336181/151‑1 (MQ=255)                                               aGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGGAAGCGGCGAACTGTTAACg                                                                                              >  2:254445/1‑151 (MQ=255)                                               aGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACg                                                                                              <  2:760943/151‑1 (MQ=255)                                                gCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGt                                                                                             <  1:396990/151‑1 (MQ=255)                                                 cATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGGGGGGGTGTGGCGGCGGGCGATTAACCTGGCGAAGGGAAGCGGGCCACTGTTCAcct                                                                                             >  1:180128/1‑148 (MQ=255)                                                        aTTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCg                                                                                      <  2:533198/150‑1 (MQ=255)                                                         ttGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGt                                                                                     <  2:592425/150‑1 (MQ=255)                                                           ggTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAg                                                                                  <  2:516209/151‑1 (MQ=255)                                                                     cgcgCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGttgagg                                                                          >  1:844311/1‑144 (MQ=255)                                                                       cgcAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTc                                                                      <  1:164642/151‑1 (MQ=255)                                                                                 atcatcCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCggag                                                            <  1:11983/151‑1 (MQ=255)                                                                                    atccatccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGGAAGCGGCGGACTGTTAACGGTGCATCGTCAGGGTAGTGGTTTCGGCCCCggagg                                                           >  1:235849/1‑149 (MQ=255)                                                                                    atccatccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCggagg                                                           >  1:93958/1‑149 (MQ=255)                                                                                      ccatccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATgg                                                       <  2:834446/151‑1 (MQ=255)                                                                                          ccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGGACTGTTAACGTTGCATCGTCAGGGTAGGGGTTTCGGCCCCGGAGggtgggtg                                                    >  1:316813/1‑150 (MQ=255)                                                                                          ccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTg                                                    <  1:529915/150‑1 (MQ=255)                                                                                              tGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTc                                                >  1:903937/1‑150 (MQ=255)                                                                                                cTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGGTGACGGAAAGCGGCGAACTGTTATCGGTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGGGCtt                                                 >  2:178544/1‑147 (MQ=255)                                                                                                  tGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGcc                                             <  1:664115/149‑1 (MQ=255)                                                                                                  tGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCg                                            <  1:248024/150‑1 (MQ=255)                                                                                                     cAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTgcg                                        <  1:788214/151‑1 (MQ=255)                                                                                                      agtagtAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTgcgc                                       <  2:755231/151‑1 (MQ=255)                                                                                                            agCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTc                                  <  1:389882/150‑1 (MQ=255)                                                                                                            agCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTc                                  <  2:551572/150‑1 (MQ=255)                                                                                                              cgcACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGAt                               <  1:439847/151‑1 (MQ=255)                                                                                                                cACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATg                              <  1:404872/150‑1 (MQ=255)                                                                                                                  ccGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGGTGACGGAAAGCGGCGAACTGTTAACGTTGCATCGGCCGGGTAGGGGTTTCGGCCCCGGCGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGt                           >  1:78553/1‑151 (MQ=255)                                                                                                                  ccGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCg                            <  1:800831/150‑1 (MQ=255)                                                                                                                  ccGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGt                           <  2:85731/151‑1 (MQ=255)                                                                                                                    gAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGtt                          >  1:393280/1‑150 (MQ=255)                                                                                                                        tATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATg                       >  2:297914/1‑149 (MQ=255)                                                                                                                         aTCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCgg                    <  1:654099/151‑1 (MQ=255)                                                                                                                           cGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTg                  <  2:754210/151‑1 (MQ=255)                                                                                                                           cGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTg                  <  1:111780/151‑1 (MQ=255)                                                                                                                            gCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTgg                 <  1:517547/151‑1 (MQ=255)                                                                                                                                gTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGGTTACGGGAAGCGGCGAACTGGTAACGGTGCATCGTCCGGGTAGGGGGTTCGGCCCCGGGGGGTGGGGGCCTCGCCGTGCGCCATTCGATGGGTTATGCGGGGGGttt             >  2:872316/1‑150 (MQ=255)                                                                                                                                  ccTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTAt            <  1:321698/150‑1 (MQ=255)                                                                                                                                   cTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTTACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGGTGGGTGCTTCGCCGTGCGCAATTCCATGTGTTATGCGGTGGATTAt            >  2:574528/1‑149 (MQ=255)                                                                                                                                       cGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGc       >  2:769763/1‑150 (MQ=255)                                                                                                                                       cGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGc       >  1:619856/1‑150 (MQ=255)                                                                                                                                        ggTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCGTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCa      <  2:609259/150‑1 (MQ=255)                                                                                                                                            aGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGa  <  1:523357/150‑1 (MQ=255)                                                                                                                                            aGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGa  >  1:600321/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                GGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGA  >  NC_000913/549553‑549841

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