Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,608,367	T→C	I484T (ATC→ACC)	hyfF →	hydrogenase 4, membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,608,367	0	T	C	100.0%	273.8 / NA	75	I484T (ATC→ACC)	hyfF	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/40);  total (35/40)

CTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCA  >  NC_000913/2608235‑2608501                                                                                                                                     |                                                                                                                                       cTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGAt                                                                                                                                  >  1:696820/1‑139 (MQ=255) cTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGAt                                                                                                                                  >  1:185993/1‑139 (MQ=255)      cgGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTgg                                                                                                                             >  1:1534835/1‑139 (MQ=255)       gggTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTG‑‑‑ACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:741379/139‑1 (MQ=255)          tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  1:1016559/139‑1 (MQ=255)          tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  1:861302/139‑1 (MQ=255)          tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:696820/139‑1 (MQ=255)          tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGACGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:371382/139‑1 (MQ=255)           gtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcgggc                                                                                                                        <  2:691733/139‑1 (MQ=255)            ttAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcgggcg                                                                                                                       >  2:825848/1‑139 (MQ=255)             tAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      <  2:1534835/139‑1 (MQ=255)              aaTGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCt                                                                                                                     <  1:325211/139‑1 (MQ=255)               aTGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCtt                                                                                                                    <  2:204443/139‑1 (MQ=255)                tGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  2:1214471/139‑1 (MQ=255)                tGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  1:825848/139‑1 (MQ=255)                tGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCCGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  1:745608/139‑1 (MQ=255)                  gCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCa                                                                                                                 <  1:1633143/139‑1 (MQ=255)                   cGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:1232556/130‑1 (MQ=255)                   cGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         >  1:1232556/1‑130 (MQ=255)                    gAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcg                                                                                                                           <  2:1595997/127‑1 (MQ=255)                    gAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcg                                                                                                                           >  1:1595997/1‑127 (MQ=255)                         cgccgcAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          <  1:1242355/139‑1 (MQ=255)                            cgcAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   >  2:528439/1‑127 (MQ=255)                            cgcAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  1:528439/127‑1 (MQ=255)                                             ggTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      >  1:869353/1‑107 (MQ=255)                                             ggTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      <  2:869353/107‑1 (MQ=255)                                                  tctcGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgcct                                                                                 >  1:1345777/1‑139 (MQ=255)                                                      ggCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          <  2:1598065/110‑1 (MQ=255)                                                      ggCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          >  1:1598065/1‑110 (MQ=255)                                                           gTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGt                                                                        >  2:248134/1‑139 (MQ=255)                                                           gTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGt                                                                        >  1:751134/1‑139 (MQ=255)                                                              gACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGc                                                                     >  2:250334/1‑139 (MQ=255)                                                                 caccTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCAcc                                                                  <  2:751134/139‑1 (MQ=255)                                                                    cTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTcc                                                                                                        >  2:1439521/1‑98 (MQ=255)                                                                    cTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTcc                                                                                                        <  1:1439521/98‑1 (MQ=255)                                                                      cGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCa                                                             >  2:212273/1‑139 (MQ=255)                                                                      cGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCa                                                             >  1:286378/1‑139 (MQ=255)                                                                              gTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                                     <  1:212273/139‑1 (MQ=255)                                                                                 aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     <  2:24005/88‑1 (MQ=255)                                                                                 aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     >  1:24005/1‑88 (MQ=255)                                                                                 aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGAt                                                                                       >  1:133846/1‑102 (MQ=255)                                                                                 aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGAt                                                                                       <  2:133846/102‑1 (MQ=255)                                                                                  ttCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  >  2:17562/1‑90 (MQ=255)                                                                                  ttCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  <  1:17562/90‑1 (MQ=255)                                                                                    ctgctgGGCATGATGCTGGCGATGGGGACCCATATTCCACAGCCTGTCACCAGGATCCTGGCGGGCGCTTCCACTAAAGACCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCtt                                               <  2:520218/139‑1 (MQ=255)                                                                                      gctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCa                                             >  2:1436211/1‑139 (MQ=255)                                                                                      gctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCa                                             <  2:214887/139‑1 (MQ=255)                                                                                      gctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCa                                             >  1:214887/1‑139 (MQ=255)                                                                                       ctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACtata                                                                                                            <  2:1336352/75‑1 (MQ=255)                                                                                       ctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACtata                                                                                                            >  1:1336352/1‑75 (MQ=255)                                                                                           tCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCTGGCAc                                        <  2:519832/139‑1 (MQ=255)                                                                                                    tGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     <  2:149362/69‑1 (MQ=255)                                                                                                    tGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     >  1:149362/1‑69 (MQ=255)                                                                                                      gCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      <  1:803314/50‑1 (MQ=255)                                                                                                      gCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      >  2:803314/1‑50 (MQ=255)                                                                                                           gggAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGgag                        >  2:1661290/1‑139 (MQ=255)                                                                                                           gggAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGgag                        >  1:1127923/1‑139 (MQ=255)                                                                                                             gAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  <  1:436192/63‑1 (MQ=255)                                                                                                             gAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  >  2:436192/1‑63 (MQ=255)                                                                                                                     atTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              <  1:1661290/139‑1 (MQ=255)                                                                                                                     atTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              <  2:463415/139‑1 (MQ=255)                                                                                                                              ccTGTCGCCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          <  1:1074193/38‑1 (MQ=255)                                                                                                                              ccTGTCGCCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          >  2:1074193/1‑38 (MQ=255)                                                                                                                               cTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAAcaca                 <  2:65422/126‑1 (MQ=255)                                                                                                                               cTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAAcaca                 >  1:65422/1‑126 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                                     >  2:59104/1‑88 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                                     <  1:59104/88‑1 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCa                                    <  2:1472706/105‑1 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCa                                    >  1:1472706/1‑105 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                             <  2:1243729/112‑1 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                             >  1:1243729/1‑112 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              <  2:654326/127‑1 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              >  1:654326/1‑127 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCa  >  2:907976/1‑139 (MQ=255)                                                                                                                                 gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCa  >  2:907977/1‑139 (MQ=255)                                                                                                                                  tcaccaGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTc                                              >  2:330218/1‑94 (MQ=255)                                                                                                                                  tcaccaGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTc                                              <  1:330218/94‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       CTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCA  >  NC_000913/2608235‑2608501

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