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%	305.1 / NA	93	*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 (44/49);  total (44/49)

AGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAC  >  NC_000913/4541424‑4541691                                                                                                                                        |                                                                                                                                     aGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                   <  1:187004/139‑1 (MQ=255)       gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             >  2:1073580/1‑139 (MQ=255)        cgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            <  1:1073580/139‑1 (MQ=255)        cgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            >  2:724744/1‑139 (MQ=255)         gACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                           <  1:724744/139‑1 (MQ=255)                ccAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAg                                                                                                                    >  2:75942/1‑139 (MQ=255)                   ggATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 <  2:130616/139‑1 (MQ=255)                     aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 >  1:1302960/1‑137 (MQ=255)                     aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 <  2:1302960/137‑1 (MQ=255)                     aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                               >  1:521550/1‑139 (MQ=255)                         tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGgtgt                                                                                                           >  1:1239822/1‑139 (MQ=255)                          cTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGtt                                                                                                          <  1:938578/139‑1 (MQ=255)                              ggCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaa                                                                                                      >  2:1523/1‑139 (MQ=255)                                cATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGtt                                                                                                          <  2:397253/133‑1 (MQ=255)                                cATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGtt                                                                                                          >  1:397253/1‑133 (MQ=255)                                   cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTAc                                                                                                 >  2:15599/1‑139 (MQ=255)                                   cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTAc                                                                                                 >  1:298720/1‑139 (MQ=255)                                    gCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                                <  1:909738/139‑1 (MQ=255)                                      aaTATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                              <  1:556411/139‑1 (MQ=255)                                       atatTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatta                                                                                             <  2:400650/139‑1 (MQ=255)                                        tatTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAACACTACattat                                                                                            <  2:176967/139‑1 (MQ=255)                                            cGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCtt                                                                                        <  1:367827/139‑1 (MQ=255)                                             gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTa                                                                                       >  2:205439/1‑139 (MQ=255)                                             gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTa                                                                                       <  1:15599/139‑1 (MQ=255)                                                    tGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATgg                                                                                                              >  1:604718/1‑109 (MQ=255)                                                    tGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATgg                                                                                                              <  2:604718/109‑1 (MQ=255)                                                    tGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGatata                                                                                >  2:561764/1‑139 (MQ=255)                                                      tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                              >  2:746321/1‑139 (MQ=255)                                                       cTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCa                                                                             >  2:687654/1‑139 (MQ=255)                                                        tGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            >  2:651409/1‑91 (MQ=255)                                                        tGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            <  1:651409/91‑1 (MQ=255)                                                         ggTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaa                                                                           <  2:521550/139‑1 (MQ=255)                                                          gTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaaa                                                                          <  1:75942/139‑1 (MQ=255)                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                           <  1:490627/84‑1 (MQ=255)                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                           >  2:490627/1‑84 (MQ=255)                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                    >  2:70572/1‑139 (MQ=255)                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                    >  1:278344/1‑139 (MQ=255)                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACACTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCt                                                                                         <  1:896121/118‑1 (MQ=255)                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACACTGATTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCt                                                                                         >  2:896121/1‑118 (MQ=255)                                                                  ccAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTaa                                                                  <  1:1523/139‑1 (MQ=255)                                                                   cAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAc                                                                 <  1:746321/139‑1 (MQ=255)                                                                    aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAt                                                                     >  2:1283212/1‑134 (MQ=255)                                                                    aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAt                                                                     <  1:1283212/134‑1 (MQ=255)                                                                     gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             <  2:423085/77‑1 (MQ=255)                                                                     gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                             >  1:423085/1‑77 (MQ=255)                                                                        aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCAttt                                                            <  2:278344/139‑1 (MQ=255)                                                                         tGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                   <  1:718957/67‑1 (MQ=255)                                                                         tGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                                   >  2:718957/1‑67 (MQ=255)                                                                         tGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATtta                                                           <  1:687654/139‑1 (MQ=255)                                                                               gCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            >  2:1257094/1‑68 (MQ=255)                                                                               gCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            <  1:1257094/68‑1 (MQ=255)                                                                                cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            <  2:1220502/67‑1 (MQ=255)                                                                                cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                            >  1:1220502/1‑67 (MQ=255)                                                                                cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaa                                                                                                       <  2:1129466/88‑1 (MQ=255)                                                                                cGTTTTTACGGCATCGGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaa                                                                                                       >  1:1129466/1‑88 (MQ=255)                                                                                   ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAt                                                                     <  2:423782/119‑1 (MQ=255)                                                                                   ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAt                                                                     >  1:423782/1‑119 (MQ=255)                                                                                   ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                 <  2:298720/139‑1 (MQ=255)                                                                                   ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                 >  1:796895/1‑139 (MQ=255)                                                                                    tttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGATGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                                <  2:796895/139‑1 (MQ=255)                                                                                        cGGCATCTGGGATAGAGCCAGAGGACGAGAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                     <  1:334452/82‑1 (MQ=255)                                                                                        cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                     >  2:334452/1‑82 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGAt                                                                                                                         >  2:1172321/1‑47 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGAt                                                                                                                         <  1:1172321/47‑1 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATAc                                                                                                                       >  1:387908/1‑49 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATAc                                                                                                                       <  2:387908/49‑1 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 <  1:386864/55‑1 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                                 >  2:386864/1‑55 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                              <  1:566958/74‑1 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                                <  1:1253293/72‑1 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                                >  2:1253293/1‑72 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAc                                                                 <  1:288240/103‑1 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAc                                                                 >  2:288240/1‑103 (MQ=255)                                                                                                       agCCAGAGGACGACAGCGCCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                              >  2:566958/1‑74 (MQ=255)                                                                                                           agagGACGACAGCGTCACGCTGTTTTATAAAGATAGtt                                                                                                                              >  2:945381/1‑38 (MQ=255)                                                                                                           agagGACGACAGCGTCACGCTGTTTTATAAAGATAGtt                                                                                                                              <  1:945381/38‑1 (MQ=255)                                                                                                            gagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                              <  1:1172070/69‑1 (MQ=255)                                                                                                            gagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACatt                                                                                              >  2:1172070/1‑69 (MQ=255)                                                                                                             agGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTa                                              >  2:1263357/1‑116 (MQ=255)                                                                                                             agGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTa                                              <  1:1263357/116‑1 (MQ=255)                                                                                                             agGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTAt                       >  2:717919/1‑139 (MQ=255)                                                                                                                   acAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                 >  2:1078780/1‑139 (MQ=255)                                                                                                                   acAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                 <  1:70572/139‑1 (MQ=255)                                                                                                                       cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                               <  2:475205/41‑1 (MQ=255)                                                                                                                       cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                               >  1:475205/1‑41 (MQ=255)                                                                                                                       cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                     <  2:458445/51‑1 (MQ=255)                                                                                                                       cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                     >  1:458445/1‑51 (MQ=255)                                                                                                                       cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                      <  2:377813/98‑1 (MQ=255)                                                                                                                       cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                      >  1:377813/1‑98 (MQ=255)                                                                                                                       cGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAc             >  2:40128/1‑139 (MQ=255)                                                                                                                                 ttttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAtt          <  1:52888/132‑1 (MQ=255)                                                                                                                                 ttttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAtt          >  2:52888/1‑132 (MQ=255)                                                                                                                                  tttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAc  <  1:566012/139‑1 (MQ=255)                                                                                                                                  tttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAc  <  2:1239822/139‑1 (MQ=255)                                                                                                                                  tttATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAc  <  1:1078780/139‑1 (MQ=255)                                                                                                                                  tttATAAAGATAGATGAGATACCAGGGATGGTGTTAAAAACTACATTATGCTTACGAGATATCAAAAATGACTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATCAAGTTAAc  <  1:352450/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     AGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAC  >  NC_000913/4541424‑4541691

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: ‑