Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,119	T→C	V131A (GTA→GCA)	fimI →	fimbrial protein involved in type 1 pilus biosynthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,119	0	T	C	100.0%	247.1 / NA	68	V131A (GTA→GCA)	fimI	fimbrial protein involved in type 1 pilus biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/32);  total (36/32)

ATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAG  >  NC_000913/4543973‑4544264                                                                                                                                                   |                                                                                                                                                  aTGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCg                                                                                                                                               <  2:896449/151‑1 (MQ=255) aTGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCg                                                                                                                                               >  2:630751/1‑151 (MQ=255)    tAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGAtt                                                                                                                                            >  1:342163/1‑151 (MQ=255)      gCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTa                                                                                                                                           >  2:939506/1‑150 (MQ=255)      gCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTa                                                                                                                                           >  2:323903/1‑150 (MQ=255)       cACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAAt                                                                                                                                         >  2:796500/1‑151 (MQ=255)        aCGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTTGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATc                                                                                                                                        <  1:675371/151‑1 (MQ=255)        aCGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATc                                                                                                                                        >  1:1135336/1‑151 (MQ=255)         cGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCg                                                                                                                                       >  2:428894/1‑151 (MQ=255)                gagtgaACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGtcctcc                                                                                                                                 >  2:560884/1‑150 (MQ=255)                gagtgaACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCa                                                                                                                                <  1:1075173/151‑1 (MQ=255)                    gaACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCa                                                                                                                             >  1:635493/1‑150 (MQ=255)                       cGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACt                                                                                                                         <  2:614421/151‑1 (MQ=255)                       cGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACt                                                                                                                         <  1:1048302/151‑1 (MQ=255)                        gtgtAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTg                                                                                                                        >  2:791686/1‑151 (MQ=255)                             ggagtgGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGaaa                                                                                                                    >  1:1171542/4‑150 (MQ=255)                               tgtgGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACgg                                                                                                                 >  1:855109/1‑151 (MQ=255)                                gtgGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACgg                                                                                                                 <  2:901622/150‑1 (MQ=255)                                gtgGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACgg                                                                                                                 <  1:1039710/150‑1 (MQ=255)                                      tttCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTAtt                                                                                                          >  1:832886/1‑151 (MQ=255)                                        tCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTc                                                                                                         >  1:400198/1‑150 (MQ=255)                                        tCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTc                                                                                                         >  1:702656/1‑150 (MQ=255)                                        tCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTc                                                                                                         >  1:8799/1‑150 (MQ=255)                                         cACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCa                                                                                                        <  2:1145573/150‑1 (MQ=255)                                          aCGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAgg                                                                                                      >  1:145979/1‑151 (MQ=255)                                               gtCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGctct                                                                                                  <  1:1114849/150‑1 (MQ=255)                                               gtCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTa                                                                                                 <  2:373167/151‑1 (MQ=255)                                                  gcgGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACtt                                                                                              >  2:711988/1‑151 (MQ=255)                                                    ggATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCg                                                                                            >  1:147921/1‑151 (MQ=255)                                                      aTGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCt                                                                                          >  1:477855/1‑151 (MQ=255)                                                       tGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCt                                                                                          >  1:719470/1‑150 (MQ=255)                                                          tAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACAt                                                                                      >  2:568940/1‑151 (MQ=255)                                                            aaaaTCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACAttt                                                                                    <  1:568936/151‑1 (MQ=255)                                                               aTCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCAt                                                                                 <  1:560881/151‑1 (MQ=255)                                                                 ccGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATc                                                                                <  1:101137/150‑1 (MQ=255)                                                                 ccGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATc                                                                                <  2:794156/150‑1 (MQ=255)                                                                 ccGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCATAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCATCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATc                                                                                <  1:726148/150‑1 (MQ=255)                                                                 ccGGATGTGCTTTCCGTGGGAGAGGGCCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATc                                                                                <  2:1020646/150‑1 (MQ=255)                                                                       gtgCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAat                                                                         >  2:529345/1‑151 (MQ=255)                                                                             tCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGt                                                                    <  1:958783/150‑1 (MQ=255)                                                                             tCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGt                                                                    <  1:1058034/150‑1 (MQ=255)                                                                              ccGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTg                                                                   <  2:228448/150‑1 (MQ=255)                                                                                gTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTa                                                                <  1:706084/151‑1 (MQ=255)                                                                                gTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTa                                                                <  2:702666/151‑1 (MQ=255)                                                                                      gagGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGc                                                          <  1:1065175/151‑1 (MQ=255)                                                                                        ggggCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGt                                                        <  1:939487/151‑1 (MQ=255)                                                                                              aGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGtt                                                  >  2:740591/1‑151 (MQ=255)                                                                                                ggATAGCCACCAATATTGGCGTAGCGTTGTTTGACTCTGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTTCTATCGGGCGTCGGGTTAc                                                >  1:781036/1‑151 (MQ=255)                                                                                                  aTAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTg                                              >  1:203380/1‑151 (MQ=255)                                                                                                          caATATCGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTAGTGGCGGCATc                                      >  2:103392/1‑151 (MQ=255)                                                                                                             tatTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGcc                                   >  1:400171/1‑151 (MQ=255)                                                                                                               ttGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCa                                  >  1:886764/1‑150 (MQ=255)                                                                                                                 ggCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATg                               >  1:298053/1‑151 (MQ=255)                                                                                                                  gCGTAGCGTTGTTTGATGATGAAGGAACCCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGcc                                   <  2:504099/146‑1 (MQ=255)                                                                                                                  gCGTAGCGTTGTTTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGc                              <  2:510210/151‑1 (MQ=255)                                                                                                                         gttgttTGATGATGAAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGcc                       <  2:477858/151‑1 (MQ=255)                                                                                                                               tgatgatgaAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTtc                 >  1:514217/1‑151 (MQ=255)                                                                                                                                gatgatgaAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTtct                >  1:410579/1‑151 (MQ=255)                                                                                                                                 atgatgaAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTtctc               <  2:125695/151‑1 (MQ=255)                                                                                                                                 atgatgaAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTtctc               >  1:396024/1‑151 (MQ=255)                                                                                                                                 atgatgaAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTtctc               <  1:347506/151‑1 (MQ=255)                                                                                                                                     tgaAGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCttt            >  2:697297/1‑150 (MQ=255)                                                                                                                                        aGGAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAAcc        <  2:400174/151‑1 (MQ=255)                                                                                                                                         ggAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCt       <  1:131835/151‑1 (MQ=255)                                                                                                                                          gAAACCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCt       >  2:698280/1‑150 (MQ=255)                                                                                                                                           aaaCCTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTAt     <  2:873580/151‑1 (MQ=255)                                                                                                                                              ccTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATc    >  2:154252/1‑149 (MQ=255)                                                                                                                                              ccTCGCACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAg  <  1:123289/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAG  >  NC_000913/4543973‑4544264

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