Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,668,330	G→A	P287S (CCG→TCG)	yhjA ←	putative cytochrome C peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,668,330	0	G	A	100.0%	280.7 / NA	71	P287S (CCG→TCG)	yhjA	putative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/37);  total (34/37)

TCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGAGTTAAATACCGT  >  NC_000913/3668183‑3668477                                                                                                                                                    |                                                                                                                                                    tCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATc                                                                                                                                                  >  2:228703/1‑151 (MQ=255)  cAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCt                                                                                                                                                 <  2:201793/151‑1 (MQ=255)     cgcATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCtttt                                                                                                                                              >  1:257438/1‑151 (MQ=255)      gcATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCttttt                                                                                                                                             >  2:360894/1‑151 (MQ=255)        aTTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTcc                                                                                                                                           >  1:879886/1‑151 (MQ=255)             cATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCtt                                                                                                                                      >  1:111503/1‑151 (MQ=255)             cATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCt                                                                                                                                       <  1:717928/150‑1 (MQ=255)             cATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCt                                                                                                                                       <  1:1052761/150‑1 (MQ=255)                  ccGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGc                                                                                                                                   <  2:629845/149‑1 (MQ=255)                  ccGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCa                                                                                                                                  <  2:30949/150‑1 (MQ=255)                  ccGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCa                                                                                                                                  >  1:404485/1‑150 (MQ=255)                   cgcgCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCaat                                                                                                                                <  2:185429/151‑1 (MQ=255)                        aaCCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAAttt                                                                                                                           <  2:775620/151‑1 (MQ=255)                         aCCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTc                                                                                                                          <  1:382860/151‑1 (MQ=255)                           cATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGt                                                                                                                        <  1:200701/151‑1 (MQ=255)                            aTTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTc                                                                                                                       >  2:737078/1‑151 (MQ=255)                            aTTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTc                                                                                                                       <  1:924116/151‑1 (MQ=255)                              ttATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTcc                                                                                                                      <  1:798956/150‑1 (MQ=255)                                   aaaTGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGAt                                                                                                                <  1:985032/151‑1 (MQ=255)                                     aTGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGT‑CCAGGATTTc                                                                                                             <  2:608637/151‑1 (MQ=255)                                         ggAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGAc                                                                                                          <  1:95337/151‑1 (MQ=255)                                            aTCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGcc                                                                                                       >  1:201296/1‑151 (MQ=255)                                            aTCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGcc                                                                                                       >  2:942293/1‑151 (MQ=255)                                                 gCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTc                                                                                                  <  2:641949/151‑1 (MQ=255)                                                  cGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCa                                                                                                 <  1:1026267/151‑1 (MQ=255)                                                     aattaatGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATc                                                                                              >  1:542034/1‑151 (MQ=255)                                                        taatGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCgg                                                                                            >  2:327412/1‑150 (MQ=255)                                                           tGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTc                                                                                        >  1:522446/1‑151 (MQ=255)                                                            gTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTc                                                                                        >  1:978988/1‑150 (MQ=255)                                                            gTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTc                                                                                        >  1:432660/1‑150 (MQ=255)                                                                     aaTTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTc                                                                              >  1:47848/1‑151 (MQ=255)                                                                       ttCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTcc                                                                             >  1:1013726/1‑150 (MQ=255)                                                                             aaTGGCATCAGTAATATTTTCGCCACTGAAACCTCGCGGATAGACTTCGAGGATCTGCGTTTTAAGCTGCGAATATTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATGGGGTTCAACGGCGGTcccc                                                                           >  2:172639/1‑144 (MQ=255)                                                                              aTGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCt                                                                     <  2:589746/151‑1 (MQ=255)                                                                                tgCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTg                                                                    >  1:934159/2‑150 (MQ=255)                                                                                   aTCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATcc                                                                >  1:833409/1‑151 (MQ=255)                                                                                      aGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGc                                                             >  1:1042256/1‑151 (MQ=255)                                                                                         aaTATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGCTTCAGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAAt                                                          >  2:830143/1‑151 (MQ=255)                                                                                         aaTATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCaa                                                           <  2:919222/150‑1 (MQ=255)                                                                                           tatTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATg                                                         >  1:740615/1‑150 (MQ=255)                                                                                            atTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGtt                                                       <  1:1066074/151‑1 (MQ=255)                                                                                            atTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGAAATGt                                                        <  2:493746/150‑1 (MQ=255)                                                                                                tCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCACCTTAGCAATAATTTCGTCCCAGGTTTTCGCCGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGc                                                     >  1:20143/1‑149 (MQ=255)                                                                                                 cGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGcca                                                   <  2:393488/150‑2 (MQ=255)                                                                                                  gCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCa                                                 <  2:833363/151‑1 (MQ=255)                                                                                                  gCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCa                                                 <  1:180483/151‑1 (MQ=255)                                                                                                   ccACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCCAATCCTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCTGTCCTCCAGCCTGATCCTGCAATGTTGcc                                                    <  2:4417/147‑1 (MQ=255)                                                                                                        gAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGAcc                                            <  1:234455/150‑1 (MQ=255)                                                                                                         aaaCCTTGCGGATAGACTTCGGGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAAGGATGCCGCACGACCa                                           >  1:130576/1‑150 (MQ=255)                                                                                                         aaaCCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCAt                                          >  1:426801/1‑151 (MQ=255)                                                                                                               tGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCa                                      <  1:396280/149‑1 (MQ=255)                                                                                                                gCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGa                                    <  2:437492/150‑1 (MQ=255)                                                                                                                 cGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAAc                                  >  1:909104/1‑151 (MQ=255)                                                                                                                     tAGACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCt                              >  1:209714/1‑151 (MQ=255)                                                                                                                       gACTTCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCa                            <  1:543294/151‑1 (MQ=255)                                                                                                                          ttCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACg                         >  1:375038/1‑151 (MQ=255)                                                                                                                           tCGAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGt                        >  2:383916/1‑151 (MQ=255)                                                                                                                             gAGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTa                      <  2:642336/151‑1 (MQ=255)                                                                                                                              aGGAACTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTaa                     <  2:809478/151‑1 (MQ=255)                                                                                                                                 aaCTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAAt                   >  1:127778/1‑150 (MQ=255)                                                                                                                                  aCTGCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATAc                 >  1:233710/1‑151 (MQ=255)                                                                                                                                     gCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTg               >  1:902212/1‑150 (MQ=255)                                                                                                                                     gCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGa              <  2:391320/151‑1 (MQ=255)                                                                                                                                     gCGTTTTAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGa              <  1:184656/151‑1 (MQ=255)                                                                                                                                         tttAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGAGtt           <  2:450055/150‑1 (MQ=255)                                                                                                                                           tAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGtt                       >  2:683252/1‑136 (MQ=255)                                                                                                                                           tAAGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGtt                       <  1:683286/136‑1 (MQ=255)                                                                                                                                            aaGCTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGAGTTaaa        <  1:628622/150‑1 (MQ=255)                                                                                                                                               cTGCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGAGTTAAATAcc    >  2:643211/1‑151 (MQ=255)                                                                                                                                                 gCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGAGTTAAATAcc    <  2:233685/149‑1 (MQ=255)                                                                                                                                                 gCGAATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGAGTTAAATACCGt  <  2:156438/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAACGGCGGTCCACCAGCCTGATCCTGCAATGTTGCCGCACGACCATCCCAGAACTGCTCAACGTTAAATACTGAGTTAAATACCGT  >  NC_000913/3668183‑3668477

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