Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,969,912	C→T	G358D (GGC→GAC)	tap ←	methyl‑accepting chemotaxis protein Tap

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,969,912	0	C	T	100.0%	222.3 / NA	68	G358D (GGC→GAC)	tap	methyl‑accepting chemotaxis protein Tap
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/28);  total (40/28)

CGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT  >  NC_000913/1969667‑1970156                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      cGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGCAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAAtt                                                                                                                                                                                                                                                  >  2:215793/1‑251 (MQ=255)  aaTCATGGTCGCGGCGGCGTTTTTCACCAGTTTCGAACCCTGCGGGACACGATTGACTGACTCTTCAAGCGCCCCTTTGATCTCTTTTGCGGCCGGCGCGCTACGGCTGGAAAGTTGGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAAttt                                                                                                                                                                                                                                                 <  1:630761/250‑1 (MQ=255)   aTCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAAtttt                                                                                                                                                                                                                                                >  1:547743/1‑251 (MQ=255)     cATGGTCGCGGCGGCGTTATTCACTAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAAtttttt                                                                                                                                                                                                                                              >  2:47877/1‑251 (MQ=255)                  gcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGc                                                                                                                                                                                                                                >  2:130137/1‑251 (MQ=255)                                tttCGAACCCTCCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCCAGGTTGCGCCCTTCACCTGCCACTACCGCAAAACCACGCACCTGCTCTCCGTCGCGAGCCGCTTCAACTACCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAATCCGTCGCTAACGCTGATAATGTCGTCAATTTTTTGAGAG‑TGTGAGGCGATCGCCTGCATgg                                                                                                                                                                                                                  >  2:222439/1‑251 (MQ=255)                                                                ctTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACAccgccc                                                                                                                                                                                  >  1:487391/1‑251 (MQ=255)                                                                    gATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCt                                                                                                                                                                              >  2:563358/1‑251 (MQ=255)                                                                    gATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCt                                                                                                                                                                              >  1:410045/1‑251 (MQ=255)                                                                           ccTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGt                                                                                                                                                                       >  1:529339/1‑251 (MQ=255)                                                                             tttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAac                                                                                                                                                                                         >  1:226419/1‑231 (MQ=255)                                                                             tttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTg                                                                                                                                                                     >  1:203166/1‑251 (MQ=255)                                                                                  tctcTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGAcc                                                                                                                                                                                              >  1:36605/1‑221 (MQ=255)                                                                                       tttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCAtt                                                                                                                                                           >  1:137049/1‑251 (MQ=255)                                                                                             gcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGc                                                                                                                                                     >  2:476037/1‑251 (MQ=255)                                                                                                tgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGcc                                                                                                                                                    >  1:598008/1‑249 (MQ=255)                                                                                                       gcggttgtacgaggctcccgctctCACCTGCCACGCCCGCAAACCCCGGCCCCTTCTCTCCGGCGCTGCCCGCTTCCACTGCCGCATCCAGGGCCAAAATATTGGTCGGGAAAGCAATTCCGTCGATAACTCTGATAGTGGCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGAc                                                                                                                                           <  2:482976/228‑1 (MQ=255)                                                                                                        aCGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACg                                                                                                                                          >  1:402154/1‑251 (MQ=255)                                                                                                               gcaagttccgccatccccccTGCCCCTCCCGCAAAACCACGCCCCTGCTCCCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAGTGGTCTGGAAAGCAATTCCGTCGATACCGCTGATATTGTCTTCATTTTTTTGCGATCTG‑GTGGCGATCTCCTGCATGGTTGGAGTCATGGTACTGACCGAAACGCCGCCCGCCTGCGCCGTTGTCGCTGCATTTTTTGCCAGTTCCGACGCCTGTCg                                                                                                                                   <  1:47629/233‑1 (MQ=255)                                                                                                                cAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCgg                                                                                                                                  >  1:272786/1‑250 (MQ=255)                                                                                                                 aaGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTAGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGGTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTg                                                                                                                                 >  1:318547/1‑251 (MQ=255)                                                                                                                   gATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCCGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACAGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCg                                                                                                                               >  1:11905/1‑251 (MQ=255)                                                                                                                         gcACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCg                                                                                                                         >  2:444994/1‑251 (MQ=255)                                                                                                                          cACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTAGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCgg                                                                                                                        >  1:343822/1‑251 (MQ=255)                                                                                                                           ccccccccTGCCCCTACCTAAAAACCACCCCCCTGCTCTCCGCCCCGATCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGCCCTGAACTCCTCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGc                                                                                                                       <  1:344130/245‑1 (MQ=255)                                                                                                                            cccccccgccccctcccccAAACCCACCCCCCTGCTCCCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTGGCGAGCTG‑GTGGCGATCTCCTGCAGGGTGTGAGTAATGGTACTGACCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCg                                                                                                                      <  1:156322/233‑1 (MQ=255)                                                                                                                            cTTCACCTGCCACTACTGCAAAAGCACGCCCCTGCTCTCCGGCGCGAGCCGCTTACACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GAGGCGATCTCCTGCATGGTGTGAGTCATGGGACGGACCTGAACACCGCCCGCCTGCGCCGTTGTCGAGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCg                                                                                                                      >  2:468283/1‑251 (MQ=255)                                                                                                                              cccccTGCCACTACCCCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGt                                                                                                                     <  1:444994/247‑1 (MQ=255)                                                                                                                                  cTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGAACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGc                                                                                                                >  2:439201/1‑251 (MQ=255)                                                                                                                                       ccTACCGCAAACCCACGCCCCGGCTCCCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTAc                                                                                                           <  1:71863/250‑1 (MQ=255)                                                                                                                                           ccGCAAAACCACGCCCCTGCTCCCGGGCGCGCGCCGCTTCCACTGCCGCATGCAGTGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATACCGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTg                                                                                                       <  2:340821/251‑1 (MQ=255)                                                                                                                                                    cACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGCATCTTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAAt                                                                                              >  2:140214/1‑251 (MQ=255)                                                                                                                                                       gCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGCTTATTGc                                                                                           >  2:212542/1‑251 (MQ=255)                                                                                                                                                                 tcCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGGTTATTGCTCCATACTgg                                                                                 >  2:192382/1‑251 (MQ=255)                                                                                                                                                                    ggCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTg                                                                                                                 <  1:615596/216‑1 (MQ=255)                                                                                                                                                                    ggCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTg                                                                                                                 >  2:615596/1‑216 (MQ=255)                                                                                                                                                                     gcgcGAGCCGCTCCCACTGCCGCATTCAGGGCCAAAATATTGTTCTGAAAAGCATTTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTggcggc                                                                             <  2:487391/251‑1 (MQ=255)                                                                                                                                                                      cgcgAGCCGCTTCCCCTGCCGCATTGGGGTCCAGAATATGGGTCTGGAAAGCAATTCCGTCATTAACGCTGAGAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCAGGGGGTGAGTCATGGTACTGACCTGAACACCTCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGTCGGcc                                                                            <  2:272786/251‑1 (MQ=255)                                                                                                                                                                       gcgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTTGCGGCCg                                                                           >  1:132923/1‑251 (MQ=255)                                                                                                                                                                                 tCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCa                                                                  <  1:47877/250‑1 (MQ=255)                                                                                                                                                                                  ccACTGCCGCATTCAGGTCCAGAATATGGGGCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCATTTTGTGGCGAGTTG‑GGGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCGGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGc                                                                <  1:140214/251‑1 (MQ=255)                                                                                                                                                                                  ccACTGCCGCATTCAGGGCAAAAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGTTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGc                                                                <  2:331142/251‑1 (MQ=255)                                                                                                                                                                                    ccTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTGGAAAACGTGGATAATGTCGTCATTTTTTTGCGAGCGG‑GGGCCGATCTCCTGCATGGTGTGAGCGAGGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTGCCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTCCCAGCGa                                                              <  2:267828/250‑1 (MQ=255)                                                                                                                                                                                             ttCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTTGCCAGCGATgccgcccc                                                     >  2:626527/1‑249 (MQ=255)                                                                                                                                                                                                  ggCCAGAAGATTGGTCTGGAAAGCAATTCCGTCGAAAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGtt                                                <  2:137049/251‑1 (MQ=255)                                                                                                                                                                                                   gcaggAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCATTTTTTGGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGTTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTc                                               <  1:62761/247‑1 (MQ=255)                                                                                                                                                                                                       gAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTa                                           <  1:215793/251‑1 (MQ=255)                                                                                                                                                                                                         attttGGTCTGGAAAGCAATTCCGTCGATAACGCGGATAATGTCGTCAATTTTTTGCGAGCGG‑GTGGCGATCCCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACg                                         <  2:537040/248‑1 (MQ=255)                                                                                                                                                                                                              ggTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTTCTCCATACTGGCGGCCGTTTGTGCCCGCGTTGCCCCCTGCTGTTCCGgacgacctt                                    >  2:26108/1‑248 (MQ=255)                                                                                                                                                                                                              ggTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGGTGCCGCCTGCTGGTCGGTACGACttt                                    >  1:274204/1‑250 (MQ=255)                                                                                                                                                                                                              ggTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACtt                                     <  1:543351/250‑1 (MQ=255)                                                                                                                                                                                                                 cTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACAccgcccg                                                                                                                                                                                 <  2:642705/107‑1 (MQ=255)                                                                                                                                                                                                                 cTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACAccgcccg                                                                                                                                                                                 >  1:642705/1‑107 (MQ=255)                                                                                                                                                                                                                       aGCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTTCATACTTGCGGGCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGGACGACTTGAGAGATCGt                           >  1:374635/1‑251 (MQ=255)                                                                                                                                                                                                                        gCAATTCCGTCGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGtt                          <  1:523427/251‑1 (MQ=255)                                                                                                                                                                                                                                  cGATAACGCTGATAATGTCTTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGTTACTGCCCTGACCACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGc                 <  2:58051/250‑1 (MQ=255)                                                                                                                                                                                                                                  cGATAACGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGTGAGATCGTTATTGCCTGCg                >  1:617298/1‑251 (MQ=255)                                                                                                                                                                                                                                   gaaaaCGCTGAAAATGTCGTAAATTTTTTGCGAGGTG‑GTGGCGAGCTCCTGCTTGGTGTGAGTCATGGTACGGACCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGa               <  2:547743/248‑1 (MQ=255)                                                                                                                                                                                                                                      aaCGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCgg                                                                                                                        <  1:55870/143‑1 (MQ=255)                                                                                                                                                                                                                                      aaCGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCgg                                                                                                                        >  2:55870/1‑143 (MQ=255)                                                                                                                                                                                                                                      aaCGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGGTTGTGCCCGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACa             >  1:439503/1‑250 (MQ=255)                                                                                                                                                                                                                                       aCGCTGATAATGTCGTCTATTTTTTGCGTGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAAt           <  2:388114/251‑1 (MQ=255)                                                                                                                                                                                                                                       aCGCTGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGCTCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACAGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGGTAATTGCTCCATTCTGGCGGCCGCTTGTGCCAGCGATGCCGCCTGTTGGTCGGTACCACCTGTGAGATCGGTATTGCCTCCGACAAt           >  2:95946/1‑251 (MQ=255)                                                                                                                                                                                                                                           tGATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCAGGGTGTGAGTCATGGTACTGCCCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAAtctc        <  1:224014/250‑1 (MQ=255)                                                                                                                                                                                                                                            gATAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGTTCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTcc       <  2:95237/250‑1 (MQ=255)                                                                                                                                                                                                                                             aTAATGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGAGGCCGCCTGCTGTTCGGTACGACTTGCGAGATCGTTATTGCCTGCGACAATCTCcgc     >  1:367096/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGGGGCCGTTTGTGCCAGCGATGCCGCCTGCGGTTCGGTACGACTTTGGAGATCGTTATTGCCTGCGACAATCTCCCCGAt  >  2:120451/1‑250 (MQ=255)                                                                                                                                                                                                                                                 tGTCGTCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGGGGCCGGTTTTTCCAGCGAGGGCGCCTGCTTTTTGGTACGACTTGGGAGATAGtttatgtcctcgacaacctccgcgca  >  2:569866/1‑226 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTG‑GTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT  >  NC_000913/1969667‑1970156

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