Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,331,551	T→C	pseudogene (2852/4605 nt)	yfaS ←	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,331,551	0	T	C	100.0%	284.5 / NA	76	pseudogene (2852/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (46/30);  total (46/30)

CCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCG  >  NC_000913/2331403‑2331697                                                                                                                                                     |                                                                                                                                                   cccATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCccgc                                                                                                                                                  <  1:834712/151‑1 (MQ=255)    aTTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGctga                                                                                                                                               <  2:668243/151‑1 (MQ=255)        tgttgCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGctgactg                                                                                                                                            <  2:992975/150‑1 (MQ=255)        tgttgCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGctgactg                                                                                                                                            <  2:1088218/150‑1 (MQ=255)          ttgCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGctgactga                                                                                                                                           >  2:738613/1‑149 (MQ=255)          ttgCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGctgactga                                                                                                                                           >  1:539905/1‑149 (MQ=255)              aGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATaa                                                                                                                                     <  2:804593/151‑1 (MQ=255)               gTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATaa                                                                                                                                     >  2:89959/1‑150 (MQ=255)                ttCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAAc                                                                                                                                   >  1:953508/1‑151 (MQ=255)                 tCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAAc                                                                                                                                   <  2:506722/150‑1 (MQ=255)                    gCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTc                                                                                                                               >  1:1149090/1‑151 (MQ=255)                     cACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTcc                                                                                                                              >  1:105186/1‑151 (MQ=255)                         tAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTgc                                                                                                                           <  2:631467/150‑1 (MQ=255)                                ccAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGtt                                                                                                                   <  2:880599/151‑1 (MQ=255)                                ccAGATTGCTGAATATTCCGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGtt                                                                                                                   >  2:1030396/1‑151 (MQ=255)                                  aGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTg                                                                                                                 >  1:274644/1‑151 (MQ=255)                                   gATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGt                                                                                                                >  2:373924/1‑151 (MQ=255)                                   gATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGt                                                                                                                <  2:1204967/151‑1 (MQ=255)                                   gATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGt                                                                                                                >  1:831793/1‑151 (MQ=255)                                    aTTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTc                                                                                                               >  1:1197766/1‑151 (MQ=255)                                      tccTGAATATTCTGCGTCAGCGAAATATAATACGCCCCTGTGTGCAGCGTCGGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGATGATAAACAGTCCTGCCGCCGGTTTg                                                                                                                 <  1:260513/145‑1 (MQ=255)                                         tGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGccc                                                                                                           >  1:225140/1‑150 (MQ=255)                                          gAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCAc                                                                                                         >  2:521317/1‑151 (MQ=255)                                          gAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCAc                                                                                                         <  1:760953/151‑1 (MQ=255)                                            atatTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCAcgc                                                                                                       >  2:19786/1‑151 (MQ=255)                                            atatTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCAcgc                                                                                                       >  2:1033028/1‑151 (MQ=255)                                             tatTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCAcgcg                                                                                                      >  1:1225848/1‑151 (MQ=255)                                             tatTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCAcgcg                                                                                                      >  1:946578/1‑151 (MQ=255)                                             tatTCTGCCTCAGCGAAATATAATTCGCCCATTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCCTCGGTTCGTTATCCCGCTGACTGAAGATAATCAGTCCTGCCGCCGGTTTgccccacac                                                                                                         >  1:802532/1‑140 (MQ=255)                                                    cGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACAcc                                                                                               <  2:113654/151‑1 (MQ=255)                                                       cAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGtt                                                                                            >  2:430989/1‑151 (MQ=255)                                                        aGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTg                                                                                           <  1:993017/151‑1 (MQ=255)                                                          cGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGc                                                                                         >  2:199338/1‑151 (MQ=255)                                                          cGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGCCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGc                                                                                         >  2:460291/1‑151 (MQ=255)                                                           gAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGc                                                                                         >  2:31276/1‑150 (MQ=255)                                                           gAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGc                                                                                         >  2:1209791/1‑150 (MQ=255)                                                           gAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGc                                                                                         >  2:1149804/1‑150 (MQ=255)                                                            aaaTATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCAt                                                                                       >  1:427969/1‑151 (MQ=255)                                                            aaaTATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCAt                                                                                       <  1:1204040/151‑1 (MQ=255)                                                              atatAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATAc                                                                                     >  2:549011/1‑151 (MQ=255)                                                                 taATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTcc                                                                                  >  1:492227/1‑151 (MQ=255)                                                                  aaTTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTcc                                                                                  >  2:786961/1‑150 (MQ=255)                                                                    ttCGCCCCTTTGTGCAGCGTCAGCGTCTGCCGCATCTCAGCGCCTGCCAATTTAGTCACCAGCGCCACCGGCTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACt                                                                               >  1:657465/1‑151 (MQ=255)                                                                      cGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTc                                                                             >  2:286434/1‑151 (MQ=255)                                                                       gCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTc                                                                             >  2:264082/1‑150 (MQ=255)                                                                        ccccTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCCCCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCa                                                                            <  1:16266/150‑1 (MQ=255)                                                                            tttGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGt                                                                         >  2:1193230/1‑149 (MQ=255)                                                                             ttGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTag                                                                       >  2:821247/1‑150 (MQ=255)                                                                             ttGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTag                                                                       <  1:734312/150‑1 (MQ=255)                                                                             ttGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTag                                                                       <  2:231405/150‑1 (MQ=255)                                                                             ttGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTag                                                                       <  2:977610/150‑1 (MQ=255)                                                                              tgtgCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTag                                                                       >  2:139526/1‑149 (MQ=255)                                                                               gtgCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAgaga                                                                    <  2:946960/151‑1 (MQ=255)                                                                                  cagcgtcagcgtCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGAttt                                                                 >  1:1214044/1‑151 (MQ=255)                                                                                  cagcgtcagcgtCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTGCATGTAGAGAttt                                                                 <  2:640998/151‑1 (MQ=255)                                                                                    gcgtcagcgtcTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGAttttt                                                               <  1:1174273/151‑1 (MQ=255)                                                                                            gtcTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAAc                                                       >  2:266541/1‑151 (MQ=255)                                                                                            gtcTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAAc                                                       >  1:354786/1‑151 (MQ=255)                                                                                               tGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGc                                                     >  2:441529/1‑150 (MQ=255)                                                                                                ggCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGc                                                     <  2:1159771/149‑1 (MQ=255)                                                                                                ggCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCa                                                    <  2:185381/150‑1 (MQ=255)                                                                                                  cgcATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGAt                                                 <  1:199335/151‑1 (MQ=255)                                                                                                              ccTGCAAATTTAGTCACCACCGCCACCGGTTCGTTACCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCt                                     <  1:1091823/151‑1 (MQ=255)                                                                                                                     aTTTAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGAc                              <  1:1207004/151‑1 (MQ=255)                                                                                                                       ttAGTCCCCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGAcc                             <  2:427999/150‑1 (MQ=255)                                                                                                                        tAGTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGCTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGAcc                             >  1:595468/1‑149 (MQ=255)                                                                                                                          gTCACCAGCGCCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGcc                         <  2:105181/151‑1 (MQ=255)                                                                                                                                 gcgcCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGcc                  <  2:1148688/151‑1 (MQ=255)                                                                                                                                  cgcCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCg                 >  1:1148301/1‑151 (MQ=255)                                                                                                                                  cgcCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCg                 >  1:876701/1‑151 (MQ=255)                                                                                                                                   gcCACCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGt                >  1:803994/1‑151 (MQ=255)                                                                                                                                      aCCGGTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTc              >  2:182522/1‑150 (MQ=255)                                                                                                                                          gTTCGTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATccc         >  1:1172476/1‑151 (MQ=255)                                                                                                                                              gTTATCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgc     <  1:846296/151‑1 (MQ=255)                                                                                                                                                 aTCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgcgcg  >  2:1241979/1‑151 (MQ=255)                                                                                                                                                 aTCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgcgcg  >  1:613254/1‑151 (MQ=255)                                                                                                                                                 aTCCCGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgcgc   <  2:125699/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCG  >  NC_000913/2331403‑2331697

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