Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,541,559	G→A	*201* (TAG→TAA)	fimB →	tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,541,559	0	G	A	100.0%	224.9 / NA	69	*201* (TAG→TAA)	fimB	tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/35);  total (34/35)

GCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCA  >  NC_000913/4541430‑4541693                                                                                                                                  |                                                                                                                                       gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                               >  1:79811/1‑139 (MQ=255) gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                               >  1:749555/1‑139 (MQ=255) gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                               <  2:560255/139‑1 (MQ=255)  cgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                              <  2:262898/139‑1 (MQ=255)        aTCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATAcc                                                                                                                        >  2:1007942/1‑139 (MQ=255)               aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                             >  2:93801/1‑127 (MQ=255)               aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                             <  1:93801/127‑1 (MQ=255)               aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                                 <  2:79811/139‑1 (MQ=255)               aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                                 <  1:1007942/139‑1 (MQ=255)                  aTCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGgtg                                                                                                              <  2:68554/139‑1 (MQ=255)                        ggCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaa                                                                                                        >  1:383706/1‑139 (MQ=255)                             cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                                    <  1:782576/106‑1 (MQ=255)                             cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                                    >  2:782576/1‑106 (MQ=255)                             cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaa                                                                                                         >  2:673489/1‑133 (MQ=255)                             cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaa                                                                                                         <  1:673489/133‑1 (MQ=255)                               cAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACat                                                                                                 <  2:606325/139‑1 (MQ=255)                                aaTATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                                <  2:749555/139‑1 (MQ=255)                                aaTATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                                <  1:540565/139‑1 (MQ=255)                                       gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGAt                                                                                                                  <  2:261262/114‑1 (MQ=255)                                       gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGAt                                                                                                                  >  1:261262/1‑114 (MQ=255)                                                tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTTACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                                >  1:668775/1‑139 (MQ=255)                                                tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                   >  2:827676/1‑104 (MQ=255)                                                tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                   <  1:827676/104‑1 (MQ=255)                                                tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                                >  1:683669/1‑139 (MQ=255)                                                tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                                >  1:686337/1‑139 (MQ=255)                                                tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTCTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                                >  2:681284/1‑139 (MQ=255)                                                 cTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                              >  1:899990/1‑92 (MQ=255)                                                 cTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                              <  2:899990/92‑1 (MQ=255)                                                 cTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCa                                                                               >  1:493123/1‑139 (MQ=255)                                                  tGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaa                                                                              >  2:464644/1‑139 (MQ=255)                                                  tGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaa                                                                              >  1:668298/1‑139 (MQ=255)                                                   ggTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaa                                                                             <  2:686337/139‑1 (MQ=255)                                                     tataCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                               >  2:921164/1‑87 (MQ=255)                                                     tataCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                               <  1:921164/87‑1 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                     <  2:998996/76‑1 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                     <  2:714353/76‑1 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                     >  1:714353/1‑76 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                     >  1:998996/1‑76 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                               >  2:159191/1‑82 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                               <  1:159191/82‑1 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAGCTACATTATTCTTACGATATATCAAAAATGAtt                                                                      >  1:855480/1‑139 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                      >  1:684214/1‑139 (MQ=255)                                                          cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                      >  2:824015/1‑139 (MQ=255)                                                              aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATg                                                                         <  1:714926/132‑1 (MQ=255)                                                              aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATg                                                                         >  2:714926/1‑132 (MQ=255)                                                              aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAcc                                                                  >  2:776830/1‑139 (MQ=255)                                                               gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCa                                                                 >  1:636285/1‑139 (MQ=255)                                                                    gCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttat                                                            <  1:776830/139‑1 (MQ=255)                                                                        ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                        <  1:464644/139‑1 (MQ=255)                                                                        ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                        <  2:668298/139‑1 (MQ=255)                                                                          cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                      <  2:493123/139‑1 (MQ=255)                                                                          cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                      <  2:684214/139‑1 (MQ=255)                                                                                       tCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaa                                                                              <  1:971800/102‑1 (MQ=255)                                                                                       tCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaa                                                                              >  2:971800/1‑102 (MQ=255)                                                                                           ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGtgtg                                     <  2:139683/139‑1 (MQ=255)                                                                                                 agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                   <  1:158140/55‑1 (MQ=255)                                                                                                 agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                   >  2:158140/1‑55 (MQ=255)                                                                                                 agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                                <  1:840625/74‑1 (MQ=255)                                                                                                 agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                                >  2:840625/1‑74 (MQ=255)                                                                                                   ccAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTaa                             >  2:741804/1‑139 (MQ=255)                                                                                                     agagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                        <  2:1051674/110‑1 (MQ=255)                                                                                                     agagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                        >  1:1051674/1‑110 (MQ=255)                                                                                                      gagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                                  >  2:169334/1‑115 (MQ=255)                                                                                                      gagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                                  <  1:169334/115‑1 (MQ=255)                                                                                                       agGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTAt                         >  1:756555/1‑139 (MQ=255)                                                                                                           cgacAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                                >  2:404029/1‑64 (MQ=255)                                                                                                           cgacAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                                <  1:404029/64‑1 (MQ=255)                                                                                                             acAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                   <  2:668775/139‑1 (MQ=255)                                                                                                             acAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                   <  2:383706/139‑1 (MQ=255)                                                                                                                 cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAc               >  2:318143/1‑139 (MQ=255)                                                                                                                       gCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAGCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg         <  2:855480/139‑1 (MQ=255)                                                                                                                        cTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGt        <  1:318143/139‑1 (MQ=255)                                                                                                                              tataAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCa  <  1:210206/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                       GCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCA  >  NC_000913/4541430‑4541693

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