Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,061	A→G	Y85C (TAT→TGT)	fimH →	minor component of type 1 fimbriae

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,061	0	A	G	100.0%	292.2 / NA	81	Y85C (TAT→TGT)	fimH	minor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/39);  total (42/39)

TGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTA  >  NC_000913/4548927‑4549193                                                                                                                                       |                                                                                                                                     tGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggc                                                                                                                                  >  2:251691/1‑139 (MQ=255)  gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 <  1:251691/139‑1 (MQ=255)  gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 <  1:1666320/139‑1 (MQ=255)      aTGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgt                                                                                                                             <  2:1813328/139‑1 (MQ=255)       tctAAACCTAGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACTCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGATGCTCGGCTTGTGGCGGCgtg                                                                                                                            <  1:1211853/137‑1 (MQ=255)                tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAAt                                                                                                                   >  2:240635/1‑139 (MQ=255)                 gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAAtt                                                                                                                  >  1:24395/1‑139 (MQ=255)                 gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAAtt                                                                                                                  >  2:449034/1‑139 (MQ=255)                 gcgcCCGTCGTGAATGT‑GGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggc                                                                                                                                  >  2:1856453/1‑122 (MQ=255)                 gcgcCCGTCGTGAATGT‑GGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggc                                                                                                                                  <  1:1856453/122‑1 (MQ=255)                       gtcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 >  1:914848/1‑118 (MQ=255)                       gtcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 <  2:914848/118‑1 (MQ=255)                        tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                           <  1:1522333/139‑1 (MQ=255)                         cgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcggc                                                                                                                               >  2:1253510/1‑118 (MQ=255)                         cgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcggc                                                                                                                               <  1:1253510/118‑1 (MQ=255)                          tttaaTGTGGGTCAAAACCTTGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         <  1:240635/136‑1 (MQ=255)                          gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         <  1:1559958/139‑1 (MQ=255)                          gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         <  2:682532/139‑1 (MQ=255)                          gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         <  2:24395/139‑1 (MQ=255)                           tGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGAc                                                                                                        >  1:1990244/1‑139 (MQ=255)                                   gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 <  1:771387/106‑1 (MQ=255)                                   gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 >  2:771387/1‑106 (MQ=255)                                   gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgt                                                                                                                             >  1:1601762/1‑110 (MQ=255)                                   gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgt                                                                                                                             <  2:1601762/110‑1 (MQ=255)                                   gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                                                                                >  1:1992649/1‑139 (MQ=255)                                           ccTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                                        <  1:449034/139‑1 (MQ=255)                                                              aCGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgtg                                                                                                                            <  2:668678/84‑1 (MQ=255)                                                              aCGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgtg                                                                                                                            >  1:668678/1‑84 (MQ=255)                                                              aCGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTa                                                                     >  1:1735400/1‑139 (MQ=255)                                                                          tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         <  1:1465218/91‑1 (MQ=255)                                                                          tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         >  2:1465218/1‑91 (MQ=255)                                                                                   gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                                <  2:1992649/139‑1 (MQ=255)                                                                                   gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                                <  2:523918/139‑1 (MQ=255)                                                                                       aTCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCa                                                                                       >  2:1288102/1‑96 (MQ=255)                                                                                       aTCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCa                                                                                       <  1:1288102/96‑1 (MQ=255)                                                                                             aaaCCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                           <  1:508264/102‑1 (MQ=255)                                                                                             aaaCCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                           >  2:508264/1‑102 (MQ=255)                                                                                                  aTTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 >  2:1681565/1‑43 (MQ=255)                                                                                                  aTTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                                 <  1:1681565/43‑1 (MQ=255)                                                                                                    tACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgt                                               >  1:87431/1‑123 (MQ=255)                                                                                                    tACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgt                                               <  2:87431/123‑1 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTc                                                                                                              <  1:67486/55‑1 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTc                                                                                                              >  2:67486/1‑55 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                           <  2:913103/58‑1 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                           >  1:913103/1‑58 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                            <  1:1789303/89‑1 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                            >  2:1789303/1‑89 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGa                          >  1:770753/1‑139 (MQ=255)                                                                                                         aCTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGa                          <  2:1735400/139‑1 (MQ=255)                                                                                                           tATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATa                        >  1:758805/1‑139 (MQ=255)                                                                                                             tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg                      >  2:1812952/1‑139 (MQ=255)                                                                                                             tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg                      >  1:1428663/1‑139 (MQ=255)                                                                                                             tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg                      >  2:853875/1‑139 (MQ=255)                                                                                                                  cacTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         >  2:1314752/1‑51 (MQ=255)                                                                                                                  cacTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         <  1:1314752/51‑1 (MQ=255)                                                                                                                  cacTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTa                                                                                                    >  2:733921/1‑56 (MQ=255)                                                                                                                  cacTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTa                                                                                                    <  1:733921/56‑1 (MQ=255)                                                                                                                   acTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                                <  2:1443744/107‑1 (MQ=255)                                                                                                                   acTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                                >  1:1443744/1‑107 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                           <  1:891264/46‑1 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                           >  2:891264/1‑46 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                             <  2:233101/60‑1 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                             >  1:233101/1‑60 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                           >  2:937366/1‑78 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                           <  1:937366/78‑1 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgaga                               <  1:1923692/122‑1 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgaga                               >  2:1923692/1‑122 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGAt                         <  2:545365/128‑1 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGAt                         >  1:545365/1‑128 (MQ=255)                                                                                                                     tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc              >  1:820060/1‑139 (MQ=255)                                                                                                                       cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCg                                                                                                      <  1:1971818/49‑1 (MQ=255)                                                                                                                       cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCg                                                                                                      >  2:1971818/1‑49 (MQ=255)                                                                                                                       cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg            >  2:147692/1‑139 (MQ=255)                                                                                                                        aaCGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGt           <  1:147692/139‑1 (MQ=255)                                                                                                                               cTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCTTTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgag                                <  1:1458566/111‑1 (MQ=255)                                                                                                                               cTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgag                                >  2:1458566/1‑111 (MQ=255)                                                                                                                               cTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCtt    >  1:1350463/1‑139 (MQ=255)                                                                                                                                 cGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         >  2:336235/1‑36 (MQ=255)                                                                                                                                 cGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                         <  1:336235/36‑1 (MQ=38)                                                                                                                                 cGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTa  >  1:239378/1‑139 (MQ=255)                                                                                                                                 cGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTa  >  2:1168857/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                     TGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTA  >  NC_000913/4548927‑4549193

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