Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,798	C→T	V179I (GTC→ATC)	ecpC ←	putative fimbrial usher protein EcpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,798	0	C	T	100.0%	214.9 / NA	64	V179I (GTC→ATC)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/17);  total (47/17)

ACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAG  >  NC_000913/308661‑308933                                                                                                                                          |                                                                                                                                        aCTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGc                                                                                                                                       >  2:732247/1‑140 (MQ=255)       cTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGt                                                                                                                                 >  2:459203/1‑140 (MQ=255)             aCCGCTACCGATCCCGTACAGCGAGCCGTTGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGc                                                                                                                           >  1:1866042/1‑140 (MQ=255)              ccGCTCCCGCCCCCGGCCGGCGCGCCGGCGGCCACCACGTGGTGTGCGAGCAGGGCATTAGCTTTATTCAGCAACAGAGGGGGGTGTGTTTTTCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCt                                                                                                                          <  2:631412/140‑1 (MQ=255)               cGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTtatctgtc                                                                                                                         >  1:1856182/1‑135 (MQ=255)               cGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCC‑AGTTATAGCTCa                                                                                                                        >  1:375668/1‑140 (MQ=255)                 cTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAg                                                                                                                       >  1:1129291/1‑140 (MQ=255)                 cTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAg                                                                                                                       >  1:1857256/1‑140 (MQ=255)                 cTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAg                                                                                                                       >  1:1718799/1‑140 (MQ=255)                 cTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGATAAAGATGCGCAAGTTATAGCTCAg                                                                                                                       >  1:668457/1‑140 (MQ=255)                                  cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGTCGTTACGCAACTGGTTGTTATAAATGCCCAAGTTATACCTCAGATTACTGCTGAGGgtgt                                                                                                      >  1:924256/1‑140 (MQ=255)                                  cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGgtgt                                                                                                      >  1:1824883/1‑140 (MQ=255)                                  cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGgtgt                                                                                                      >  1:1912015/1‑140 (MQ=255)                                  cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGgtgt                                                                                                      >  1:368400/1‑140 (MQ=255)                                  cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAAATGCCCAAGTTGTAGCTCAGAATACTGCTGAGGgtgt                                                                                                      >  1:506452/1‑140 (MQ=255)                                  cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACAGGTTGTTATAGATGCCAAAGTTATACCTCACATTACTGCTGAGGgtgt                                                                                                      >  2:584006/1‑140 (MQ=255)                                  cGAGCCGCCGACCACCACAGGATGTCCGCGCAGTGAAGTAACTTTATTCAGCGAAAGATGGCTGGTTGTGTTCCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGgtgt                                                                                                      <  2:1129291/140‑1 (MQ=255)                                    aGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTa                                                                                                    >  2:87036/1‑140 (MQ=255)                                                 cacaTGATGTCCGCGAAGTACAGTAACTTTATTCAGCAAAAGATAGCGGGTTTTTTTCCCCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGcc                                                                                       <  2:1718799/140‑1 (MQ=255)                                                     tgatgTCCGCGAGGTGCATTATTTTTATTCAGCGACAGATAGTTGGTTGTTTGGCCCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGAt                                                                                   <  2:1824883/140‑1 (MQ=255)                                                        tttgcgcgCAGTGAAGTACTTTTATTCAGCGAAAGATGGTGGGTTGTTTTGCCCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTc                                                                                <  2:1857256/136‑1 (MQ=255)                                                           tCGCGCAGTGCAGTAACGTTATTCAGTGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGAGTGGTAACACTGGACTGCCCGATGTCTTc                                                                             >  1:657803/1‑140 (MQ=255)                                                           tCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCTACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTc                                                                             >  1:1248087/1‑140 (MQ=255)                                                           tCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTc                                                                             >  1:298850/1‑140 (MQ=255)                                                           tCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTc                                                                             >  2:251418/1‑140 (MQ=255)                                                            cgcgCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCg                                                                            >  2:1341703/1‑140 (MQ=255)                                                            cgcgCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCg                                                                            >  2:1559470/1‑140 (MQ=255)                                                            cgcgAAGTGCATTAATTTTATCCAGCGAAAGATGGGTGGTTGTTTTCCCCCCGCCGTTACGCTACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCg                                                                            <  2:1248087/140‑1 (MQ=255)                                                              cgcAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGTTGAGGGTCTTAACATTGGACTACCCGATGTCTTCGCt                                                                          >  1:686274/1‑140 (MQ=255)                                                                  gtgcagtATTTTTATTCCGTGAAAGATGGTGGGTTTTTTTCCCCCCGCGGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGt                                                                      <  2:657803/140‑1 (MQ=255)                                                                  gtgcagtAACGTTATTCCGCGATAGATGGCTCGTTTTGTTGCTCCCGCTGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGt                                                                      <  2:506452/140‑1 (MQ=255)                                                                             ttATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTAc                                                           >  1:1159786/1‑140 (MQ=255)                                                                             ttATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAAGGTGTTAAAACTAGACGGCCCGATGTCTTCGCTACGTGAGCGTAGTAc                                                           >  2:752084/1‑140 (MQ=255)                                                                               aTTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGATGTTAACACTGGACTGGCCGATGTCGTCGCAACGTGAGCGTAGTACgg                                                         >  2:499248/1‑140 (MQ=255)                                                                                          gATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCt                                              >  1:463982/1‑140 (MQ=255)                                                                                             aGCTGGATGTGTTGCTGACGTCGTTACGCAGCTGGTTGTAATAGATGCTCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTcg                                           <  2:924256/140‑1 (MQ=255)                                                                                                tGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCgcgc                                        >  1:1311074/1‑140 (MQ=255)                                                                                                 ggATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCt                                       <  1:459203/140‑1 (MQ=255)                                                                                                   aggtgtTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTg                                     >  2:1367764/3‑140 (MQ=255)                                                                                                    tgtgTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGa                                    >  2:1794876/1‑140 (MQ=255)                                                                                                     gtgtTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGAc                                   >  1:250145/1‑140 (MQ=255)                                                                                                       gtTGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCa                                 >  2:126146/1‑140 (MQ=255)                                                                                                        ttGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCa                                 <  1:1562843/139‑1 (MQ=255)                                                                                                        ttGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCa                                 >  2:1562843/1‑139 (MQ=255)                                                                                                        ttGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCAc                                >  2:1850250/1‑140 (MQ=255)                                                                                                        ttGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCAc                                >  1:1233853/1‑140 (MQ=255)                                                                                                        ttGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCAc                                >  1:1876748/1‑140 (MQ=255)                                                                                                         tGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACt                               >  1:1019996/1‑140 (MQ=255)                                                                                                         tGCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACt                               <  1:1367764/140‑1 (MQ=255)                                                                                                          gCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTa                              >  1:1162518/1‑140 (MQ=255)                                                                                                          gCTCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAACGTAGTAGGGTACCTAGCGCTACGCGCTTGACCACTa                              >  1:582771/1‑140 (MQ=255)                                                                                                          ccTCCCGCCGTTACGCAACTGGCTGCTATAGATGCCCCAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTa                              <  1:732247/139‑1 (MQ=255)                                                                                                            tCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGc                            <  1:1559470/140‑1 (MQ=255)                                                                                                            tCCCGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGc                            >  1:1873444/1‑140 (MQ=255)                                                                                                            ccccGCCGTTACGCAACGGGTTTTTATAGATCCCCCAGTGATAGCTCAAATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGc                            <  2:587210/139‑1 (MQ=255)                                                                                                             cccGCCGTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCt                           >  1:1320276/1‑140 (MQ=255)                                                                                                                gccgTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTgca                        >  1:211377/1‑140 (MQ=255)                                                                                                                 ccgTTACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTgcag                       >  2:310815/1‑140 (MQ=255)                                                                                                                     tACGCAACTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGGGCCCAGCGCTTGGCGCTTGACAACTAGCTgcagcagc                   >  1:587210/1‑140 (MQ=255)                                                                                                                            cTGGTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGgc            >  2:1226921/1‑140 (MQ=255)                                                                                                                              ggTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGgcgc          <  2:1320276/140‑1 (MQ=255)                                                                                                                              ggTTGTTATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGgcgc          <  2:1233853/140‑1 (MQ=255)                                                                                                                                 tgttATAGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAg       >  2:1756832/1‑140 (MQ=255)                                                                                                                                      taGATGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTgag  >  1:20691/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        ACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAG  >  NC_000913/308661‑308933

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