Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	562,993	C→T	P220S (CCT→TCT)	sfmH →	FimA homolog, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	562,993	0	C	T	100.0%	350.1 / NA	92	P220S (CCT→TCT)	sfmH	FimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/48);  total (44/48)

TTGGTTACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGT  >  NC_000913/562849‑563138                                                                                                                                                 |                                                                                                                                                  ttGGTTACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATctgctg                                                                                                                                             >  1:54101/1‑151 (MQ=255)     ttACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGgtgt                                                                                                                                         >  4:858329/1‑151 (MQ=255)          ccTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATg                                                                                                                                     >  2:759388/1‑150 (MQ=255)             gTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGc                                                                                                                                 >  4:157081/1‑151 (MQ=255)             gTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGc                                                                                                                                 >  3:360388/1‑151 (MQ=255)             gTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGc                                                                                                                                 <  2:1590100/151‑1 (MQ=255)               ttATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCATATTGTTGAATTTGATTTTGGTGATAGCGGCGCATCGTTATTTAGTGCGGCAGGGGCAGGTAATCGATCTGCTGGTGTCATGCCGCa                                                                                                                                <  2:89298/150‑1 (MQ=255)               ttATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCa                                                                                                                                >  4:720197/1‑150 (MQ=255)               ttATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCa                                                                                                                                >  3:265460/1‑150 (MQ=255)                tataCCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCaa                                                                                                                               >  4:1059608/1‑150 (MQ=255)                 ataCCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCaaa                                                                                                                              <  4:1110134/150‑1 (MQ=255)                  taCCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAAcc                                                                                                                            >  4:723354/1‑151 (MQ=255)                     cATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAag                                                                                                                         <  4:1175442/151‑1 (MQ=255)                      aTTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAaga                                                                                                                        <  2:495995/151‑1 (MQ=255)                        tAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGc                                                                                                                      >  2:1216516/1‑151 (MQ=255)                        tAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGGGGCGTGCCGCAAACCAAGAGc                                                                                                                      >  3:26170/1‑151 (MQ=255)                         aGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGc                                                                                                                      <  3:180334/150‑1 (MQ=255)                            tttGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCAtt                                                                                                                   <  2:1525730/150‑1 (MQ=255)                            tttGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCAtt                                                                                                                   <  2:307780/150‑1 (MQ=255)                             ttGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTg                                                                                                                  <  2:224057/150‑1 (MQ=255)                                  ggCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                                             <  1:196992/150‑1 (MQ=255)                                  ggCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCa                                                                                                            <  2:1074997/151‑1 (MQ=255)                                     cGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCaaa                                                                                                          >  4:1604815/1‑150 (MQ=255)                                     cGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAAt                                                                                                         <  3:333164/151‑1 (MQ=255)                                      gggTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAAt                                                                                                         >  2:142240/1‑150 (MQ=255)                                       ggTTGAAGTTCCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGt                                                                                                       <  2:1593568/151‑1 (MQ=255)                                       ggTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGt                                                                                                       <  4:84904/151‑1 (MQ=255)                                        gTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTa                                                                                                      <  2:112272/151‑1 (MQ=255)                                          tGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACg                                                                                                    >  2:37053/1‑151 (MQ=255)                                           gAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACg                                                                                                    >  4:592680/1‑150 (MQ=255)                                           gAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGa                                                                                                   <  2:946325/151‑1 (MQ=255)                                            aaGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGaa                                                                                                  <  2:368511/151‑1 (MQ=255)                                              gTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATg                                                                                                <  4:1259934/151‑1 (MQ=255)                                              gTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATg                                                                                                <  2:903370/151‑1 (MQ=255)                                                aCCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGtt                                                                                              >  2:1568975/1‑151 (MQ=255)                                                aCCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGtt                                                                                              >  2:1033555/1‑151 (MQ=255)                                                 ccGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTtg                                                                                             <  2:784487/151‑1 (MQ=255)                                                  cGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTtgc                                                                                            >  4:220285/1‑151 (MQ=255)                                                   gCAAAACTGCGAGTTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGGAATCGATCTGCTGGGGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTtgc                                                                                            >  4:737452/1‑150 (MQ=255)                                                      aaaCTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTgcg                                                                                        >  3:4284/1‑151 (MQ=255)                                                      aaaCTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTgcg                                                                                        >  3:367556/1‑151 (MQ=255)                                                       aaCTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTgcgc                                                                                       >  4:582057/1‑151 (MQ=255)                                                        aCTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCa                                                                                      <  4:932211/151‑1 (MQ=255)                                                        aCTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCa                                                                                      <  2:1403673/151‑1 (MQ=255)                                                         cTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTgcgc                                                                                       >  2:613075/1‑149 (MQ=255)                                                          tGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAgg                                                                                    >  2:471206/1‑151 (MQ=255)                                                            cGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCt                                                                                  >  3:69510/1‑151 (MQ=255)                                                             gAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCtt                                                                                 <  2:1443914/151‑1 (MQ=255)                                                               aTTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTAt                                                                               <  4:473939/151‑1 (MQ=255)                                                                    aTGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTaa                                                                           >  2:649049/1‑150 (MQ=255)                                                                      gCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACa                                                                         >  4:77508/1‑150 (MQ=255)                                                                      gCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACa                                                                         <  2:1006431/150‑1 (MQ=255)                                                                        cGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCTTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTa                                                                            <  1:408785/145‑1 (MQ=255)                                                                        cGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTa                                                                            >  3:408802/1‑145 (MQ=255)                                                                          ggCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCg                                                                    >  2:631611/1‑151 (MQ=255)                                                                           gCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGt                                                                   <  2:742677/151‑1 (MQ=255)                                                                               aTTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTgg                                                               <  4:885204/151‑1 (MQ=255)                                                                                 tgttgAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGa                                                              <  2:27249/150‑1 (MQ=255)                                                                                  gttgAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGa                                                              >  4:1360770/1‑149 (MQ=255)                                                                                  gttgAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAg                                                            <  4:1527904/151‑1 (MQ=255)                                                                                     gAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGcc                                                          <  1:172791/150‑1 (MQ=255)                                                                                      aaTTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGcc                                                          >  4:1469421/1‑149 (MQ=255)                                                                                      aaTTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGcc                                                          >  2:394106/1‑149 (MQ=255)                                                                                         ttGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTg                                                      <  4:174496/150‑1 (MQ=255)                                                                                           gATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATTCGTCTGGAAGCCAGTGcc                                                    >  4:1251674/1‑150 (MQ=255)                                                                                            aTTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGt                                                  >  2:1029740/1‑151 (MQ=255)                                                                                            aTTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGt                                                  >  2:1043155/1‑151 (MQ=255)                                                                                             ttttGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGtt                                                 >  4:1473655/1‑151 (MQ=255)                                                                                                  gTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTgg                                            <  2:402856/151‑1 (MQ=255)                                                                                                    gATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGt                                           <  2:212001/150‑1 (MQ=255)                                                                                                      tatCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTc                                          >  4:144999/1‑149 (MQ=255)                                                                                                      tatCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTc                                          >  4:739393/1‑149 (MQ=255)                                                                                                         cGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGc                                      <  4:806961/150‑1 (MQ=255)                                                                                                                tCGTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGgtg                               >  4:743379/1‑150 (MQ=255)                                                                                                                  gTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTc                             <  4:947412/150‑1 (MQ=255)                                                                                                                  gTTATTTAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCg                            >  4:428309/1‑151 (MQ=255)                                                                                                                       ttAGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACa                        <  2:1660142/150‑1 (MQ=255)                                                                                                                         aGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAAt                      <  3:413016/150‑1 (MQ=255)                                                                                                                         aGTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAAt                      <  4:1658/150‑1 (MQ=255)                                                                                                                          gTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAAt                      >  3:172195/1‑149 (MQ=255)                                                                                                                          gTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATc                     >  4:1508699/1‑150 (MQ=255)                                                                                                                          gTGCGGCAGGGCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCa                    <  4:230180/151‑1 (MQ=255)                                                                                                                            gCGGCAGGGCCGGGTAATCGATCTGCTGGGGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAgg                  >  2:336951/1‑151 (MQ=255)                                                                                                                                  gggCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAg            <  4:1119324/151‑1 (MQ=255)                                                                                                                                   ggCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAgg           <  2:701943/151‑1 (MQ=255)                                                                                                                                    gCCGGGTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGt          <  4:75271/151‑1 (MQ=255)                                                                                                                                         gTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTAt     <  4:702492/151‑1 (MQ=255)                                                                                                                                         gTAATCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTAt     <  3:80501/151‑1 (MQ=255)                                                                                                                                           aaTCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTAtt    <  4:1440285/150‑1 (MQ=255)                                                                                                                                            aTCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGt  <  2:1406529/151‑1 (MQ=255)                                                                                                                                            aTCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGt  >  2:1589620/1‑151 (MQ=255)                                                                                                                                             tCGATCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGt  <  2:736466/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TTGGTTACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGT  >  NC_000913/562849‑563138

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