Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,407,641	A→G	V52A (GTC→GCC)	ydaM ←	diguanylate cyclase, csgD regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,407,641	0	A	G	100.0%	231.4 / NA	64	V52A (GTC→GCC)	ydaM	diguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/27);  total (37/27)

CACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCG  >  NC_000913/1407493‑1407788                                                                                                                                                     |                                                                                                                                                    cACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCg                                                                                                                                                   >  2:220256/1‑151 (MQ=255) cACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCg                                                                                                                                                   <  2:840073/151‑1 (MQ=255)  aCCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCgg                                                                                                                                                  >  2:704528/1‑151 (MQ=255)       ttctgtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAAt                                                                                                                                               >  1:373871/1‑149 (MQ=255)           gtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTc                                                                                                                                         <  2:1186483/151‑1 (MQ=255)           gtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTc                                                                                                                                         <  2:893049/151‑1 (MQ=255)                             aaaaCAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTa                                                                                                                        >  2:338878/1‑150 (MQ=255)                                   aGCCGATAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTg                                                                                                                  >  1:1120667/1‑150 (MQ=255)                                     ccGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAgg                                                                                                               >  2:265502/1‑151 (MQ=255)                                     ccGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAgg                                                                                                               >  1:606825/1‑151 (MQ=255)                                           gCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTg                                                                                                          >  1:1160456/1‑150 (MQ=255)                                             tCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCAt                                                                                                       >  2:225320/1‑151 (MQ=255)                                               aaTGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATc                                                                                                      >  1:701604/1‑150 (MQ=255)                                               aaTGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATc                                                                                                      >  1:632905/1‑150 (MQ=255)                                                 tGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGGTCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAAc                                                                                                   >  2:37076/1‑151 (MQ=255)                                                    agTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGt                                                                                                >  2:358927/2‑151 (MQ=255)                                                      tttCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGGGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGttt                                                                                              >  1:108239/1‑151 (MQ=255)                                                      tttCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGtt                                                                                               <  2:780908/150‑1 (MQ=255)                                                       ttCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGttt                                                                                              <  2:315779/150‑1 (MQ=255)                                                       ttCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGttt                                                                                              <  1:940558/150‑1 (MQ=255)                                                        tCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTc                                                                                             <  1:419341/150‑1 (MQ=255)                                                          ctctTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcgcc                                                                                          >  2:912145/1‑151 (MQ=255)                                                             ttcttTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCcgcg                                                                                       >  2:74055/1‑151 (MQ=255)                                                              tcttTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCcgcg                                                                                       <  2:309470/150‑1 (MQ=255)                                                               cttTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCt                                                                                     >  1:138159/1‑151 (MQ=255)                                                                tttACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTACACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCt                                                                                     >  1:828373/1‑150 (MQ=255)                                                                    cGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTg                                                                                 <  1:1182699/150‑1 (MQ=255)                                                                       tGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCAt                                                                             >  2:840335/1‑151 (MQ=255)                                                                        gAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCAt                                                                             <  2:972355/150‑1 (MQ=255)                                                                        gAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATa                                                                            <  2:1178881/151‑1 (MQ=255)                                                                           cAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAAt                                                                         <  1:566399/151‑1 (MQ=255)                                                                                 ggATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTg                                                                   <  2:1160715/151‑1 (MQ=255)                                                                                  gATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGt                                                                  >  2:886950/1‑151 (MQ=255)                                                                                   aTTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGtt                                                                 <  1:747764/151‑1 (MQ=255)                                                                                    tttCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTc                                                                >  2:1074033/1‑151 (MQ=255)                                                                                      tCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTcc                                                               <  2:1105621/150‑1 (MQ=255)                                                                                         ataataTCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCaa                                                           <  2:820845/151‑1 (MQ=255)                                                                                         ataataTCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCa                                                            >  2:39843/1‑150 (MQ=255)                                                                                             tatCGTGTTGGTTTTGCAGATCGGGGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAggg                                                       >  1:1096528/1‑151 (MQ=255)                                                                                             tatCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAggg                                                       >  1:983836/1‑151 (MQ=255)                                                                                              atCGTGTTGGTTTTGCAGATCTGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGATGCGGAATAGCGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCAACAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGc                                                      <  1:727721/151‑1 (MQ=255)                                                                                               tCGTGTTTGTTTTGCAGAACGGAGAGGTATTTCGGCAGTTCTTTTTGTGAAGAGGCCGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAgg                                                        <  2:1107681/148‑1 (MQ=255)                                                                                                cGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCg                                                     >  2:134404/1‑150 (MQ=255)                                                                                                   gtTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGaaa                                                  <  1:693044/150‑1 (MQ=255)                                                                                                     tGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACg                                                <  2:98677/150‑1 (MQ=255)                                                                                                        ttttGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATcc                                            >  1:306509/1‑151 (MQ=255)                                                                                                         tttGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATcc                                            >  2:1045073/1‑150 (MQ=255)                                                                                                          ttGCAGATCGGAGAGGTATTACGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCa                                           <  1:558708/150‑1 (MQ=255)                                                                                                               gATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAgg                                     >  1:434777/1‑151 (MQ=255)                                                                                                               gATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAgg                                     >  2:877219/1‑151 (MQ=255)                                                                                                                 tCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGAc                                   >  1:148589/1‑151 (MQ=255)                                                                                                                 tCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGAc                                   >  1:337025/1‑151 (MQ=255)                                                                                                                      gagGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTg                              <  2:437801/151‑1 (MQ=255)                                                                                                                        ggTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAg                            <  1:901930/151‑1 (MQ=255)                                                                                                                             ttCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAg                        <  1:144863/150‑1 (MQ=255)                                                                                                                             ttCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGt                       >  1:521338/1‑151 (MQ=255)                                                                                                                              tCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGt                       <  1:357693/150‑1 (MQ=255)                                                                                                                               cGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTc                      <  1:1090845/150‑1 (MQ=255)                                                                                                                                  cAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAgg                  >  2:158211/1‑151 (MQ=255)                                                                                                                                   aGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAggg                 >  1:296916/1‑151 (MQ=255)                                                                                                                                      tCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTAt              >  1:673932/1‑151 (MQ=255)                                                                                                                                      tCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTAt              >  1:137462/1‑151 (MQ=255)                                                                                                                                              tgAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtg      >  2:35750/1‑151 (MQ=255)                                                                                                                                                 ggagGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtgtg    <  2:828561/150‑1 (MQ=255)                                                                                                                                                  gaggCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCg  <  1:766761/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCG  >  NC_000913/1407493‑1407788

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