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%	238.3 / NA	73	*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 (37/36);  total (37/36)

ATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAC  >  NC_000913/4541426‑4541691                                                                                                                                      |                                                                                                                                     aTACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                   <  1:73927/137‑1 (MQ=255) aTACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                   >  2:73927/1‑137 (MQ=255) aTACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATa                                                                                                                                 >  1:1207323/1‑139 (MQ=255)     gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             >  2:946334/1‑139 (MQ=255)     gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             >  1:115113/1‑139 (MQ=255)     gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             <  1:946334/139‑1 (MQ=255)       gACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                           >  1:767133/1‑139 (MQ=255)        aCTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgaga                                                                                                                          <  1:185594/139‑1 (MQ=255)         cTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGAt                                                                                                                         <  1:1505609/139‑1 (MQ=255)                 ggATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCATCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 <  1:1047813/139‑1 (MQ=255)                 ggATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 <  2:1392339/139‑1 (MQ=255)                 ggATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 <  2:1207323/139‑1 (MQ=255)                      aTCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGgtg                                                                                                            >  2:373612/1‑139 (MQ=255)                       tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGt                                                                                                                               >  2:1546422/1‑119 (MQ=255)                       tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGt                                                                                                                               <  1:1546422/119‑1 (MQ=255)                       tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGgtgt                                                                                                           >  2:500882/1‑139 (MQ=255)                         ttGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTa                                                                                                         >  1:1527476/1‑139 (MQ=255)                           gggCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaa                                                                                                       <  1:373612/139‑1 (MQ=255)                            ggCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaa                                                                                                      >  2:1613083/1‑139 (MQ=255)                             gCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                     <  2:767133/139‑1 (MQ=255)                                tCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             <  1:1010835/112‑1 (MQ=255)                                tCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             >  2:1010835/1‑112 (MQ=255)                                 cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTAc                                                                                                 >  2:333578/1‑139 (MQ=255)                                 cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTAc                                                                                                 >  1:1261547/1‑139 (MQ=255)                                   cAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACat                                                                                               <  1:333578/139‑1 (MQ=255)                                           gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCt                                                                                         <  2:219032/137‑1 (MQ=255)                                           gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCt                                                                                         >  1:453525/1‑137 (MQ=255)                                           gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCt                                                                                         <  2:453525/137‑1 (MQ=255)                                           gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCt                                                                                         >  1:219032/1‑137 (MQ=255)                                                  tGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGatata                                                                                >  2:1460987/1‑139 (MQ=255)                                                    tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                              >  2:1387268/1‑139 (MQ=255)                                                    tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                              >  1:454523/1‑139 (MQ=255)                                                    tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                              >  1:921375/1‑139 (MQ=255)                                                       ggTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaa                                                                           <  1:1375603/139‑1 (MQ=255)                                                        gTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaaa                                                                          >  1:239491/1‑139 (MQ=255)                                                           taCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                                  <  1:1114000/80‑1 (MQ=255)                                                           taCCGCCAGCAATCCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                                  >  2:1114000/1‑80 (MQ=255)                                                             ccGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAt                                                                     <  2:1261547/139‑1 (MQ=255)                                                              cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                    <  1:918369/139‑1 (MQ=255)                                                              cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                    >  2:1264824/1‑139 (MQ=255)                                                              cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                    <  2:921375/139‑1 (MQ=255)                                                              cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                    >  1:1614241/1‑139 (MQ=255)                                                                 cAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAc                                                                 <  2:1527476/139‑1 (MQ=255)                                                                  aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAcc                                                                >  2:804955/1‑139 (MQ=255)                                                                   gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                                  >  2:527622/1‑72 (MQ=255)                                                                   gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                                  <  1:527622/72‑1 (MQ=255)                                                                   gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgaga                                                                                                                          >  1:215053/1‑80 (MQ=255)                                                                   gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgaga                                                                                                                          <  2:215053/80‑1 (MQ=255)                                                                   gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCa                                                               >  2:1294733/1‑139 (MQ=255)                                                                             gCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                     >  1:269716/1‑91 (MQ=255)                                                                             gCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                     <  2:269716/91‑1 (MQ=255)                                                                              cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACat                                                                                               <  2:433091/96‑1 (MQ=255)                                                                              cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACat                                                                                               >  1:433091/1‑96 (MQ=255)                                                                              cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                    >  2:1418359/1‑139 (MQ=255)                                                                              cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                    <  2:454523/139‑1 (MQ=255)                                                                                  tttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                                <  2:239491/139‑1 (MQ=255)                                                                                            cTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgt                                      <  1:1294733/139‑1 (MQ=255)                                                                                                     agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGt                                                                                                                               >  2:955558/1‑41 (MQ=255)                                                                                                     agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGt                                                                                                                               <  1:955558/41‑1 (MQ=255)                                                                                                     agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCtt                             >  2:178702/1‑139 (MQ=255)                                                                                                     agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCtt                             >  1:336542/1‑139 (MQ=255)                                                                                                     agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCtt                             >  2:1147059/1‑139 (MQ=255)                                                                                                       ccAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATAcc                                                                                                                      <  1:37740/48‑1 (MQ=255)                                                                                                       ccAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATAcc                                                                                                                      >  2:37740/1‑48 (MQ=255)                                                                                                       ccAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGat                                                                                   >  2:91767/1‑83 (MQ=255)                                                                                                       ccAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGat                                                                                   <  1:91767/83‑1 (MQ=255)                                                                                                        cAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAAt                          <  1:799044/139‑1 (MQ=255)                                                                                                         agagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAgggg                                                                                                                 >  1:340405/1‑50 (MQ=39)                                                                                                         agagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 <  2:340405/51‑1 (MQ=255)                                                                                                          gagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                        <  2:1614241/139‑1 (MQ=255)                                                                                                          gagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                        <  2:336542/139‑1 (MQ=255)                                                                                                                               ttttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGat                                                                                   >  2:1322371/1‑59 (MQ=255)                                                                                                                               ttttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGat                                                                                   <  1:1322371/59‑1 (MQ=255)                                                                                                                                tttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAc  <  1:1147059/139‑1 (MQ=255)                                                                                                                                tttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAc  <  1:1460987/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     ATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAC  >  NC_000913/4541426‑4541691

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