Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,969,982	C→T	A335T (GCG→ACG)	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,982	0	C	T	100.0%	235.8 / NA	70	A335T (GCG→ACG)	tap	methyl‑accepting chemotaxis protein Tap
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/26);  total (44/26)

CAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCT  >  NC_000913/1969737‑1970227                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       cAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTgcg                                                                                                                                                                                                                                                   >  1:185747/1‑251 (MQ=255)          gATCTTTTTTGCTCCCGGCGCGCCACGGGTGGCAAGTTGGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTc                                                                                                                                                                                                                                           <  1:66714/250‑1 (MQ=255)                   tgctgcCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGCGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGACCACCGCCCGTCTGCGCCGTTGTCGCGGCAtttt                                                                                                                                                                                                                                 >  1:432694/1‑251 (MQ=255)                         ctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGcc                                                                                                                                                                                                                            >  1:176167/1‑250 (MQ=255)                          tgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACAACAAAAACACGCCCATGCTCT‑CAGGCGCAACACGCTTCCCCTGACGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTCGCGAGCTGGCGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGAATTTTTTGca                                                                                                                                                                                                                            >  1:43485/1‑248 (MQ=255)                               gcTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCg                                                                                                                                                                                                                    >  1:44199/1‑251 (MQ=255)                                cTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCCTGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCgc                                                                                                                                                                                                                   >  2:414475/1‑250 (MQ=255)                                     gCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCACCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGACAGT‑TCCGACGCCt                                                                                                                                                                                                              >  1:514993/1‑251 (MQ=255)                                      cTGGCAAGCTTGCGCACTTCAGCTGCCACTACTTCAAAACCACGCCCCTGCTCT‑CCGGCGCGCGCCGCTTCCACTGCCGCATTCAGGGCCAGAATCTTGGTCTAGAAAGCAATTCCGCCGATACCGCTGCGAATGACGCCATTTTTTTGTGAGCTGGTGGCGATCTCCTGCATAGTGTGAGTCATGGTACTGCCCTGACGACCGCCCGTCTGCGCCGTTGTCACGGCATTCTT‑GCCAATGTCCGACGccgg                                                                                                                                                                                                             >  2:320883/1‑249 (MQ=255)                                       tGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTAT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGCGCCAGCATATTGGTCTGGAAAGCTATTCCGTCGATAACGCTGCTAATGTCGCCACTTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGCGTCATGGTACTGACCTGGACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCATGt                                                                                                                                                                                                            >  1:296247/1‑251 (MQ=255)                                        cctgaaaggcccctcttccgcTGCCCCTACCGCACAACCACGCCCCGGCTCC‑CCGGCGCGAGCCGCGTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCTTCGATAACGCTGATAATTTCGCCAATTTTTTGCGATCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGTCCTGACCACCGCCCGTCTGCGCCGTTGACGCGGCATTTTTTGCCAGT‑TCCGATGCCTGTc                                                                                                                                                                                                           <  2:74608/231‑1 (MQ=255)                                            atatggcgccctccACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGc                                                                                                                                                                                                       <  1:114478/238‑1 (MQ=255)                                            agttggcccccccccccTGCCCCTACCGCAAACCCACGCCCCTGC‑CTGCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAGCACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGGCAGG‑TCCGACGCCtcgcgtgc                                                                                                                                                                                                       <  2:78865/236‑8 (MQ=255)                                                    cACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCAGCATTCAGGGCCCGAATATTGGTCTGGAAAGCAATTCCGTAGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCgg                                                                                                                                                                                               >  1:556511/1‑251 (MQ=255)                                                     ccccccccTGCCACTCCCGCAAAACCCCGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGc                                                                                                                                                                                              <  2:166822/245‑1 (MQ=255)                                                        tCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGc                                                                                                                                                                                              >  1:166822/1‑248 (MQ=255)                                                             tGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGcc                                                                                                                                                                                      >  1:21761/1‑251 (MQ=255)                                                                       gCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGc                                                                                                                                                                            <  1:160928/251‑1 (MQ=255)                                                                        cAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCg                                                                                                                                                                           >  1:474974/1‑251 (MQ=255)                                                                         aaaaCCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCgg                                                                                                                                                                          >  1:9640/1‑251 (MQ=255)                                                                          aaaCCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCAATCAGGGCCAGAATCTTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTCCCGTGGCGGt                                                                                                                                                                         >  2:425894/1‑251 (MQ=255)                                                                                  ccccTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATCTTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCACTTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCACGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTTAtttct                                                                                                                                                                 >  2:574259/1‑248 (MQ=255)                                                                                    ccTGCTCC‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCTACGCCTGTCGTGCGTTATCGGCGTTTTGCCGTACCGTCGCGTTTAATTGCTcc                                                                                                                                                               <  2:381508/251‑1 (MQ=255)                                                                                      tGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTccct                                                                                                                                                             >  2:435918/1‑249 (MQ=255)                                                                                               ggCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTggcgg                                                                                                                                                     <  2:109035/251‑1 (MQ=255)                                                                                                  gcgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATACCGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGCTCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTTCTCCCTACCGGCGGCCg                                                                                                                                                  >  2:70121/1‑251 (MQ=255)                                                                                                   cgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATACTGTCGCCAATTTTTTGCGAGCTGGTGGCGCTCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGCCGt                                                                                                                                                 >  1:357976/1‑251 (MQ=255)                                                                                                        cGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGGTGATTGCTCCATACTGGCGGCCGtttctg                                                                                                                                            >  2:10440/1‑248 (MQ=255)                                                                                                          cTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTTCTCCATACTGGCGGCCGTTTGTGcc                                                                                                                                          >  1:208739/1‑251 (MQ=255)                                                                                                               aCTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATTCTGGCGGCCGGTTGTGCCAGCgg                                                                                                                                     >  2:450089/1‑250 (MQ=255)                                                                                                                   ccGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATgcc                                                                                                                                 <  1:548077/251‑1 (MQ=255)                                                                                                                       aTTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTTATTGCTCCATAATGGCGGCCGTTTTTGCCAGCGATTCCGCct                                                                                                                             >  2:178133/1‑251 (MQ=255)                                                                                                                             ggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCCGCTGtt                                                                                                                       >  1:219536/1‑251 (MQ=255)                                                                                                                              gCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCCGTTc                                                                                                                      >  2:422083/1‑251 (MQ=255)                                                                                                                                         ggTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATGTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTg                                                                                                           <  2:185747/251‑1 (MQ=255)                                                                                                                                         ggTCTGGAAAGCAATTCCGTCGATAACGCGGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACGCCGCCCGTCTGCGCCGTTGTCGCGGCAGTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACtt                                                                                                            <  1:380960/250‑1 (MQ=255)                                                                                                                                         ggTCTGGAAAGCAAGTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACtt                                                                                                            <  1:414475/250‑1 (MQ=255)                                                                                                                                            cgggAAAGCAAGTCCGTCGAGAACGCTGATAATGTCGCCAATTTTTTGCGAGCGGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCCGT‑TCCGACGCCTGTCGTGCTTTATCGGCGTTTTGCCCTACCGTGGCGGGTAATTGCTCCAGACTGGCGGCCGTTTGTGCCAGCGATGCCGCTTGCTGTTCGGTACGACTTgaga                                                                                                        <  2:77448/249‑1 (MQ=255)                                                                                                                                              ggAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTgag                                                                                                         <  1:275433/248‑1 (MQ=255)                                                                                                                                              ggAAAGCAAGTCCGTCGATAACGCTGATAATGTCGCCAATGTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTgag                                                                                                         <  1:210222/248‑1 (MQ=255)                                                                                                                                               gAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAg                                                                                                                                                                                                                          <  2:204299/135‑1 (MQ=255)                                                                                                                                               gAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAg                                                                                                                                                                                                                          >  1:204299/1‑135 (MQ=255)                                                                                                                                                aaaGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGGCTTgag                                                                                                         >  2:268758/1‑246 (MQ=255)                                                                                                                                                   gCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGt                                                                                                                                                                                                            >  1:518196/1‑144 (MQ=255)                                                                                                                                                   gCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGt                                                                                                                                                                                                            <  2:518196/144‑1 (MQ=255)                                                                                                                                                     aaTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTTGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGCCTTGAGAGATCGTTAt                                                                                               >  2:439977/1‑251 (MQ=255)                                                                                                                                                      aTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATCCTGGCGGCCGTTTGTGCCAGCGGAGCCGTCTGCTGTTCGGTACGACTTGGGAGATCGTTAtt                                                                                              >  1:7154/1‑251 (MQ=255)                                                                                                                                                      aTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATTCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTAtt                                                                                              >  1:149134/1‑251 (MQ=255)                                                                                                                                                      aTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTAtt                                                                                              <  1:268758/251‑1 (MQ=255)                                                                                                                                                               aTAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTTTCCAGCCATGCCGCCTTCTTTTCGGTACGCCtttagtgcttgttttttcctgcggc                                                                                     >  2:461209/1‑228 (MQ=255)                                                                                                                                                               aTAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGTCTGCGAc                                                                                     <  1:141706/251‑1 (MQ=255)                                                                                                                                                                   cGCTGATAATGTCGCCACTTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGCTGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAAtc                                                                                 >  1:184518/1‑251 (MQ=255)                                                                                                                                                                      tGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGt                                                                                                                                                                                            >  2:267267/1‑141 (MQ=255)                                                                                                                                                                       gATAATGGCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGTCCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATGTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTcc                                                                              <  2:474879/250‑1 (MQ=255)                                                                                                                                                                         tAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCCGCGATGCCGCCTGCTGTTCGGTACGACTTGAGGGCTCGTTATTGCCTGCTACAATCTCcgcg                                                                           >  1:4713/1‑251 (MQ=255)                                                                                                                                                                           aTGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAt                                                                         <  1:439754/251‑1 (MQ=255)                                                                                                                                                                                gCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTTGCCAGCGATGCCGCCTGCTGTTCGGTACGACCTGAGAGATCGTTATTGCCCGCGACAATCTCCCCGATACCa                                                                     >  2:16552/1‑250 (MQ=255)                                                                                                                                                                                gCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCCGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACCATCTCCGCGATACCa                                                                     >  1:111067/1‑250 (MQ=255)                                                                                                                                                                                        tttGCGTGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTACTTGCTCCATACTGGCGGCCGTTTTTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTTAGAGATCGTTTTTTCCTGCGACAATCTCCCCGCTACCCATGTGCAtc                                                            >  1:60655/1‑251 (MQ=255)                                                                                                                                                                                              aGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTTGCCAGCGATGCCGCCTGCTTTTCGGTACGGCTTTAGAGATCGTTATTGCCTGCGACAATCTCCCCGTTACCCACGTGCCTCTCCTgg                                                      >  1:565875/1‑251 (MQ=255)                                                                                                                                                                                                     ggCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCt                                               <  1:439029/251‑1 (MQ=255)                                                                                                                                                                                                          cctcCTGCAGGGTGTGATTCATGGTACTTACCTGAACACCTCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGGTTGTGCCAGCGATGCCTCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCGCCTCGATACCAATGTGCATCTCCTGGCTTCCCTTACGc                                          <  2:73826/250‑1 (MQ=255)                                                                                                                                                                                                               tGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGGTAATTTCTCCCTACTGGCGGCCGTTTGGGGCAGCGATGCCGCCTGCTGTTTGGTACGGCTTGAGAGAACGTTATTGCCTGGGCCAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACggccatc                                     >  1:233221/1‑245 (MQ=255)                                                                                                                                                                                                                  aTGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTAGCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACt                                  <  2:543206/251‑1 (MQ=255)                                                                                                                                                                                                                   tGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGGGGTCGTTATTGCCTGCGACAATCTCCGCGATACCAACGTGCATCTCCTGGCTTCCCTTCCGCACATCCCtt                                 >  1:503851/1‑251 (MQ=255)                                                                                                                                                                                                                    ggTGTGAGTCATTGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTa                                <  2:44199/251‑1 (MQ=255)                                                                                                                                                                                                                             cATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCAc                       <  1:69436/251‑1 (MQ=255)                                                                                                                                                                                                                              aTGGTACTGACATGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTAGCGGTTAATTGCTCCATACTGGCGGCCGGTTGTGCCAGCGGTGCCGCCTGCTGTTCGGTACGACTTGGGAGATCGTTATTGCCTGCGACAATCTCCGAGATACCAATGGGCATCTCCTGGCTTCCCTTACGCACATCACGTACCGTcccccg                      >  1:453324/1‑248 (MQ=255)                                                                                                                                                                                                                                       aCCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGGGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCCAAGCctg             >  1:254980/1‑251 (MQ=255)                                                                                                                                                                                                                                         cTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCctgct           >  1:337849/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gcccgtcTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGACGTTTTTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTAACGTCCCACCCAAAGCCTGCTGCCTggcct  >  1:560061/1‑248 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       CAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCT‑CCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGT‑TCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCT  >  NC_000913/1969737‑1970227

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: ‑