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%	210.4 / NA	61	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 (27/34);  total (27/34)

GCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAAC  >  NC_000913/549554‑549843                                                                                                                                                 |                                                                                                                                                  gCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTgg                                                                                                                                             >  2:69098/1‑151 (MQ=255)         ggCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGcggcgg                                                                                                                                       <  2:351678/149‑1 (MQ=255)         ggCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCggcgggc                                                                                                                                     <  2:265573/151‑1 (MQ=255)          gCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCggcgggcg                                                                                                                                    <  2:659119/151‑1 (MQ=255)           cGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGa                                                                                                                                   <  2:707957/151‑1 (MQ=255)             tGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGCATTGTATGAGCAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGCTAGCGGGGGTGTGGGggggggc                                                                                                                                     >  1:77464/1‑147 (MQ=255)              gCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTa                                                                                                                                <  2:716160/151‑1 (MQ=255)               cTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGGGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGgtgt                                                                                                                                               <  2:370403/135‑1 (MQ=255)                 ttGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTaa                                                                                                                               <  1:382289/149‑1 (MQ=255)                   gAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCt                                                                                                                            <  1:218993/150‑1 (MQ=255)                              ggAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGaaa                                                                                                                 <  2:476808/150‑1 (MQ=255)                                  gCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAgcggc                                                                                                            <  2:803981/151‑1 (MQ=255)                                     cTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCCGATGACGGAAAgcggcg                                                                                                           >  1:435998/1‑149 (MQ=255)                                         cTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACt                                                                                                       >  1:458953/1‑149 (MQ=255)                                           gAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGtt                                                                                                    >  2:696012/1‑150 (MQ=255)                                           gAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGt                                                                                                     <  2:640884/149‑1 (MQ=255)                                           gAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGt                                                                                                     <  1:822142/149‑1 (MQ=255)                                            aaaGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTaa                                                                                                  <  1:445868/151‑1 (MQ=255)                                                 aTGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTg                                                                                             <  2:416417/151‑1 (MQ=255)                                                   gggTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGc                                                                                            >  2:380961/1‑150 (MQ=255)                                                    ggTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCAt                                                                                          <  2:436032/151‑1 (MQ=255)                                                       aTTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAAACTGGTGACGGGAAGCGGCGAAATGGTAACGGTGGATCg                                                                                        >  1:510253/1‑150 (MQ=255)                                                           gTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGGCGGAAAGCGGCGGACTGTTTACGGTGCATCGTCAgg                                                                                   >  1:550762/1‑151 (MQ=255)                                                           gTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAgg                                                                                   >  2:587067/1‑151 (MQ=255)                                                             tGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAggg                                                                                  >  2:115735/1‑150 (MQ=255)                                                             tGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAggg                                                                                  >  2:518540/1‑150 (MQ=255)                                                                    cgcgCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGtt                                                                          <  2:849218/151‑1 (MQ=255)                                                                    cgcgCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGtt                                                                          <  2:340186/151‑1 (MQ=255)                                                                    cgcgCAATGACGATCATCCATCCTCTGCTCGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCGGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGCTAACGTTGCATCGTCAGGGTAGTGGt                                                                           <  1:503945/150‑1 (MQ=255)                                                                         aaTGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGGACTGTTAACCTTGCATCGTCAGGGTAGTGGTTTCgg                                                                      >  1:25596/1‑150 (MQ=255)                                                                         aaTGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAAATGTTAACGTTGCATCGTCAGGGTAGTGGTTTCgg                                                                      >  2:236717/1‑150 (MQ=255)                                                                              cGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCgg                                                                <  2:38190/151‑1 (MQ=255)                                                                                  catccatcCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCggagg                                                             <  1:265054/150‑1 (MQ=255)                                                                                   atccatccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCggagga                                                            >  2:731885/1‑150 (MQ=255)                                                                                   atccatccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGAt                                                           <  2:498357/151‑1 (MQ=255)                                                                                         ccTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGc                                                     <  2:780961/151‑1 (MQ=255)                                                                                           tctGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGGTTAACCTGATGACGGAAAGGGGCGGACTGTTAACGTTGCATCGTCAGGGGAGTGGTTTCGGCCCCGGAGGGTGGGTGCtt                                                   >  2:174186/1‑151 (MQ=255)                                                                                                    cAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTgcg                                          >  2:29732/1‑151 (MQ=255)                                                                                                     agtagtAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGGAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTgcg                                          >  2:325021/1‑150 (MQ=255)                                                                                                      gtagtaGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGGACTGTTAACGTTGCCTTGTCAGGGGAGTGGGTTCGGCCCCGGAGGATGGGGGGTTCGcc                                               >  1:537404/1‑144 (MQ=255)                                                                                                        agtagCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGGAAGCGGCGGACTGTTAAAGTTGCATCGTCAGGGGAGTGGTTTCGGCCCCGGGGGATGGGGGCTTCGCCGGGCGCaa                                       >  1:611381/1‑150 (MQ=255)                                                                                                        agtagCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCaa                                       <  2:114649/150‑1 (MQ=255)                                                                                                        agtagCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAAt                                      <  2:463772/151‑1 (MQ=255)                                                                                                           agCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTc                                    <  2:423676/150‑1 (MQ=255)                                                                                                           agCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCg                                   <  2:287003/151‑1 (MQ=255)                                                                                                                    aaTTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTAt                          <  2:423344/151‑1 (MQ=255)                                                                                                                     aTTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTAt                          >  2:435177/1‑150 (MQ=255)                                                                                                                     aTTATCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTAt                          >  2:136162/1‑150 (MQ=255)                                                                                                                        aTCGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCgg                      <  1:742913/151‑1 (MQ=255)                                                                                                                          cGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGGAAGCGGGGAACTGTTAAAGTTGCATCGTCAGGGTAGGGGTTTCGGCCCCGGAGGATGGGGGCTTTGCCGTGCGCAATTTGATGCGGTATGCgggg                    >  1:460713/1‑149 (MQ=255)                                                                                                                          cGCCAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTg                    <  1:163049/151‑1 (MQ=255)                                                                                                                           gCCAGTCCTGTCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCCGATGACGGAGAGTGGGGGACCGTTAACGTTGCATCGTCCGGGGAGTGGGTTCGGCCCCGGAGGATGGGGGCCTTGCCGTGTGCAATTCGATGCGCTATGCGgggg                   >  1:275938/1‑151 (MQ=255)                                                                                                                             cAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGgggg                   >  2:445305/1‑149 (MQ=255)                                                                                                                             cAGTCCTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGAt                 <  2:63062/151‑1 (MQ=255)                                                                                                                                 ccTGGCGGTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGTAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCTGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTAt              <  1:548987/150‑1 (MQ=255)                                                                                                                                        gTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATggg                                                        >  2:168493/1‑101 (MQ=255)                                                                                                                                        gTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGt                                                       <  1:168470/102‑1 (MQ=255)                                                                                                                                        gTTAGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAAt      <  2:170370/151‑1 (MQ=255)                                                                                                                                           aGCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGa    >  2:856313/1‑150 (MQ=255)                                                                                                                                            gCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAAc  <  2:36286/151‑1 (MQ=255)                                                                                                                                            gCGGGGGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGGGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAAc  >  1:593921/1‑151 (MQ=255)                                                                                                                                              ggggGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAAt      >  2:269031/1‑145 (MQ=255)                                                                                                                                              ggggGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAAt      <  1:269006/145‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGGCGAACTGTTAACGTTGCATCGTCAGGGTAGTGGTTTCGGCCCCGGAGGATGGGTGCTTCGCCGTGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAAC  >  NC_000913/549554‑549843

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