Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,430,519	T→C	S491P (TCA→CCA)	stfR →	Rac prophage; putative tail fiber protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,430,519	0	T	C	100.0%	276.2 / NA	76	S491P (TCA→CCA)	stfR	Rac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/47);  total (29/47)

CCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTAACG  >  NC_000913/1430373‑1430666                                                                                                                                                   |                                                                                                                                                    ccGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGc                                                                                                                                                 <  2:148527/151‑1 (MQ=255)   gagaCAGCGGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCt                                                                                                                                               <  1:436802/151‑1 (MQ=255)   gagaCAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCt                                                                                                                                               <  1:233545/151‑1 (MQ=255)     gaCAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTAt                                                                                                                                             <  2:473473/151‑1 (MQ=255)     gaCAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTAt                                                                                                                                             >  1:832609/1‑151 (MQ=255)      aCAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATg                                                                                                                                            >  2:429749/1‑151 (MQ=255)       cagcagCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGa                                                                                                                                           <  1:57256/151‑1 (MQ=255)         gcagcTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACa                                                                                                                                         <  2:120376/151‑1 (MQ=255)          cagcTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACa                                                                                                                                         >  1:543036/1‑150 (MQ=255)          cagcTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACa                                                                                                                                         >  1:336173/1‑150 (MQ=255)           agcTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAAt                                                                                                                                       <  2:858462/151‑1 (MQ=255)              tAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAAt                                                                                                                                       <  1:244942/148‑1 (MQ=255)                   gggCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGaa                                                                                                                               <  2:268420/151‑1 (MQ=255)                   gggCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGaa                                                                                                                               <  1:479476/151‑1 (MQ=255)                         agGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCt                                                                                                                         <  2:109179/151‑1 (MQ=255)                            atatTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGc                                                                                                                       >  2:412772/1‑150 (MQ=255)                            atatTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGc                                                                                                                       >  2:747489/1‑150 (MQ=255)                                tGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                                                                                   >  1:409078/1‑150 (MQ=255)                                 gCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAg                                                                                                                 <  2:711400/151‑1 (MQ=255)                                   aTCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCATTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAc                                                                                                               >  2:739231/1‑151 (MQ=255)                                   aTCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAc                                                                                                               >  2:641354/1‑151 (MQ=255)                                     ccgccgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAcc                                                                                                              <  2:340141/150‑1 (MQ=255)                                      cgccgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAcc                                                                                                              >  2:310264/1‑149 (MQ=255)                                        ccgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGaa                                                                                                          <  1:298657/151‑1 (MQ=255)                                        ccgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGa                                                                                                           >  2:609684/1‑150 (MQ=255)                                         cgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAAc                                                                                                         <  2:415148/151‑1 (MQ=255)                                            ggCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGgc                                                                                                      <  2:212454/151‑1 (MQ=255)                                             gcgcTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGgcg                                                                                                     <  1:539454/151‑1 (MQ=255)                                               gcTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCt                                                                                                   >  1:521566/1‑151 (MQ=255)                                               gcTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCt                                                                                                   <  2:74854/151‑1 (MQ=255)                                                cTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTg                                                                                                  >  1:719877/1‑151 (MQ=255)                                                   gAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGat                                                                                                <  1:429717/150‑1 (MQ=255)                                                     ggATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGAtata                                                                                             <  2:445550/151‑1 (MQ=255)                                                       aTGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATAcc                                                                                           >  2:848132/1‑151 (MQ=255)                                                         gCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATAccc                                                                                          <  2:660243/150‑1 (MQ=255)                                                               acgacgAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAg                                                                                    <  1:747639/150‑1 (MQ=255)                                                                 gacgaAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAggg                                                                                  >  1:189993/1‑150 (MQ=255)                                                                  acgaAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGAt                                                                                >  2:495049/1‑151 (MQ=255)                                                                   cgaAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATg                                                                               >  2:237306/1‑151 (MQ=255)                                                                    gaAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGc                                                                              >  1:60532/1‑151 (MQ=255)                                                                      aaaaGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCtt                                                                            <  1:237282/151‑1 (MQ=255)                                                                       aaaGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCt                                                                             >  2:356150/1‑149 (MQ=255)                                                                        aaGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGcgtt                                                                                                                                              <  1:847105/83‑4 (MQ=255)                                                                        aaGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTcc                                                                          >  1:521356/1‑151 (MQ=255)                                                                         aGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCt                                                                         >  2:656606/1‑151 (MQ=255)                                                                                 gTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGaacaac                                                                  >  2:272174/1‑150 (MQ=255)                                                                                   aCAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACAtt                                                               <  2:546009/151‑1 (MQ=255)                                                                                    cAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACAtt                                                               >  2:241577/1‑150 (MQ=255)                                                                                         cagcagTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcgcg                                                         <  2:218489/151‑1 (MQ=255)                                                                                         cagcagTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcgcg                                                         >  2:463333/1‑151 (MQ=255)                                                                                         cagcagTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcgc                                                          <  1:24103/150‑1 (MQ=255)                                                                                         cagcagTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcg                                                           <  1:105950/149‑1 (MQ=255)                                                                                            cagTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTc                                                      <  2:205583/151‑1 (MQ=255)                                                                                             agTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCa                                                     <  1:556284/151‑1 (MQ=255)                                                                                                 cGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTaa                                                 <  2:526605/151‑1 (MQ=255)                                                                                                  gACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTaaa                                                <  1:546548/151‑1 (MQ=255)                                                                                                       aCAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGa                                           <  2:196726/151‑1 (MQ=255)                                                                                                         aGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACt                                         <  1:102819/151‑1 (MQ=255)                                                                                                            aCGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACttt                                       <  1:475255/150‑1 (MQ=255)                                                                                                            aCGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTg                                      <  2:678105/151‑1 (MQ=255)                                                                                                             cGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGc                                     <  1:471235/151‑1 (MQ=255)                                                                                                             cGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGc                                     <  1:527248/151‑1 (MQ=255)                                                                                                              gTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCt                                    <  2:683962/151‑1 (MQ=255)                                                                                                              gTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCt                                    >  1:826664/1‑151 (MQ=255)                                                                                                                cTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTg                                   >  2:604670/1‑150 (MQ=255)                                                                                                                   aaaCGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATaa                               <  2:675501/151‑1 (MQ=255)                                                                                                                      cGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCg                            >  1:354802/1‑151 (MQ=255)                                                                                                                           gcggcAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTa                       <  1:364399/151‑1 (MQ=255)                                                                                                                                     cAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTAtgt             <  2:630562/151‑1 (MQ=255)                                                                                                                                         ggCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCgg         >  2:118048/1‑151 (MQ=255)                                                                                                                                         ggCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCgg         <  2:629081/151‑1 (MQ=255)                                                                                                                                          gCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCggg        <  1:444001/151‑1 (MQ=255)                                                                                                                                           cAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGt       <  1:15923/151‑1 (MQ=255)                                                                                                                                             gTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAGGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTa     <  2:101239/151‑1 (MQ=255)                                                                                                                                              tAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTa     >  2:333118/1‑150 (MQ=255)                                                                                                                                                aaaCCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTAACg  <  1:184384/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTAACG  >  NC_000913/1430373‑1430666

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