Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,636,877	T→C	V41A (GTC→GCC)	ynfO →	uncharacterized protein, Qin prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,636,877	0	T	C	100.0%	236.6 / NA	64	V41A (GTC→GCC)	ynfO	uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/35);  total (29/35)

AGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCC  >  NC_000913/1636730‑1637020                                                                                                                                                    |                                                                                                                                                aGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCAt                                                                                                                                              >  1:944918/1‑151 (MQ=14) aGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCAt                                                                                                                                              >  2:718732/1‑151 (MQ=14)   cGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTg                                                                                                                                            >  2:828660/1‑151 (MQ=14)     acacTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGaa                                                                                                                                          <  2:70546/151‑1 (MQ=14)     acacTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGaa                                                                                                                                          <  2:952907/151‑1 (MQ=14)       acTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGc                                                                                                                                        >  2:407955/1‑151 (MQ=14)             tACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTgccg                                                                                                                                  >  1:893375/1‑151 (MQ=14)                     aaCCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGgaga                                                                                                                           >  1:407028/1‑150 (MQ=21)                          cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACtgt                                                                                                                       <  1:480144/149‑1 (MQ=38)                          cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGc                                                                                                                     >  2:80456/1‑151 (MQ=39)                           aTGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACtgt                                                                                                                       <  2:79812/148‑1 (MQ=38)                             gTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTc                                                                                                                   <  2:633550/150‑1 (MQ=255)                               cAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAg                                                                                                                 <  1:374298/150‑1 (MQ=255)                               cAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGa                                                                                                                <  1:91931/151‑1 (MQ=255)                               cAACCAAGAACAGAACCCGCAGAACAACAACCCGCAACATCGGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAg                                                                                                                 <  2:514765/150‑1 (MQ=255)                                    aagaaCAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGaaaaaa                                                                                                           <  1:50416/151‑1 (MQ=255)                                     agaacagaacCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGaaaaa                                                                                                            >  1:258257/1‑149 (MQ=255)                                      gaacagaacCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAaga                                                                                                         <  1:299955/151‑1 (MQ=255)                                                  cAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGaa                                                                                             <  1:261178/151‑1 (MQ=255)                                                  cAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGa                                                                                              <  2:430548/150‑1 (MQ=255)                                                     aacaacaacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGca                                                                                          <  1:959249/151‑1 (MQ=255)                                                      acaacaacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGca                                                                                          >  1:705437/1‑150 (MQ=255)                                                        aacaacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaac                                                                                       >  2:154663/1‑151 (MQ=255)                                                         acaacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaaca                                                                                      <  2:764992/151‑1 (MQ=255)                                                         acaacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaac                                                                                       >  2:762076/1‑150 (MQ=255)                                                          caacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaacaa                                                                                     <  2:314036/151‑1 (MQ=255)                                                           aacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaacaaa                                                                                    <  1:3259/151‑1 (MQ=255)                                                           aacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaacaaa                                                                                    <  1:913168/151‑1 (MQ=255)                                                           aacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaacaaa                                                                                    <  2:944834/151‑1 (MQ=255)                                                            acCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaacaa                                                                                     >  1:18730/1‑149 (MQ=255)                                                            acCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaag                                                                                   <  2:123884/151‑1 (MQ=255)                                                             cccGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaacaa                                                                                     >  2:437182/1‑148 (MQ=255)                                                             cccGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaaga                                                                                  <  1:903943/151‑1 (MQ=255)                                                                 cAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaaa                                                                               <  2:292566/150‑1 (MQ=255)                                                                 cAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaa                                                                                >  1:727285/1‑149 (MQ=255)                                                                   aCATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAagag                                                                            <  1:504231/151‑1 (MQ=255)                                                                     acccGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTg                                                                          >  1:826629/3‑151 (MQ=255)                                                                     aTCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTg                                                                          <  2:330386/151‑1 (MQ=255)                                                                      tCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTg                                                                          >  1:818055/1‑150 (MQ=255)                                                                        cGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACa                                                                       <  2:179917/151‑1 (MQ=255)                                                                           tttCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACAtt                                                                     >  2:990659/1‑150 (MQ=255)                                                                            ttCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTa                                                                    <  1:632210/150‑1 (MQ=255)                                                                              ccTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTg                                                                 <  2:525680/151‑1 (MQ=255)                                                                                tAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGa                                                                >  1:535572/1‑150 (MQ=255)                                                                                 aaCCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGaa                                                               <  2:779334/150‑1 (MQ=255)                                                                                    cAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGc                                                           >  1:621243/1‑151 (MQ=255)                                                                                    cAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGc                                                           >  2:711051/1‑151 (MQ=255)                                                                                             gAACAAATTAACATCTCTCTTGAGCAAAAAGGGTCCGGGAATGTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGaaaa                                                   <  1:219770/150‑1 (MQ=255)                                                                                              aaCAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAAc                                                  <  1:298621/150‑1 (MQ=255)                                                                                               aCAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAAc                                                  >  2:560796/1‑149 (MQ=255)                                                                                                      aaCATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGAc                                         <  2:703644/151‑1 (MQ=255)                                                                                                       aCATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGAcc                                        >  1:612254/1‑151 (MQ=255)                                                                                                           cGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGa                                     <  2:494870/150‑1 (MQ=255)                                                                                                           cGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGAt                                    <  2:51490/151‑1 (MQ=255)                                                                                                               ctTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGa                                 >  2:514741/1‑150 (MQ=255)                                                                                                                ttGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGa                                 >  1:272061/1‑149 (MQ=255)                                                                                                                       aaaaGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAAttt                         >  2:314930/1‑150 (MQ=255)                                                                                                                                    aaTTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGc             >  1:980778/1‑149 (MQ=255)                                                                                                                                     aTTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGc             >  2:969265/1‑148 (MQ=255)                                                                                                                                      tttCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGAc         >  1:266629/1‑151 (MQ=255)                                                                                                                                        tctcAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACt        <  1:558225/150‑1 (MQ=255)                                                                                                                                          tcAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACg     >  1:639486/1‑151 (MQ=255)                                                                                                                                           cAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACg     <  2:449560/150‑1 (MQ=255)                                                                                                                                             gCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGccc  <  2:607088/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                AGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCC  >  NC_000913/1636730‑1637020

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