Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,611,401	T→C	T152A (ACT→GCT)	dgcF ←	putative diguanylate cyclase DgcF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,611,401	0	T	C	100.0%	237.4 / NA	64	T152A (ACT→GCT)	dgcF	putative diguanylate cyclase DgcF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/28);  total (36/28)

CATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGTCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACATTTGCGGAAT  >  NC_000913/1611268‑1611533                                                                                                                                      |                                                                                                                                     caTTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCa                                                                                                                                    >  2:364325/1‑136 (MQ=255) caTTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGaaa                                                                                                                                >  1:638234/1‑140 (MQ=255)      aTGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCaa                                                                                                                           <  2:638234/140‑1 (MQ=255)       tGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCaaa                                                                                                                          >  2:35508/1‑140 (MQ=255)            aTAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCgg                                                                                                                     >  1:1024151/1‑140 (MQ=255)            aTAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCgg                                                                                                                     >  1:350460/1‑140 (MQ=255)                ttATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCg                                                                                                                 >  1:101995/1‑140 (MQ=255)                  aTCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCa                                                                                                               <  1:1086001/140‑1 (MQ=255)                   tCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCg                                                                                                                 <  1:1108168/137‑1 (MQ=255)                   tCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCg                                                                                                                 >  2:1108168/1‑137 (MQ=255)                        aaTGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCaa                                                                                                         >  2:1070386/1‑140 (MQ=255)                          tGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACt                                                                                                       >  2:535974/1‑140 (MQ=255)                          tGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACt                                                                                                       >  2:393580/1‑140 (MQ=255)                          tGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACt                                                                                                       <  1:314594/140‑1 (MQ=255)                          tGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACt                                                                                                       >  1:765878/1‑140 (MQ=255)                          tGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAG‑AAATCAAAGTCGGCTCGCAGCGCAACtt                                                                                                      >  1:798905/1‑140 (MQ=255)                          tGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTAAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCg                                                                                                                 <  2:422811/130‑1 (MQ=255)                          tGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTAAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCg                                                                                                                 >  1:422811/1‑130 (MQ=255)                              tttGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCt                                                                                                   >  2:203235/1‑140 (MQ=255)                              tttGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCt                                                                                                   >  1:743560/1‑140 (MQ=255)                               ttGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCtt                                                                                                  <  1:1070386/140‑1 (MQ=255)                               ttGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCtt                                                                                                  <  1:535974/140‑1 (MQ=255)                                      aGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGc                                                                                           >  1:723272/1‑140 (MQ=255)                                        tCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGa                                                                                         >  1:939420/1‑140 (MQ=255)                                        tCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCG‑AACTTGCTTCATTAGCGaa                                                                                        <  1:423856/140‑1 (MQ=255)                                         cGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGaa                                                                                        <  2:1024151/140‑1 (MQ=255)                                              tCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGAt                                                                                   <  2:939420/140‑1 (MQ=255)                                                   cagcaTGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAG‑AAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGc                                                                             <  2:798905/140‑1 (MQ=255)                                                    agcaTGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGc                                                                             <  2:743560/140‑1 (MQ=255)                                                      caTGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGccg                                                                           >  2:726965/1‑139 (MQ=255)                                                      caTGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCa                                                                           >  2:774960/1‑140 (MQ=255)                                                      caTGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCa                                                                           >  1:386443/1‑140 (MQ=255)                                                        tGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACg                                                                         <  1:203235/140‑1 (MQ=255)                                                          aCAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCt                                                                       <  1:774960/140‑1 (MQ=255)                                                           cAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTa                                                                      >  1:36161/1‑140 (MQ=255)                                                              tCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTaaaa                                                                   <  1:35508/140‑1 (MQ=255)                                                                      tgtgTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGc                                                                        <  1:657564/127‑1 (MQ=255)                                                                      tgtgTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGc                                                                        >  2:657564/1‑127 (MQ=255)                                                                        tgtCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGATCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGGAGCGCAACTTTCTTCATTAGCGAATTGATCGCTGCTACGCTAAAAGAAACCATaa                                                         >  2:770247/1‑140 (MQ=255)                                                                        tgtCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATaa                                                         >  2:757727/1‑140 (MQ=255)                                                                          tCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAAtt                                                       >  2:365800/1‑140 (MQ=255)                                                                            tGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTgg                                                     >  2:339139/1‑140 (MQ=255)                                                                            tGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGCAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTgg                                                     >  2:587679/1‑140 (MQ=255)                                                                               ttCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCa                                                  >  2:41561/1‑140 (MQ=255)                                                                                    cGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCg                                             <  1:757727/140‑1 (MQ=255)                                                                                      aTGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGcc                                           <  2:101995/140‑1 (MQ=255)                                                                                       tGGACTCTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCa                                          >  2:562251/1‑140 (MQ=255)                                                                                        ggACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCGCGCTAAAAGAAACCATAATTGGGCATATCGCCAt                                         <  1:726965/140‑1 (MQ=255)                                                                                        ggACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCAt                                         <  1:339139/140‑1 (MQ=255)                                                                                        ggACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCAt                                         <  2:103804/140‑1 (MQ=255)                                                                                               tCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCa                                  <  1:41561/140‑1 (MQ=255)                                                                                                  gcgcCTCATAAAGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATa                               <  2:723272/140‑1 (MQ=255)                                                                                                            aaGACCGGACCGTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTg                     <  1:393580/140‑1 (MQ=255)                                                                                                                 cggaccgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGgag                >  2:651264/1‑140 (MQ=255)                                                                                                                 cggaccgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGgag                >  2:704302/1‑140 (MQ=255)                                                                                                                   gaccgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGaa              >  1:159417/1‑140 (MQ=255)                                                                                                                   gaccgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGaa              >  1:451201/1‑140 (MQ=255)                                                                                                                   gaccgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGaa              >  1:236307/1‑140 (MQ=255)                                                                                                                    accgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAAc             <  2:950546/140‑1 (MQ=255)                                                                                                                      cgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACAt           <  2:765455/140‑1 (MQ=255)                                                                                                                      cgTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACAt           >  2:988147/1‑140 (MQ=255)                                                                                                                       gTGAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACAtt          >  1:373825/1‑140 (MQ=255)                                                                                                                         gAGTAAACCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACATTTg        <  2:87542/140‑1 (MQ=255)                                                                                                                              aaCCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACATTTGCGGaa   <  1:770247/140‑1 (MQ=255)                                                                                                                               aCCTGAGCCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACATTTGCGGAAt  >  1:243638/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                     CATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGACAGTCAGATGTTGTGTCTGTTTCAGCGATGGACTTTTCAGCGCCTCATAAAGACCGGACCGTGAGTAAACCTGAGTCAGAAAATCAAAGTCGGCTCGCAGCGCAACTTGCTTCATTAGCGAATTGATCGCTGCCACGCTAAAAGAAACCATAATTGGGCATATCGCCATCGTGGCAATACCGAGACGTGCGGAGAACATTTGCGGAAT  >  NC_000913/1611268‑1611533

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