Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,408,612	C→T	S653S (AGC→AGT)	rnr →	exoribonuclease R, RNase R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,408,612	0	C	T	100.0%	328.3 / NA	85	S653S (AGC→AGT)	rnr	exoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (46/39);  total (46/39)

CGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGC  >  NC_000913/4408467‑4408756                                                                                                                                                  |                                                                                                                                                 cGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCa                                                                                                                                              <  1:325722/150‑1 (MQ=255)    tGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACtgg                                                                                                                                          >  1:891436/1‑151 (MQ=255)     ggAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACtggc                                                                                                                                         >  2:383664/1‑151 (MQ=255)          agaTGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACtggctttg                                                                                                                                     <  2:548177/150‑1 (MQ=255)              gctgcAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGcttc                                                                                                                                >  1:797412/1‑151 (MQ=255)               ctgcAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGcttct                                                                                                                               >  1:520016/1‑151 (MQ=255)                 gcAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCttcttc                                                                                                                             >  1:1156731/1‑151 (MQ=255)                  cAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTcg                                                                                                                            >  2:757804/1‑151 (MQ=255)                   aaCTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTcgt                                                                                                                           >  1:535649/1‑151 (MQ=255)                    aCTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTcgtc                                                                                                                          >  1:337046/1‑151 (MQ=255)                    aCTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTcgt                                                                                                                           >  2:918174/1‑150 (MQ=255)                        ggTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTcgtccgtc                                                                                                                      >  2:509093/1‑151 (MQ=255)                        ggTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTcgtccgtc                                                                                                                      <  1:416239/151‑1 (MQ=255)                           cAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTgg                                                                                                                   >  2:1129161/1‑151 (MQ=255)                           cAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTgg                                                                                                                   >  1:384428/1‑151 (MQ=255)                           cAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTgg                                                                                                                   >  1:408850/1‑151 (MQ=255)                             gCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGgac                                                                                                                 <  1:557261/151‑1 (MQ=255)                                    tCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTg                                                                                                           <  2:350360/150‑1 (MQ=255)                                       aTGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCa                                                                                                       <  2:115333/151‑1 (MQ=255)                                         ggCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCAtt                                                                                                     >  2:817574/1‑151 (MQ=255)                                         ggCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCAtt                                                                                                     >  2:21099/1‑151 (MQ=255)                                           cGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGa                                                                                                   <  1:734574/151‑1 (MQ=255)                                           cGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGa                                                                                                   <  2:805527/151‑1 (MQ=255)                                            ggAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGAt                                                                                                  <  2:38829/151‑1 (MQ=255)                                             gAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGAtg                                                                                                 <  2:588294/151‑1 (MQ=255)                                               aCGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGAtggt                                                                                               <  1:644157/151‑1 (MQ=255)                                                 gtcgtGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGAtggtct                                                                                             >  2:1154507/1‑151 (MQ=255)                                                    gtGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGAtggtctggt                                                                                          <  1:995000/151‑1 (MQ=255)                                                          acgaAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGt                                                                                    <  1:397235/151‑1 (MQ=255)                                                           cgaAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGtc                                                                                   <  2:755276/151‑1 (MQ=255)                                                            gaAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGtct                                                                                  <  1:26371/151‑1 (MQ=255)                                                            gaAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGtct                                                                                  <  1:746399/151‑1 (MQ=255)                                                             aaGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGtctc                                                                                 >  2:389268/1‑151 (MQ=255)                                                                cAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCtt                                                                               <  1:354908/150‑1 (MQ=255)                                                                cAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTc                                                                              <  2:999693/151‑1 (MQ=255)                                                                cAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTc                                                                              <  1:508963/151‑1 (MQ=255)                                                                 aaCGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCg                                                                             <  1:693649/151‑1 (MQ=255)                                                                 aaCGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCg                                                                             <  2:1014512/151‑1 (MQ=255)                                                                    gcgcgATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTg                                                                          <  2:891642/151‑1 (MQ=255)                                                                      gcgATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTgg                                                                         >  1:1178810/1‑150 (MQ=255)                                                                      gcgATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGa                                                                        <  2:1175052/151‑1 (MQ=255)                                                                      gcgATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGa                                                                        <  1:385367/151‑1 (MQ=255)                                                                        gATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGAc                                                                       >  2:460127/1‑150 (MQ=255)                                                                            tgGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATg                                                                   >  1:1051757/1‑150 (MQ=255)                                                                             ggCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGa                                                                  <  1:1078359/150‑1 (MQ=255)                                                                              gCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGact                                                                >  1:453996/1‑151 (MQ=255)                                                                               cTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGact                                                                >  1:1024714/1‑150 (MQ=255)                                                                               cTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGact                                                                >  2:287395/1‑150 (MQ=255)                                                                               cTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTATTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGact                                                                >  2:607858/1‑150 (MQ=255)                                                                                tGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGactac                                                              >  1:914691/1‑151 (MQ=255)                                                                                 gACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGactact                                                             <  1:757632/151‑1 (MQ=255)                                                                                 gACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGactac                                                              <  2:1122799/150‑1 (MQ=255)                                                                                     ggCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATc                                                          <  2:733118/150‑1 (MQ=255)                                                                                      gCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGc                                                        >  1:1167727/1‑151 (MQ=255)                                                                                       cTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCt                                                       >  1:291989/1‑151 (MQ=255)                                                                                          aaGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTg                                                    >  1:783587/1‑151 (MQ=255)                                                                                           aGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGa                                                   <  1:304946/151‑1 (MQ=255)                                                                                              gtgACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGAcc                                                 >  2:928388/1‑150 (MQ=255)                                                                                              gtgACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCa                                                <  1:83012/151‑1 (MQ=255)                                                                                                gACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAg                                               <  2:977853/150‑1 (MQ=255)                                                                                                 aCTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAgg                                              >  1:205343/1‑150 (MQ=255)                                                                                                  cTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGt                                             >  1:404548/1‑150 (MQ=255)                                                                                                    tCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAg                                           <  2:561250/150‑1 (MQ=255)                                                                                                     cATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCCCTTCGCTGGACAATGACTACTATCGCTTTGACTAGGTa                                            >  2:826698/1‑148 (MQ=255)                                                                                                      aTGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGc                                        >  2:912596/1‑151 (MQ=255)                                                                                                      aTGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGc                                        >  1:1186491/1‑151 (MQ=255)                                                                                                        gCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCa                                       >  1:491920/1‑150 (MQ=255)                                                                                                          tCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACg                                    >  2:847490/1‑151 (MQ=255)                                                                                                           cGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGc                                   >  1:964921/1‑151 (MQ=255)                                                                                                             aCCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCt                                 >  1:249189/1‑151 (MQ=255)                                                                                                                aggtaggtAACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGAt                              <  1:618016/151‑1 (MQ=255)                                                                                                                 ggtaggtaACGTCGTTAAAGGCGTAATTTCCAGTGTTACTGGCTGTGGCTTCTTCGTCCGTCTGGACGACTTTTTCATTGATGGTCTGGTCCATGTCTCTTCGATGGACAATGACTACTATCGCGTTGACCAGGTAGGGCAACGCATGAt                              <  1:821527/150‑1 (MQ=255)                                                                                                                      gtaACGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATggggg                         >  1:1030420/1‑150 (MQ=255)                                                                                                                        aaCGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAAt                      >  1:1076109/1‑151 (MQ=255)                                                                                                                        aaCGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAAt                      >  1:998696/1‑151 (MQ=255)                                                                                                                          cGTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATcc                    >  2:60679/1‑151 (MQ=255)                                                                                                                           gTCTTTAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCa                   <  2:1152243/151‑1 (MQ=255)                                                                                                                               ttAAAGGCGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCgg               >  2:680750/1‑151 (MQ=255)                                                                                                                                      cGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGCCt        >  1:792446/1‑151 (MQ=255)                                                                                                                                      cGTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACt        <  1:903226/151‑1 (MQ=255)                                                                                                                                       gTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACtt       >  2:185876/1‑151 (MQ=255)                                                                                                                                       gTAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACtt       <  2:445675/151‑1 (MQ=255)                                                                                                                                        tAATTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACtt       <  1:917962/150‑1 (MQ=255)                                                                                                                                         aaTTTCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTAt     <  1:21097/151‑1 (MQ=255)                                                                                                                                            ttCCAGTGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGc  >  1:1099807/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGC  >  NC_000913/4408467‑4408756

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