Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,615,786	T→G	E427A (GAA→GCA)	yjjI ←	DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,615,786	0	T	G	100.0%	218.8 / NA	60	E427A (GAA→GCA)	yjjI	DUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (24/36);  total (24/36)

ATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACT  >  NC_000913/4615638‑4615933                                                                                                                                                     |                                                                                                                                                    atatTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTgcg                                                                                                                                                   <  1:439922/151‑1 (MQ=255)  tatTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTgcg                                                                                                                                                   >  1:270395/1‑150 (MQ=255)   atTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTgcgtc                                                                                                                                                 <  1:223984/151‑1 (MQ=255)        tCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGc                                                                                                                                             <  1:205752/150‑1 (MQ=255)            aaTCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAg                                                                                                                                        <  2:110668/151‑1 (MQ=255)                 gACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGt                                                                                                                                   <  2:375052/151‑1 (MQ=255)                    aaaCGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGc                                                                                                                                <  1:261377/151‑1 (MQ=255)                     aaCGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCt                                                                                                                               <  2:508512/151‑1 (MQ=255)                       cGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGa                                                                                                                             <  2:377402/151‑1 (MQ=255)                        gCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGa                                                                                                                             >  1:875359/1‑150 (MQ=255)                        gCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGAt                                                                                                                            >  2:551828/1‑151 (MQ=255)                        gCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGAt                                                                                                                            <  1:420322/151‑1 (MQ=255)                         cACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATg                                                                                                                           <  2:178830/151‑1 (MQ=255)                         cACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATg                                                                                                                           <  1:572235/151‑1 (MQ=255)                           ccATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGcc                                                                                                                         <  2:438250/151‑1 (MQ=255)                             atatAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCgg                                                                                                                       <  1:242914/151‑1 (MQ=255)                               ataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGaa                                                                                                                     <  2:239621/151‑1 (MQ=255)                               ataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGaa                                                                                                                     <  1:862708/151‑1 (MQ=255)                                 aaCCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGaat                                                                                                                    >  1:379657/1‑150 (MQ=255)                                 aaCCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGaat                                                                                                                    <  1:878210/150‑1 (MQ=255)                                  aCCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGaata                                                                                                                   >  1:500580/1‑150 (MQ=255)                                    cGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataat                                                                                                                 >  1:89269/1‑150 (MQ=255)                                       tAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataataat                                                                                                              <  1:663759/150‑1 (MQ=255)                                         aCGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAataataata                                                                                                             >  2:773336/1‑149 (MQ=255)                                         aCGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAg                                                                                                           <  1:36773/151‑1 (MQ=255)                                             gAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatg                                                                                                       <  2:548191/151‑1 (MQ=255)                                              aaCCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatg                                                                                                       >  1:339950/1‑150 (MQ=255)                                                ccAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatgatg                                                                                                    <  2:529033/151‑1 (MQ=255)                                                   gATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCCCCGAGGAAAAGCTGTACCAGTGCCTCCGGGTTACGTTTGATGGTTGGGTCGAGCGTCAGAATCTCGCTGATGCCCGGATAATAATAAGCATGATGCgg                                                                                                 >  2:798728/1‑151 (MQ=255)                                                      cGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCgggg                                                                                               <  1:878588/150‑1 (MQ=255)                                                        ttACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCg                                                                                             >  2:10807/1‑150 (MQ=255)                                                            cGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGAGCCTGCGGGGTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAg                                                                                         >  1:843618/1‑150 (MQ=255)                                                                gACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTg                                                                                     <  1:814888/150‑1 (MQ=255)                                                                gACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTg                                                                                     >  1:712090/1‑150 (MQ=255)                                                                        cGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGt                                                                            <  1:771917/151‑1 (MQ=255)                                                                          gTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTg                                                                           >  2:496447/1‑150 (MQ=255)                                                                           tAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGAt                                                                         >  2:251377/1‑151 (MQ=255)                                                                            aaaTTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATc                                                                        <  2:622621/151‑1 (MQ=255)                                                                               ttCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCgg                                                                      >  1:302045/1‑150 (MQ=255)                                                                                 cACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGAt                                                                    <  1:798961/150‑1 (MQ=255)                                                                                         ccGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCAt                                                           >  2:263189/1‑151 (MQ=255)                                                                                              tttAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGCGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGc                                                        >  2:476923/1‑149 (MQ=255)                                                                                                 aaaGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATAc                                                   >  1:828288/1‑151 (MQ=255)                                                                                                  aaGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACg                                                  <  2:205974/151‑1 (MQ=255)                                                                                                      cACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCaa                                              <  1:731764/151‑1 (MQ=255)                                                                                                        ccGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAAc                                            <  2:730616/151‑1 (MQ=255)                                                                                                            ggCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                        <  2:497216/151‑1 (MQ=255)                                                                                                              cAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                      >  1:352965/1‑151 (MQ=255)                                                                                                              cAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                      >  1:264078/1‑151 (MQ=255)                                                                                                              cAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                      <  2:474738/151‑1 (MQ=255)                                                                                                              cAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                      <  2:662204/151‑1 (MQ=255)                                                                                                              cAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                      >  1:757599/1‑151 (MQ=255)                                                                                                                  gCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAGTAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCgg                                  <  1:197390/151‑1 (MQ=255)                                                                                                                  gCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCgg                                  >  1:238532/1‑151 (MQ=255)                                                                                                                     gTACCAGTGCCTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGc                                 <  2:213122/149‑1 (MQ=255)                                                                                                                              ccTGCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                      <  1:651600/151‑1 (MQ=255)                                                                                                                                 gCGGGTTACGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATc                   <  2:84183/151‑1 (MQ=255)                                                                                                                                        aCGTTTGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTg            <  1:649105/151‑1 (MQ=255)                                                                                                                                            ttGATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATc         <  2:302037/150‑1 (MQ=255)                                                                                                                                              gATGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccc      >  1:195636/1‑151 (MQ=255)                                                                                                                                               aTGGTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg     >  2:247282/1‑151 (MQ=255)                                                                                                                                                 ggTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGAc   >  2:252341/1‑151 (MQ=255)                                                                                                                                                  gTTGCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACt  <  2:270390/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACT  >  NC_000913/4615638‑4615933

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