Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,002	C→T	C65C (TGC→TGT)	fimH →	type 1 fimbriae D‑mannose specific adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,002	0	C	T	100.0%	214.1 / NA	60	C65C (TGC→TGT)	fimH	type 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/35);  total (25/35)

GTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACC  >  NC_000913/4548867‑4549128                                                                                                                                        |                                                                                                                               gTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTTTGTAAACCCTGCGCCCCTCGTGGATGTGGGGGAAAAACCGGTCGTGGGTCTTTTGACGCCAAACTTTTTTCATa                                                                                                                            >  2:490394/1‑140 (MQ=255)     ttCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGa                                                                                                                        <  1:176540/140‑1 (MQ=255)            gTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTATGTAAAACTTTCGCCCCTCGGGGATGTGGGGGAAAAACTGGGCGTGGATCTTTCGACGCAAAACTTTTTTCATAACGATTAttcc                                                                                                                 >  2:31377/1‑137 (MQ=255)            gTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCg                                                                                                                 <  1:490394/140‑1 (MQ=255)             tAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGGCCCCATGTTTTTTTAAAACATGAGCCCCCCGCGGCAGTGGGGGCAAACCCGGTGGCGGATGTTTCTCCGCAAAACTTTTTTCATTAAGAATATTCgg                                                                                                                >  2:1138610/1‑140 (MQ=255)              aaaaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGa                                                                                                               <  2:259775/140‑1 (MQ=255)                aaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGaaa                                                                                                             <  1:1183536/140‑1 (MQ=255)                aaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGaaa                                                                                                             <  1:461383/140‑1 (MQ=255)                 cccGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAAc                                                                                                            <  2:627632/139‑1 (MQ=255)                      cAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTa                                                                                                       <  2:660038/140‑1 (MQ=255)                       aaTGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCAAGTTTTTGTAAACCCTGCGCCCCTCGTGGATGTGGGGGAAAACCTCGGCGGGGATCTTTCGACGCAAATCTTTTTTCATAAAGAATATTCCGAAACCAtttc                                                                                                      >  2:630653/1‑138 (MQ=255)                        aTGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACa                                                                                                     <  1:1218652/140‑1 (MQ=255)                        aTGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTCCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACa                                                                                                     <  2:1074844/140‑1 (MQ=255)                          ggTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGa                                                                                                   <  2:490526/140‑1 (MQ=255)                          ggTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGa                                                                                                   <  1:34325/140‑1 (MQ=255)                          ggTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGa                                                                                                   <  1:2391/140‑1 (MQ=255)                           gTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGAc                                                                                                  <  1:532425/140‑1 (MQ=255)                               cGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATg                                                                                              <  1:179258/140‑1 (MQ=255)                                   aTCCCTATTGGCGGTGGGAGCGCCAATGTTTTTGTAAACCCTGCGCCCCTCGTGGATGTGGGGGAAAACCTGGTCGTGGGTCTTTTGACGCAAATCTTTTTTTATAACGATTATCCGGAAACCATTACAGACTAttttac                                                                                          >  2:164791/1‑135 (MQ=255)                                    tCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTTCGCCCCTCGTGGATGTGGGGGAAAACCTGGTCGTGGATCTTTTGACGCAAATCTTTTTTCATAAAGATTATCCGGAAAACATTACAGACTATGTcaca                                                                                         >  2:2391/1‑140 (MQ=255)                                        tATTGGCGGTGGCAGCGCCCATGTTTATTTAAACCTTGCGCCCGTCGTGGATGTGGGGCAAAAACTGGGTGTGGATCTTTTGACGCAAATCTTTTTTCATAACGATTATCCCGAAACCATTACAGACTATGTCACACTGc                                                                                     >  2:424286/1‑140 (MQ=255)                                        tATTGGCGGTGGCAGCGCCAATGTTTATGTAAAACCTTCGCCCGTCGTGAATGTGGGGCAAAACCCGGGCGTGGGTCTTTTGACGCAAAACTTTTTTCATAAAGATTATCCGGAAACCATTACAAACTATTTTACACTGc                                                                                     >  2:18789/1‑140 (MQ=255)                                         aTTGGCGGTGGGCGCGCCAAAGTTTTTTTAAAACTTTCGCCCCCCGTGGATGTGGGGGAAAACCTGGTCGGGGGTCTTTCTGCGCCAAACTTTTGTCATAAAGATTATTCCGAAACCATTACAGACTATGTCACACTGcc                                                                                    >  2:1173708/1‑139 (MQ=255)                                                gTGGCAGCGCCAATGTTTTTGTAAACCCTTCCCCCGTCGTGAATGTGGGGGAAAACCCGGTCGTGGGTCTTTTGACGCAAAACTTTTTTCATAACGATTATCCCGGAACCATTACAGACTATGTTACACTGCAAAGAggg                                                                             >  2:742259/1‑139 (MQ=255)                                                gTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGc                                                                             <  1:566718/140‑1 (MQ=255)                                                 tGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                            <  1:18789/140‑1 (MQ=255)                                                       cgcCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGc                                                                             <  2:878649/133‑1 (MQ=255)                                                       cgcCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGc                                                                             >  1:878649/1‑133 (MQ=255)                                                       cgcCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCtt                                                                      >  2:738813/1‑140 (MQ=255)                                                            aTGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTa                                                                     >  2:179612/1‑136 (MQ=255)                                                            aTGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTa                                                                     <  1:179612/136‑1 (MQ=255)                                                                       aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATc                                                      >  2:629537/1‑140 (MQ=255)                                                                        aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCt                                                     >  2:340844/1‑140 (MQ=255)                                                                        aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCt                                                     >  2:963542/1‑140 (MQ=255)                                                                         ccTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTa                                                    >  1:262508/1‑140 (MQ=255)                                                                           ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAt                                                  <  1:517913/140‑1 (MQ=255)                                                                            tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATATAAtt                                                 >  2:79223/1‑140 (MQ=255)                                                                             gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAttt                                                >  1:341293/1‑140 (MQ=255)                                                                             gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAttt                                                >  1:551588/1‑140 (MQ=255)                                                                             gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAttt                                                >  1:553151/1‑140 (MQ=255)                                                                             gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCCAAttt                                                >  1:581003/1‑140 (MQ=255)                                                                              cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAtttt                                               >  1:215582/1‑140 (MQ=255)                                                                                   gtcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                          <  1:638261/140‑1 (MQ=255)                                                                                   gtcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                          <  1:742259/140‑1 (MQ=255)                                                                                    tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGGTTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCggg                                         >  2:214238/1‑140 (MQ=255)                                                                                     cgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                        <  2:553151/140‑1 (MQ=255)                                                                                     cgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                        <  1:738813/140‑1 (MQ=255)                                                                                     cgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                        <  1:31377/140‑1 (MQ=255)                                                                                       tGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                      <  2:243218/140‑1 (MQ=255)                                                                                               gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                        >  2:737606/1‑130 (MQ=255)                                                                                               gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                        <  1:737606/130‑1 (MQ=255)                                                                                               gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAata                              >  2:213287/1‑140 (MQ=255)                                                                                               gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAata                              >  1:380107/1‑140 (MQ=255)                                                                                               gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGAAAAata                              >  2:942754/1‑140 (MQ=255)                                                                                                 gCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                            <  1:164791/140‑1 (MQ=255)                                                                                                 gCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                            >  1:995952/1‑140 (MQ=255)                                                                                                   aaaaCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGt                          >  1:876702/1‑140 (MQ=255)                                                                                                      aCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                       <  2:334973/140‑1 (MQ=255)                                                                                                      aCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                       <  1:214238/140‑1 (MQ=255)                                                                                                          ggTCGTGGATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTa                   <  1:309281/140‑1 (MQ=255)                                                                                                                ggATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc             >  1:671443/1‑140 (MQ=255)                                                                                                                 gATCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCCGCTTATCTATTTTTTTCGGGACCGTGAAATATAGGCGCAGTAgagacac            >  2:312763/1‑134 (MQ=255)                                                                                                                  aTCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa           <  2:215582/140‑1 (MQ=255)                                                                                                                   tCTTTCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt          <  2:995952/140‑1 (MQ=255)                                                                                                                      ttCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGc                 >  1:970502/1‑130 (MQ=255)                                                                                                                      ttCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGc                 <  2:970502/130‑1 (MQ=255)                                                                                                                       tCGACGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTcc      >  2:151155/1‑140 (MQ=255)                                                                                                                          aCGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTac   <  2:380107/140‑1 (MQ=255)                                                                                                                          aCGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTac   >  1:86508/1‑140 (MQ=255)                                                                                                                           cGCAAATCTTTTGTCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc  <  1:312763/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                               GTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACC  >  NC_000913/4548867‑4549128

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