Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,042,089	T→C	S708G (AGC→GGC)	etk ←	tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,042,089	0	T	C	100.0%	376.3 / NA	98	S708G (AGC→GGC)	etk	tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (43/55);  total (43/55)

CTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGA  >  NC_000913/1041941‑1042235                                                                                                                                                     |                                                                                                                                                   cTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCgg                                                                                                                                                  >  2:213210/1‑151 (MQ=255)   gAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGgc                                                                                                                                                 <  4:138194/150‑1 (MQ=255)       gCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGtt                                                                                                                                            <  4:152769/151‑1 (MQ=255)       gCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGt                                                                                                                                             <  2:624806/150‑1 (MQ=255)          ctGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGttt                                                                                                                                           <  1:770622/149‑1 (MQ=255)          ctGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGttt                                                                                                                                           <  1:303778/149‑1 (MQ=255)          ctGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTa                                                                                                                                          >  2:338741/1‑150 (MQ=255)          ctGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTa                                                                                                                                          >  1:954571/1‑150 (MQ=255)            gAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTaa                                                                                                                                         <  1:957449/149‑1 (MQ=255)              ggCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATca                                                                                                                                      <  1:95429/150‑1 (MQ=255)                cTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCacac                                                                                                                                   >  1:182795/1‑151 (MQ=255)                   atgaCAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTACCTGTAAGCGGGGCCGGCGCGTTTAATCACACCa                                                                                                                                 >  3:3033/1‑150 (MQ=255)                     gaCAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCAtt                                                                                                                               <  1:251267/150‑1 (MQ=255)                     gaCAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCAtt                                                                                                                               <  3:102608/150‑1 (MQ=255)                     gaCAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTg                                                                                                                              <  2:703931/151‑1 (MQ=255)                     gaCAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTg                                                                                                                              <  2:713550/151‑1 (MQ=255)                      aCAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTg                                                                                                                              <  1:813500/150‑1 (MQ=255)                       cAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTg                                                                                                                              >  3:250902/1‑149 (MQ=255)                       cAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTg                                                                                                                              >  4:171850/1‑149 (MQ=255)                        aaaCGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAgg                                                                                                                           >  2:402987/1‑151 (MQ=255)                        aaaCGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAgg                                                                                                                           <  2:364039/151‑1 (MQ=255)                        aaaCGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAgg                                                                                                                           >  2:954324/1‑151 (MQ=255)                         aaCGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAgg                                                                                                                           >  2:318894/1‑150 (MQ=255)                           cGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGAt                                                                                                                         <  3:351453/150‑1 (MQ=255)                            gaagaaCTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAg                                                                                                                       <  2:648613/151‑1 (MQ=255)                               gaaCTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAgcgc                                                                                                                    >  2:850394/1‑151 (MQ=255)                                 aCTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCt                                                                                                                  >  2:1254815/1‑151 (MQ=255)                                 aCTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCt                                                                                                                  >  4:439659/1‑151 (MQ=255)                                    gTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCttt                                                                                                                <  3:218407/150‑1 (MQ=255)                                    gTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTa                                                                                                               <  4:387537/151‑1 (MQ=255)                                       tAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAat                                                                                                             >  4:171006/1‑150 (MQ=255)                                        aaTGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAatat                                                                                                           <  1:634703/151‑1 (MQ=255)                                        aaTGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAatat                                                                                                           >  1:1154458/1‑151 (MQ=255)                                         aTGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATAtt                                                                                                          >  3:5564/1‑151 (MQ=255)                                         aTGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATAtt                                                                                                          >  1:136203/1‑151 (MQ=255)                                           gCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTg                                                                                                         <  2:844253/150‑1 (MQ=255)                                           gCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGa                                                                                                        <  4:113562/151‑1 (MQ=255)                                            cGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGAc                                                                                                       <  1:905260/151‑1 (MQ=255)                                                gACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgcct                                                                                                   <  2:653387/151‑1 (MQ=255)                                                      aaaaGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGt                                                                                             >  1:680239/1‑151 (MQ=255)                                                           gCGTCCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTCATATTGACGCCTGCCTGTTcc                                                                                          >  1:327796/1‑149 (MQ=255)                                                            cGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCTTTTAATCACCCCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                                                       <  2:1217890/151‑1 (MQ=255)                                                             gTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                                                       <  1:1125361/150‑1 (MQ=255)                                                                cACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGc                                                                                    <  4:198949/150‑1 (MQ=255)                                                                cACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCt                                                                                   <  1:1162311/151‑1 (MQ=255)                                                                 aCGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTg                                                                                  <  1:981426/151‑1 (MQ=255)                                                                  cgccgcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGc                                                                                 <  1:893092/151‑1 (MQ=255)                                                                  cgccgcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGc                                                                                 <  1:371915/151‑1 (MQ=255)                                                                     cgcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAtt                                                                              <  4:304450/151‑1 (MQ=255)                                                                      gcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTg                                                                             <  1:302202/151‑1 (MQ=255)                                                                        aTCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACTCCTGCCTGTTCCAGACGCTGCATTGAc                                                                           <  1:1247090/151‑1 (MQ=255)                                                                         tCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGAc                                                                           >  1:396446/1‑150 (MQ=255)                                                                          ccGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACa                                                                          <  1:916255/150‑1 (MQ=255)                                                                            ggCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAAc                                                                       >  4:153698/1‑151 (MQ=255)                                                                              ccACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTc                                                                     >  1:438163/1‑151 (MQ=255)                                                                              ccACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTc                                                                     >  1:833180/1‑151 (MQ=255)                                                                              ccACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTc                                                                     >  2:1245288/1‑151 (MQ=255)                                                                              ccACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTc                                                                     >  2:1204895/1‑151 (MQ=255)                                                                                aCTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCAc                                                                   <  2:713145/151‑1 (MQ=255)                                                                                  tttCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACct                                                                 >  1:1040132/1‑151 (MQ=255)                                                                                   ttCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACctc                                                                >  1:323992/1‑151 (MQ=255)                                                                                     cAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCtt                                                              <  1:633807/151‑1 (MQ=255)                                                                                      aGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCtt                                                              <  1:831798/150‑1 (MQ=255)                                                                                       gTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTg                                                            <  1:170801/151‑1 (MQ=255)                                                                                         tttACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTgg                                                           <  3:27432/150‑1 (MQ=255)                                                                                          ttACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGc                                                          <  4:307687/150‑1 (MQ=255)                                                                                           tACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCg                                                         <  1:850566/150‑1 (MQ=255)                                                                                              tctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtgt                                                     >  4:451364/1‑151 (MQ=255)                                                                                                 ttCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTc                                                   >  3:85726/1‑150 (MQ=255)                                                                                                  tctcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCa                                                  >  2:1049670/1‑150 (MQ=255)                                                                                                       gAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaa                                            <  1:1133169/151‑1 (MQ=255)                                                                                                         gTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaaa                                          <  2:367146/151‑1 (MQ=255)                                                                                                         gTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaa                                           <  2:1261687/150‑1 (MQ=255)                                                                                                         gTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaa                                           >  1:477779/1‑150 (MQ=255)                                                                                                           aaCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcg                                        >  4:344587/1‑151 (MQ=255)                                                                                                            aCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc                                       <  2:1271926/151‑1 (MQ=255)                                                                                                            aCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc                                       >  1:702940/1‑151 (MQ=255)                                                                                                                 aaCCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAAc                                  >  3:334138/1‑151 (MQ=255)                                                                                                                 aaCCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAAc                                  >  4:402083/1‑151 (MQ=255)                                                                                                                     gTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagc                              <  1:340840/151‑1 (MQ=255)                                                                                                                     gTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagc                              >  1:373739/1‑151 (MQ=255)                                                                                                                     gTAATAGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagc                              >  1:543952/1‑151 (MQ=255)                                                                                                                          aGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGc                          <  1:1193386/150‑1 (MQ=255)                                                                                                                          aGTTATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCt                         <  2:182771/151‑1 (MQ=255)                                                                                                                            ttATAGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTgg                       >  3:393438/1‑151 (MQ=255)                                                                                                                              ataGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTg                     <  2:954381/151‑1 (MQ=255)                                                                                                                               taGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGc                    >  4:342552/1‑151 (MQ=255)                                                                                                                                aGCCGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGcc                   >  4:105129/1‑151 (MQ=255)                                                                                                                                   cGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAAc                >  1:998018/1‑151 (MQ=255)                                                                                                                                   cGTAACTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAAc                <  1:1165210/151‑1 (MQ=255)                                                                                                                                      aaCTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAg              <  2:1034072/150‑1 (MQ=255)                                                                                                                                      aaCTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAg              <  2:1187215/150‑1 (MQ=255)                                                                                                                                       aCTGTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAAAAGa             >  4:398575/1‑150 (MQ=255)                                                                                                                                          gTAAGCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACg          <  1:1068207/150‑1 (MQ=255)                                                                                                                                              gCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCa     <  3:388214/151‑1 (MQ=255)                                                                                                                                              gCGGTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCa     <  3:296988/151‑1 (MQ=255)                                                                                                                                                ggTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACg   >  1:92743/1‑151 (MQ=255)                                                                                                                                                 gTGCCGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGa  <  3:345446/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGA  >  NC_000913/1041941‑1042235

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