Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,540,984	(A)7→6	coding (28/603 nt)	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,540,978	0	A	.	100.0%	327.3 / NA	78	coding (22/603 nt)	fimB	tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (42/36);  total (42/36)

TTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGT  >  NC_000913/4540848‑4541112                                                                                                                                   |                                                                                                                                       ttagttaTTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaagga                                                                                                                                <  2:1256055/139‑4 (MQ=255)     ttattaAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACtt                                                                                                                           >  1:1005495/1‑139 (MQ=255)     ttattaAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACtt                                                                                                                           >  1:589531/1‑139 (MQ=255)      tattaAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTc                                                                                                                          <  2:143768/139‑1 (MQ=255)      tattaAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTc                                                                                                                          <  2:589531/139‑1 (MQ=255)           aGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGAc                                                                                                                     >  1:498970/1‑139 (MQ=255)                aCGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATGAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATa                                                                                                                >  1:794511/1‑139 (MQ=255)                              aaaCATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  <  2:1005495/139‑1 (MQ=255)                                  aTAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTcc                                                                                              >  2:942045/1‑139 (MQ=255)                                    aaTCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGt                                                                                                              <  2:910400/121‑1 (MQ=255)                                    aaTCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGt                                                                                                              >  1:910400/1‑121 (MQ=255)                                    aaTCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGa                                                                                                            >  1:731936/1‑123 (MQ=255)                                    aaTCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGa                                                                                                            <  2:731936/123‑1 (MQ=255)                                       cAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTaaa                                                                                         <  2:125933/139‑1 (MQ=255)                                        aGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAag                                                                                        <  1:942045/139‑1 (MQ=255)                                          gATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  <  1:909780/127‑1 (MQ=255)                                          gATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  <  2:1181201/127‑1 (MQ=255)                                          gATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  >  1:1181201/1‑127 (MQ=255)                                          gATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  >  2:909780/1‑127 (MQ=255)                                                    ttGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCg                                                                            >  1:36517/1‑139 (MQ=255)                                                     tGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  <  2:498308/116‑1 (MQ=255)                                                     tGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  >  1:498308/1‑116 (MQ=255)                                                     tGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCgg                                                                           <  2:794511/139‑1 (MQ=255)                                                         ttattGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCt                                                                       >  2:425580/1‑139 (MQ=255)                                                                tAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcag                                                                <  1:229753/139‑1 (MQ=255)                                                                tAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcag                                                                <  2:777059/139‑1 (MQ=255)                                                                 aaCCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  >  2:1047530/1‑104 (MQ=255)                                                                 aaCCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  <  1:1047530/104‑1 (MQ=255)                                                                    ccAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACg                                                            <  2:36517/139‑1 (MQ=255)                                                                      aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaaggaa                                                                                                                               >  1:410964/1‑67 (MQ=255)                                                                      aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaaggaa                                                                                                                               <  2:410964/71‑5 (MQ=255)                                                                      aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    >  1:50701/1‑97 (MQ=255)                                                                      aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    <  2:50701/97‑1 (MQ=255)                                                                      aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTa                                                          >  1:1200680/1‑139 (MQ=255)                                                                      aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGC‑‑ATACCGGGCCTCATGCAGCACGTAat                                                        >  2:530701/1‑139 (MQ=255)                                                                        cacaGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaa                                                                                                                                    <  2:33400/64‑1 (MQ=255)                                                                        cacaGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaa                                                                                                                                    >  1:33400/1‑64 (MQ=255)                                                                             cTAGTGTGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑ACAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTg                                                   >  2:1131434/1‑139 (MQ=255)                                                                             cTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    <  2:97898/90‑1 (MQ=255)                                                                             cTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    >  1:97898/1‑90 (MQ=255)                                                                             cTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTg                                                   <  1:886562/139‑1 (MQ=255)                                                                              tAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAgcagca                                                                                   >  2:721012/1‑106 (MQ=255)                                                                              tAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAgcagca                                                                                   <  1:721012/106‑1 (MQ=255)                                                                                 tGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTg                                               >  1:1208383/1‑139 (MQ=255)                                                                                     cgTCTGTAATTATAAGGGAAAAACGATGAAGAATGAGGCTGATAAC‑AAACAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATAcc                                                                             >  2:850243/1‑105 (MQ=255)                                                                                     cgTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATGGAATCACTCCTTAAAGCAGCAAATAcc                                                                             <  1:850243/105‑1 (MQ=255)                                                                                         tGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGActttgct                                       <  2:179712/139‑1 (MQ=255)                                                                                           tAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTaaa                                                                                         <  2:463358/87‑1 (MQ=255)                                                                                           tAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTaaa                                                                                         >  1:463358/1‑87 (MQ=255)                                                                                           tAATTATAAGGG‑AAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGActttgctttg                                    >  2:1177397/1‑139 (MQ=255)                                                                                             aTTATAAGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGt                                   <  2:1208383/139‑1 (MQ=255)                                                                                                  aaGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAag                                                                                        >  2:235306/1‑81 (MQ=255)                                                                                                  aaGGGAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAag                                                                                        <  1:235306/81‑1 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    >  2:1289150/1‑66 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    <  1:1289150/66‑1 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  <  1:972276/68‑1 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  >  2:972276/1‑68 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcag                                                                >  2:744668/1‑102 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGActttgct                                       >  1:1059844/1‑127 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGActttgct                                       <  2:1059844/127‑1 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCa                           <  1:972020/139‑1 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCa                           >  1:235245/1‑139 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCa                           <  2:498970/139‑1 (MQ=255)                                                                                                     ggAAAAACGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACGTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcag                                                                <  1:744668/102‑1 (MQ=255)                                                                                                         aaaCGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcagca                                                              <  1:413485/100‑1 (MQ=255)                                                                                                         aaaCGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcagca                                                              >  2:413485/1‑100 (MQ=255)                                                                                                         aaaCGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGt                       >  1:672348/1‑139 (MQ=255)                                                                                                         aaaCGATGAAGAATAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGt                       >  1:242543/1‑139 (MQ=255)                                                                                                                 aagaaTAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGc               >  2:130276/1‑139 (MQ=255)                                                                                                                    aaTAAGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAg            >  2:1087083/1‑139 (MQ=255)                                                                                                                        aGGCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGaa        >  1:25161/1‑139 (MQ=255)                                                                                                                         ggCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGtt                                  <  2:401134/112‑1 (MQ=255)                                                                                                                         ggCTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGtt                                  >  1:401134/1‑112 (MQ=255)                                                                                                                           cTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    <  2:1053579/44‑1 (MQ=255)                                                                                                                           cTGATAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                                    >  1:1053579/1‑44 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  <  2:444603/43‑1 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                                  >  1:444603/1‑43 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATAcc                                                                             >  2:1012136/1‑64 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATAcc                                                                             <  1:1012136/64‑1 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGttt                                 >  2:202107/1‑108 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGttt                                 <  1:202107/108‑1 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGTTTTCCGGGCGAGTGAAATTTGt  >  2:1284449/1‑139 (MQ=255)                                                                                                                              aTAAC‑AAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGt  >  1:1234591/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                       TTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGT  >  NC_000913/4540848‑4541112

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