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%	153.3 / NA	43	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 (9/34);  total (9/34)

CAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCA‑TCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCT  >  NC_000913/4223392‑4223877                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   cAGTTCGCCAACCTGACTCACGAAACCCTGCTGTTGCA‑TCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCTCCAGTTCCGTTAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTTTCCAGTCGCTGTTGCGGTGGCAACGGCGTGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCctttt                                                                                                                                                                                                                                                <  1:255199/251‑1 (MQ=255)      cGCCAACCTGACGCCCGAAACCCTGCTGTTGCA‑TCGTGGTTAGCAGGCGGTGGGTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTGTT‑CTGCCGCTTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTGTTCTGt                                                                                                                                                                                                                                           <  2:390041/251‑1 (MQ=255)          aaCCTGACGCACGAAACCCTGCTGGTGCA‑ACGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGTATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑g                                                                                                                                                                                                                                      <  2:48446/251‑1 (MQ=255)            ccTGACGCACGAAACCCTGCTGTTGCA‑TCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTg                                                                                                                                                                                                                                    <  1:238239/251‑1 (MQ=255)                  ccccGAAACCCTGCTGTTGCA‑TCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAAtc                                                                                                                                                                                                                              <  2:364004/248‑1 (MQ=255)                      gaagcccTGCTGTT‑CACTCGTGGTTAGCAGTCGGTGGGTCTTGGAATTGGGTAACCCTGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATTGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAAtcattt                                                                                                                                                                                                                          <  2:282505/247‑1 (MQ=255)                        aaCCCTGCTGTTGCA‑TCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAAtcattttc                                                                                                                                                                                                                        <  1:358887/251‑1 (MQ=255)                               cgtgtGCA‑TCGTGGTTAGCCGGCGGTGGGTCGTGGAATTGGGGAACCCGGCTTGTTGCGCCAGTTCCGTGAGGGCCCCACCGCCATTGGATTCGGCAATCCACCCCAGTAATTTCAGGCCCCGCGTTAAAGGCTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTTCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTa                                                                                                                                                                                                                 <  2:165678/247‑1 (MQ=255)                                               tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAGTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGc                                                                                                                                                                                                  <  1:242929/250‑1 (MQ=255)                                                  cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATa                                                                                                                                                                                              <  1:125177/251‑1 (MQ=255)                                                        gtggttcgtgtgaTTGGGTAACCCGGCTTGTTGCGCGAGTTCCGTGAGGGCTACTCTGCCATTTGATTCGGCAATCCACGCCATTAATTTCATGCCACGCGTTAAAGACTGAACCTGTCCAGTCTCTGGTGCGGTTGCAACGGCGGTTTTT‑CTGCCGCGTTTCGCGGTAATGGGTGCGATCATGAC‑AGTCTCCTTTTTCTTTTTC‑GGGGAAATCATTTTCATTTTTTTGGTTAGCTAA‑TGCAATAGTTact                                                                                                                                                                                        <  2:162282/239‑1 (MQ=255)                                                              cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTactgaact                                                                                                                                                                                   <  1:136984/250‑1 (MQ=255)                                                                            aCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGTGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCTACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATGTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGGACTGATCCGATGAGTTaa                                                                                                                                                                    <  1:402689/251‑1 (MQ=255)                                                                                          gcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTTCCGCGGTTCGCGGGGATGGGGGCGACCATGAC‑AGTCTCCTTTTTCTGTTTC‑GTGGGAATCATTTTTATTTTTTTTTTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGGGTTAATGTTGAACAAAtct                                                                                                                                                      >  1:208073/1‑251 (MQ=255)                                                                                            cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTT‑CTGCCGCGGTTTGCGGGGATGGGTGGGACCATGAC‑CGTCTCCCTTTTTTGTTTT‑GTGGGAATACTTTTTATTTTTATTTTTTGCTAA‑TGCAATAATTACTGAAATGATCCGAGGAGTTAATGTTGAACAAAtctc                                                                                                                                                     >  1:185989/1‑250 (MQ=255)                                                                                            cAGTTCCGTGAGTGCCAAACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACTCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCa                                                                                                                                                    <  2:401495/251‑1 (MQ=255)                                                                                             aGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑TTGCCCCGTTTCGCGGGGATGGGGGCGGCCCTTAC‑CATCTCCCTTTTCTGTATC‑GTGGGAATCATTTTTATTTTTTTTTTTTGCTTA‑TTTCATAGTTTCTGAAATGGTACGATGAGTTTAAGTTTGACAAAtctcct                                                                                                                                                   >  1:344911/1‑249 (MQ=255)                                                                                              gTTCCGTTAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAAGAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATg                                                                                                                                                  <  2:128042/251‑1 (MQ=255)                                                                                                         gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCCGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtgg‑tg                                                                                                                                       <  1:203688/250‑1 (MQ=255)                                                                                                                 gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCTCGCGTTAAAGACCGAACCTCTCCAGTCTCTGCTTCTGTGTCAACGGCGTGTTTT‑CTTCCTTGTTTCGCGTGAATGGGTGCGGCCATTAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCTATACTTACTGAACTGATCCTATGAATTAATTTTGAACAAAACTCATATTGCatggatgtctctttt                                                                                                                               <  2:168845/251‑12 (MQ=255)                                                                                                                  ccATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTCT‑CCGCGTTTCGCGGGAATAGGTGCGACCATGGC‑ACTCTCCCTTTT‑TTTATTTGTGTAAATCATTTTTTTTTTTTTTTTTAGCTTA‑TGCAATAGGTTCTGAACTTCTCCCCTGAGTTCATCATGAACAAAACACACGTTGCGTGG‑TGGTCGCttttt                                                                                                                             >  1:181758/1‑250 (MQ=255)                                                                                                                             ggCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTCCCACAGATGCg                                                                                                                  <  2:243004/251‑1 (MQ=255)                                                                                                                               cAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTT‑CCGCCGCGTTTTGCGGGGATGGGGGCGACCATGAC‑AGGCTCCTTTTTTTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTTGATAA‑TGCAATAGTTACTGAAATGACCCGATGAGTTAATGTTTAACAAATATCATGTTGTGTGG‑GGGGCGCtttttccccaaaggcgtt                                                                                                                >  1:332473/1‑237 (MQ=255)                                                                                                                                 aTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTT‑TTCCCCAGTTTCGCGGGGATGGGGGCGACCATGACCAG‑CTCCTTTTTTTGTATC‑ATGGGAATCATTTTTCTTTTTTTTTTTTGCTAA‑TGCAATAGGTTCTGAACTGATACGCTGAGGTAATATTGAACAAAACTCATGTTGCGTGG‑GGGTGGCTTTTACCACAAATGGGtttt                                                                                                              >  1:248020/1‑250 (MQ=255)                                                                                                                                   ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCa                                                                                                                                                                                                 <  1:80107/167‑1 (MQ=255)                                                                                                                                       cccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CGGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTGAGCTAA‑TGCAATAGTTACTGAACTGATCCGAGGAGTTAATGTTGAACAAATCTCATGTTTCGTGG‑TGGTCGCGTTTACCACAGATGCGTTTATGCCAg                                                                                                        <  1:184870/250‑1 (MQ=255)                                                                                                                                                  cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgtt                                                                                             <  1:322692/251‑1 (MQ=255)                                                                                                                                                  cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTGGCGGGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGTTttgtt                                                                                             <  1:262814/251‑1 (MQ=255)                                                                                                                                                          gcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGATTT‑CTGCCGGGGTTCGTGGGGATGGGGGGGACCCAGAC‑AGTCTCCTTTTTTTGTTTC‑GTGGAAAACATTTTTATTTTTATTGTTACATAAATG‑AATAATTACTGAACCGGTACGGTGAGTTAATATTGAAAAAATATCATGTTGCGTGGGTGG‑CGCTTTTTACAAAGGACCGCTTATGCCAGTATGGTTTGTTGAAttttt                                                                                     >  2:121390/1‑251 (MQ=255)                                                                                                                                                            gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttat                                                                                   <  2:355969/251‑1 (MQ=255)                                                                                                                                                                      gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTgg                                                                          <  1:5968/250‑1 (MQ=255)                                                                                                                                                                       aaCCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGAGGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGt                                                                        <  2:301131/251‑1 (MQ=255)                                                                                                                                                                                   cGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCg                                                             <  2:344911/250‑1 (MQ=255)                                                                                                                                                                                           cGGTGGCAACGGCGTGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTTAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGGTGAATTTTTATTAAATCTGTGTTGAGCGTGTCGGGAGCAAgt                                                    <  1:284725/251‑1 (MQ=255)                                                                                                                                                                                               ggCAACGGCGGGGTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCGCCTTTTTCTGTATC‑GTGGAAAGCATTTTCATTGTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAAGCTCATGTTGCGTGG‑TTGTCGCTTTTGCCACATAGGCGTTTATGCGAGTATGGGTTGTTGAAGTTTTATGAAATCTGGGTTGAGCGGGTCGGGAGCAATTGTGa                                                <  2:383614/251‑1 (MQ=255)                                                                                                                                                                                                         gggtgt‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                      <  2:130555/246‑1 (MQ=255)                                                                                                                                                                                                            ttt‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg                                    <  2:208073/250‑1 (MQ=255)                                                                                                                                                                                                                      cgTTTCGCGGGAATGGGAGAGCCCATGAC‑CGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑AGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcg                          >  2:178772/1‑251 (MQ=255)                                                                                                                                                                                                                            gcgGGAATGGGTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAAGGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                    <  2:161107/251‑1 (MQ=255)                                                                                                                                                                                                                                    gggTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt             <  1:385067/250‑1 (MQ=255)                                                                                                                                                                                                                                      gTGCGACCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa          <  1:191046/251‑1 (MQ=255)                                                                                                                                                                                                                                         cGACCATGAC‑CGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTACATGAACGt       >  2:250644/1‑251 (MQ=255)                                                                                                                                                                                                                                           aCCATGAC‑AGTCTCCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAAGCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCGGGGTTGAGCGTGTCGGGAGCAAGTGTGGGCAGCAAAGTGGAACAACTGCGTGCGCAGGTAAATGAACGTAt     <  1:311188/251‑1 (MQ=255)                                                                                                                                                                                                                                              aTGAC‑AGTCTTCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCt  <  2:299931/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   CAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCA‑TCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAA‑TGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGG‑TGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCT  >  NC_000913/4223392‑4223877

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