Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,638	T→C	intergenic (‑10/‑190)	iclR ← / → metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,638	0	T	C	100.0%	134.4 / NA	45	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/28);  total (17/28)

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAT  >  NC_000913/4223393‑4223874                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                aGTTCGCCAACCTGACGCACGAAACCCGGCTGTTGCATCGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTATCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAAGCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCcttttt                                                                                                                                                                                                                                            <  1:568905/251‑1 (MQ=255)      gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTa                                                                                                                                                                                                                                       <  1:120414/251‑1 (MQ=255)      gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGTTTTGCGGGGATGGGGGGGACCATGTCAGTATCCTTTTTCTGtt                                                                                                                                                                                                                                       >  2:332567/1‑250 (MQ=255)                  cACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCGCCTTTTTCTGTATCTTGGAAAt‑ca                                                                                                                                                                                                                          <  2:566500/251‑1 (MQ=255)                    cGAATCCCGGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCTGCTTTTTGCGCCCGTTCCGTGCGTTCCACACTGCCATTGGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCTACCATGACAGTCTCCTTTTTCTTTAGCGTGGAAAt‑catt                                                                                                                                                                                                                        <  2:318731/251‑1 (MQ=255)                      aaaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCCGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGTTGCGGTTTGCGGGGATGGGGGCGAACCTGACAGTCTCCTTTTTTTTTATCTGGGGAAt‑catttt                                                                                                                                                                                                                      >  1:195689/1‑251 (MQ=255)                                  tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTTGCGGGGATGGGGGCGGCCCTGACAGGCTCCTTTTTCTTTATCGTGGGAAT‑CATTTTCAtttttttttt                                                                                                                                                                                                          >  1:482575/1‑246 (MQ=255)                                    cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑a                                                                                                                                                                                                       <  1:254447/251‑1 (MQ=255)                                        gTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACGCCAGTAATGTCAGGCCACGCGTTAAAGACTGACCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGTCAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTGCATTTTTATTGTT‑AGCTa                                                                                                                                                                                                   <  1:109747/251‑1 (MQ=255)                                         gggTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCGGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTCTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTaa                                                                                                                                                                                                  <  1:341246/250‑1 (MQ=255)                                             tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTACC‑CGGTTCGCGGGGGAGGGGGGGAACAAGACAGTCTCCTTTTTTTGTATTGGGGGAAT‑CCTTTTCACTTTTTTTGtttggttaagcc                                                                                                                                                                                               >  1:19744/1‑241 (MQ=255)                                                  ggCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGt                                                                                                                                                                                         <  2:482575/251‑1 (MQ=255)                                                  ggCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGGGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGGCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTGCATTTTTATTGTT‑AGCTAATGCAATAGt                                                                                                                                                                                         <  2:298314/251‑1 (MQ=255)                                                                    tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGAt                                                                                                                                                                       <  2:461684/250‑1 (MQ=255)                                                                    tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTGTTCTGTATCGTGGAAAG‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGAt                                                                                                                                                                       <  1:43524/250‑1 (MQ=255)                                                                     gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACATTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTGTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGa                                                                                                                                                                     <  2:465468/251‑1 (MQ=255)                                                                           cccGGCTTGTTGCGCCCGTTCCGTGAGTGCCACACTGCCATTGGAGTCGGCAATCCACTCAAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAAt                                                                                                                                                               <  1:343253/251‑1 (MQ=255)                                                                              ggCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGTTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGtt                                                                                                                                                            <  1:414599/251‑1 (MQ=255)                                                                                        gcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAAtct                                                                                                                                                  <  1:163898/251‑1 (MQ=255)                                                                                          cAGTTCCGTGAGTGCCCCACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGACACGCGTTAAAGACTGAACCTGGTCAGTCGCTGGTGCGGTGGCCACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGGTCATGACAGTCTCCTGTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCa                                                                                                                                                <  1:261920/251‑1 (MQ=255)                                                                                            gTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATg                                                                                                                                              <  2:267469/251‑1 (MQ=255)                                                                                                       gCCACACTGCCATTGGATTCGGCAATCCTCTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGGCGCGTTTCGCGTGAATGGGTGCGATCAGGACAGGCTGCTTTTTCTGGTTCGTTGAAAT‑CATTTTCATTTTTAGTGGT‑AGTTAATGCAATAGTTACGGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGtggtgg                                                                                                                                   <  1:537990/251‑1 (MQ=255)                                                                                                        ccACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCCCGGTTCCCGGGGATGGGGGGGACCCAGGCAGTCTCCTTTTTTTGTATAGGGGAAATTCTTTTTC‑TTTTTTTTTTT‑AACTAATTGAAAAGGTTCTTAAATGCT‑CCGCTGAGTTATAGTTGTACAAATTTCATATTGCgtgggtgg                                                                                                                                  >  1:208163/1‑247 (MQ=255)                                                                                                               gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCtttt                                                                                                                           <  1:141499/251‑1 (MQ=255)                                                                                                                    tGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACca                                                                                                                       <  1:205778/250‑1 (MQ=255)                                                                                                                          cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGGTTCGCGGGGATGGGGGCGACCATGACAGTCTCCTTTTTCTTTATCGGGGGAAA‑CATTTTCATTTTTTTTTTT‑TGGTAATGCAATAGTTTATGAACTGAT‑CCGATGGGTTAATGTTTAACAAATCTCCTTTTTCGTGGGGGTCGCTTTTTCCACCGATGc                                                                                                                >  1:212092/1‑251 (MQ=255)                                                                                                                              aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTGCCGCGGTTCGCGGGGATAGGTGCGACCCTGACAGTCTCCTTTTTTTTTATCGTGGGAAT‑AATTTTCCTTTTTTTTGTT‑AGCTAATGAAAAAGATACTTAACTGAT‑CCGCTGAGTTAATGTTGTACAAATCTCATGTTGCGTGGTGGGCGCTTTTTCCACAGATGCggtt                                                                                                            >  1:58059/1‑248 (MQ=255)                                                                                                                                tccgccctagtcaTTTCAGGCCACGCGTTTAAGACTGACCCTGATCAGTCGCTGGTGCGGTGCCAACGGCGGG‑TTTTCTGCCGCGTGTCGCGTGAATGGGTTCGACCCTGACAGTCTACTTTTTCTGTATCGTTGAAAT‑CATTTTCATTTTTATTGTT‑AGTAAATGCAATAGTTATTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGTTGGTCGCTTTTACCACAGATGCGTTTa                                                                                                           <  2:422864/238‑1 (MQ=255)                                                                                                                                      ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCGCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGt                                                                                                    <  2:314066/251‑1 (MQ=255)                                                                                                                                            aTTTAAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGGCGCGTTTAGCGGGGATGGGGGGGAACATGACAGTCTCCTTTTTCTGTGTCCTGGGAAT‑TATTTTCATTTTTTTTGTT‑TACTAATTCAATAGTTACTGAACTGAT‑CCGACGAGTTAATGTTTAACAAATCTCTTTTTTCTGGGTGGGCGCTTTTTCCACCGATGCGGTTATGCCCGAATGGtt                                                                                              >  1:18833/1‑251 (MQ=255)                                                                                                                                                   gCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCCTTTTGCGGGGAAGGGGGCGACCATGAAAGTTTCCTTTTTTTTTTTCGTGGGAAT‑CATTTTTATTTTTATTTTTTGGTTAAT‑AAATTGTTACTGAACTGATTCCG‑TTAGTTAATGTTTAAAAAAACTCATGTTTGGTGGTGGGCGCTTTTACCAAAGATGCGGTTATGGCAGTATTGTTTGTTGa                                                                                        >  1:91793/1‑250 (MQ=255)                                                                                                                                                               aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCACCGGCGGG‑TTTTCTGCCGCCTTTTGCGGGGATGGGTGGGAACATGACAGGCTCCTTTTTCTGTATCGTGGGAAT‑TCTTTTCATTTTTATTGTT‑AGCTAATGCAATAATTACTGAACTGAT‑CCGCTGAGTTTATGTTGAACAAATATCATTTTGCGTGGGGGTCGCTTTTACCCCAGATGCGTTTATGCCAGTATGGTTTTTTGAATTTTTTTTAAAt                                                                           >  2:52777/1‑251 (MQ=255)                                                                                                                                                                                  gCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCATGACAATCTCCTTTTTCTTTATCGTGGAAAT‑CATTTTCATTTTTTTTGTT‑TGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTTATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCggg                                                        >  1:287105/1‑251 (MQ=255)                                                                                                                                                                                             ggCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCGAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGTTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGtgtg                                              <  1:541963/250‑1 (MQ=255)                                                                                                                                                                                              gCAACGGCGGG‑TTTGCTGCCGCGTTTCGCGTGAATGGGTGCTACCATGACAGTCTCCTTTTTCTGGATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTGTTACCACAGATGCGTGTATGCCAGTATGGTTTGTTGAAGTGTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGAGTGag                                            <  2:556112/251‑1 (MQ=255)                                                                                                                                                                                                     cggG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaa                                      >  2:484578/1‑250 (MQ=255)                                                                                                                                                                                                            ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGt                                                                                                                                                             <  1:486075/125‑1 (MQ=255)                                                                                                                                                                                                            ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGt                                                                                                                                                             >  2:486075/1‑125 (MQ=255)                                                                                                                                                                                                                gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaaca                             <  1:413978/249‑1 (MQ=255)                                                                                                                                                                                                                   gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgc                        <  1:303901/251‑1 (MQ=255)                                                                                                                                                                                                                            gggAATGGGTGCGACCATGACAGTCTACTTTTTCTGTATCGTGGCAAT‑CATTTTCATTTTTATTGTT‑AGCTATTGCAATAGTTAATGAACTGAT‑ACGATGCGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTTTTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCACAGTGGAACAACtgcgtgagccgt               >  2:206798/1‑245 (MQ=255)                                                                                                                                                                                                                               aaTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGTACTGAT‑CCGATGAGTTAATGTTGAACAAATATCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGAGCGCAGTTa             >  1:415578/1‑250 (MQ=255)                                                                                                                                                                                                                               aaTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCACAGTGGAACAACTGCGTGCGCAGTTa             >  2:13309/1‑250 (MQ=255)                                                                                                                                                                                                                                aTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa            >  1:142594/1‑250 (MQ=255)                                                                                                                                                                                                                                       cGACCATGACATGCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGAAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTGAAATGAACGt    <  2:465847/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt  <  2:554286/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAT‑CATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAT  >  NC_000913/4223393‑4223874

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