Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,742,794	G→A	E138K (GAA→AAA)	dgcN →	diguanylate cyclase DgcN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,742,794	0	G	A	100.0%	181.4 / NA	47	E138K (GAA→AAA)	dgcN	diguanylate cyclase DgcN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/19);  total (28/19)

CTTACAGTCTGGAAGCAGCTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAA  >  NC_000913/2742564‑2743015                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              cTTACAGTCTGGAAGCAGCTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTg                                                                                                                                                                                                                      <  1:33415/240‑1 (MQ=255)     cAGTCTGGAAGCAGCTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGc                                                                                                                                                                                                                    >  2:27367/1‑238 (MQ=255)          tGGAAGCAGCTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTAc                                                                                                                                                                                                             >  2:228028/1‑240 (MQ=255)                  gcTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCAt                                                                                                                                                                                                                        >  2:564167/1‑221 (MQ=255)                  gcTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAAcc                                                                                                                                                                                                     >  1:167998/1‑240 (MQ=255)                  gcTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAAcc                                                                                                                                                                                                     >  2:531395/1‑240 (MQ=255)                  gcTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAACTACGCTTAAcc                                                                                                                                                                                                     >  1:453809/1‑240 (MQ=255)                  gcTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGAAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAAcc                                                                                                                                                                                                     >  1:377302/1‑240 (MQ=255)                         tttttGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTcgcg                                                                                                                                                                                              <  1:117976/240‑1 (MQ=255)                          ttttGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGa                                                                                                                                                                                             <  2:137320/240‑1 (MQ=255)                                          gcagcaACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCa                                                                                                                                                                              <  1:531938/239‑1 (MQ=255)                                           cagcaACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCAt                                                                                                                                                                             <  1:212030/239‑1 (MQ=255)                                                aCTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTAt                                                                                                                                                                        >  2:217832/1‑239 (MQ=255)                                                aCTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTAt                                                                                                                                                                        >  2:207805/1‑239 (MQ=255)                                                              aGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcgc                                                                                                                                                         >  1:50513/1‑240 (MQ=255)                                                              aGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcg                                                                                                                                                          <  1:411733/239‑1 (MQ=255)                                                                cgcTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGcc                                                                                                                                                        <  2:342319/239‑1 (MQ=255)                                                                  cTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGt                                                                                                                                                      >  1:134063/1‑239 (MQ=255)                                                                     ggCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTg                                                                                                                                                  <  2:417600/240‑1 (MQ=255)                                                                        cAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGAc                                                                                                                                                <  1:299990/239‑1 (MQ=255)                                                                                         aaCTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTgg                                                                                                                               >  2:267445/1‑239 (MQ=255)                                                                                          aCTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGAc                                                                                                                                                <  1:61130/221‑1 (MQ=255)                                                                                                  aGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCAt                                                                                                                     <  1:406326/240‑1 (MQ=255)                                                                                                        tGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAAt                                                                                                               >  1:142117/1‑240 (MQ=255)                                                                                                                  caAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCAcc                                                                                                     >  1:342016/1‑240 (MQ=255)                                                                                                                                aTCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATCACCCTCACCCGCCAt                                                                                               <  1:561190/232‑1 (MQ=255)                                                                                                                                 tCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAAt                                                                                      >  1:253676/1‑240 (MQ=255)                                                                                                                                         ttACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGcaca                                                                                        <  1:211716/230‑1 (MQ=255)                                                                                                                                                 gTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCCCAAt                                                                                      >  2:361548/1‑224 (MQ=255)                                                                                                                                                     ggATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGaa                                                                  <  1:156224/240‑1 (MQ=255)                                                                                                                                                           aGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAata                                                             >  2:279944/1‑239 (MQ=255)                                                                                                                                                                aTACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGt                                                                                                                                         >  2:576471/1‑158 (MQ=255)                                                                                                                                                                 ttcTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCCTCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCAt                                                                                                                     <  2:274345/175‑1 (MQ=255)                                                                                                                                                                    tttCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGt                                                   <  1:93342/240‑1 (MQ=255)                                                                                                                                                                      tCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTc                                                  >  1:540511/1‑239 (MQ=255)                                                                                                                                                                       cAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCg                                                 <  2:563401/239‑1 (MQ=255)                                                                                                                                                                               tCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACatgat                                         >  1:41183/1‑239 (MQ=255)                                                                                                                                                                                  tAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTcg                                                                                                                                                                                                >  1:218048/1‑85 (MQ=255)                                                                                                                                                                                               ctTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCAc                                                         <  2:578900/207‑1 (MQ=255)                                                                                                                                                                                                ttCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAAt                       >  1:330477/1‑240 (MQ=255)                                                                                                                                                                                                    cTGCCCCCATCATTCAGCCGATTCGTCACAATGGTAAACCCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctc                                                                                                                                                           >  1:600303/1‑104 (MQ=255)                                                                                                                                                                                                       cccccATCATTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTccc                >  2:23653/1‑240 (MQ=255)                                                                                                                                                                                                            atcatTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg            >  1:102709/1‑239 (MQ=255)                                                                                                                                                                                                            atcatTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg            >  1:312611/1‑239 (MQ=255)                                                                                                                                                                                                            atcatTCAGCCGATTCGTCACAATGGTAAAAACATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGaa                                                                          >  2:618048/1‑177 (MQ=255)                                                                                                                                                                                                              catTCAGCCGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAg          <  1:27368/239‑1 (MQ=255)                                                                                                                                                                                                                     ccGATTCGTCACAATGGTAAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaa  >  2:496264/1‑240 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              CTTACAGTCTGGAAGCAGCTGTCGTTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAA  >  NC_000913/2742564‑2743015

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