Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,237,614	A→G	T240A (ACG→GCG)	yjbG →	extracellular polysaccharide export OMA protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,237,614	0	A	G	100.0%	241.1 / NA	64	T240A (ACG→GCG)	yjbG	extracellular polysaccharide export OMA protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/38);  total (26/38)

ACGGACGCACGCAAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACG  >  NC_000913/4237468‑4237759                                                                                                                                                   |                                                                                                                                                  acggacgCACGCAAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCa                                                                                                                                               <  1:522704/151‑1 (MQ=255)  cggacgCACGCAAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGgcgc                                                                                                                                                >  2:530615/1‑149 (MQ=255)    gacgCACGCAAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCg                                                                                                                                            <  2:49648/151‑1 (MQ=255)       gcacgcaAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGAt                                                                                                                                         <  2:311263/151‑1 (MQ=255)          cgcaAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATAcc                                                                                                                                      <  1:371513/151‑1 (MQ=255)           gcaAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCt                                                                                                                                     <  2:78393/151‑1 (MQ=255)            caAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTc                                                                                                                                    <  2:215922/151‑1 (MQ=255)            caAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTc                                                                                                                                    <  2:55486/151‑1 (MQ=255)               aaaGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCaat                                                                                                                                 >  2:350613/1‑151 (MQ=255)                 aGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCaataa                                                                                                                               <  2:162296/151‑1 (MQ=255)                   cgcCGGTGGCTTACTGGACCAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATg                                                                                                                             >  2:156464/1‑151 (MQ=255)                   cgcCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATg                                                                                                                             <  1:309660/151‑1 (MQ=255)                        gTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGaaaaa                                                                                                                        <  1:323763/151‑1 (MQ=255)                          ggCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGaaaaaa                                                                                                                       >  1:431419/1‑150 (MQ=255)                           gCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGaaaaaa                                                                                                                       >  1:353694/1‑149 (MQ=255)                             ttACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGAc                                                                                                                   >  2:511413/1‑151 (MQ=255)                              tACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGa                                                                                                                    <  2:35964/149‑1 (MQ=255)                              tACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGa                                                                                                                    <  2:488243/149‑1 (MQ=255)                              tACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGAc                                                                                                                   <  2:35631/150‑1 (MQ=255)                               aCTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACAt                                                                                                                 <  1:23267/151‑1 (MQ=255)                               aCTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACAt                                                                                                                 <  2:202982/151‑1 (MQ=255)                                  ggAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCt                                                                                                               <  2:335450/150‑1 (MQ=255)                                  ggAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCt                                                                                                               >  2:371259/1‑150 (MQ=255)                                  ggAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCt                                                                                                               <  1:507798/150‑1 (MQ=255)                                  ggAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTg                                                                                                              <  2:551349/151‑1 (MQ=255)                                       aaGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGACGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTgctt                                                                                                           >  1:65235/1‑149 (MQ=255)                                             cACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACt                                                                                                   <  1:456755/151‑1 (MQ=255)                                              aCGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTg                                                                                                  >  2:265031/1‑151 (MQ=255)                                              aCGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTg                                                                                                  <  1:604010/151‑1 (MQ=255)                                                  agagCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGgcgc                                                                                              >  2:629480/1‑151 (MQ=255)                                                  agagCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGgcg                                                                                               <  2:626488/150‑1 (MQ=255)                                                      ccGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTggg                                                                                          >  2:33046/1‑151 (MQ=255)                                                            cccGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAg                                                                                    <  2:430796/151‑1 (MQ=255)                                                               ggCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGcac                                                                                 <  1:541060/151‑1 (MQ=255)                                                                 cagcaTTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGcac                                                                                 >  2:383913/1‑149 (MQ=255)                                                                     attattTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTg                                                                           <  2:83679/151‑1 (MQ=255)                                                                         tttatGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTAc                                                                       >  2:166539/1‑151 (MQ=255)                                                                         tttatGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTAc                                                                       >  2:563325/1‑151 (MQ=255)                                                                            atGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACgg                                                                     <  1:154789/150‑1 (MQ=255)                                                                             tGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCg                                                                   >  1:479607/1‑151 (MQ=255)                                                                                 ggCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGaaa                                                                <  1:55266/150‑1 (MQ=255)                                                                                       gcgGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATcc                                                         >  2:294100/1‑151 (MQ=255)                                                                                                   tGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTgg                                             >  2:654743/1‑151 (MQ=255)                                                                                                    ggAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTgg                                             <  2:352398/150‑1 (MQ=255)                                                                                                    ggAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTg                                              <  2:475884/149‑1 (MQ=255)                                                                                                       gTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCg                                         <  1:234599/151‑1 (MQ=255)                                                                                                       gTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCg                                         <  2:479595/151‑1 (MQ=255)                                                                                                            aCGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCa                                    <  2:295968/151‑1 (MQ=255)                                                                                                            aCGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCa                                    <  1:301715/151‑1 (MQ=255)                                                                                                              gCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGt                                  <  1:276818/151‑1 (MQ=255)                                                                                                                    tGCCCTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGAtt                            >  1:576597/1‑151 (MQ=255)                                                                                                                       ccTTAACGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCg                         >  1:372486/1‑151 (MQ=255)                                                                                                                           aaCGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTgg                     <  1:113105/151‑1 (MQ=255)                                                                                                                           aaCGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTg                      <  1:580587/150‑1 (MQ=255)                                                                                                                           aaCGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTg                      <  2:328019/150‑1 (MQ=255)                                                                                                                            aCGCCGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGc                    >  1:294688/1‑151 (MQ=255)                                                                                                                               ccGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGt                  >  1:47283/1‑150 (MQ=255)                                                                                                                               ccGACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGt                  >  2:391547/1‑150 (MQ=255)                                                                                                                                 gACATTCTTCAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAgg               >  1:353575/1‑151 (MQ=255)                                                                                                                                         tcAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCTGATTTCGGTGGCGTAGGattat          <  1:71176/148‑1 (MQ=255)                                                                                                                                          cAGACTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTAc       >  2:224349/1‑150 (MQ=255)                                                                                                                                             aCTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAAc   >  1:53662/1‑151 (MQ=255)                                                                                                                                             aCTCTGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAAc   <  2:578654/151‑1 (MQ=255)                                                                                                                                              ctctGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACg  >  2:555144/1‑151 (MQ=255)                                                                                                                                              ctctGGCGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACg  >  2:272401/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ACGGACGCACGCAAAAAGCGCCGGTGGCTTACTGGAACAAGCGTCACGTAGAGCCGATGCCCGGCAGCATTATTTATGTTGGCCTCGCGGACTCCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCCCATTGGTCCGTCGCAGTCGGATTTCGGTGGCGTAGGATTATTACAAACG  >  NC_000913/4237468‑4237759

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