Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,459,225	T→C	21.9%	G415G (GGA→GGG)	fadI ←	beta‑ketoacyl‑CoA thiolase, anaerobic, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,459,225	0	T	C	21.9%	124.2 / 6.7	73	G415G (GGA→GGG)	fadI	beta‑ketoacyl‑CoA thiolase, anaerobic, subunit
Reads supporting (aligned to +/- strand):  ref base T (37/15);  new base C (2/15);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 2.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTCGCAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCGATCAATGTCAG  >  NC_000913/2458997‑2459461                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                cGCCTCGCAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCAcccccg                                                                                                                                                                                                                                        >  1:145231/1‑238 (MQ=255)       gCAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACCAAACC‑AAA‑TCCCCCCCCGCCC‑CGGc                                                                                                                                                                                                                                >  1:1081276/1‑239 (MQ=255)         ccttctTTTTTTTCACGGTGTTTGTTATATATTGCGCCCACCTGCGGGGCAATCCACCCCTTAAGGGTTTTCATTTCCCCCCCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGACCATTAAGGGTAAACGCTGATGCCATTTCCATTATTCCGCCACCAGAA‑CCATGGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAACCC‑AAA‑TCCACCGCCGCGA‑CGGCGa                                                                                                                                                                                                                              <  2:690489/236‑1 (MQ=255)           ctctTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGAGCCAAGCCCACCGGCAGCACAGGCGGAAACTAAAAC‑AAA‑ACCACCCCCCCGA‑CGGCGAAg                                                                                                                                                                                                                            >  2:408064/1‑239 (MQ=255)                tGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAACCCCCCGGCAGAACAGGCGGCACATAAAAC‑AAA‑TCCACCCCCGCCACCGGCG‑AGTTC‑ac                                                                                                                                                                                                                      >  2:925794/1‑238 (MQ=255)                     tCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCCCCAAGCCCACCGGCAGCACAGGCGGTAACCAAACA‑CAA‑ACCCCCCCCCCGA‑CGCCCAAGGTC‑CAGCaaa                                                                                                                                                                                                                 >  2:1170456/1‑238 (MQ=255)                      cACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAACCCCCCGGCCGCACAGGCGGGAAACAAACC‑AAA‑CACACCCCCCCGC‑CGGCGAAGTTA‑ATGCaaag                                                                                                                                                                                                                >  2:723094/1‑237 (MQ=255)                       cgggattgtgttAAAAAGGGCGCGCCCCGGCGAGGCAACCCCCGCCTTCGGGTTATTCATTTTCCCCCCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGt                                                                                                                                                                                                               <  1:203167/229‑1 (MQ=255)                            tttGGTTAATAAGGGCGCGCCCCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTgg                                                                                                                                                                                                           <  1:367819/237‑1 (MQ=255)                            gTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAAACAAAAC‑AAA‑CCCCCCCCCCCGC‑CCGCCAAAGTC‑AAGCaaattcaggg                                                                                                                                                                                                          >  2:1241311/1‑231 (MQ=255)                             ttGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCCCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGt                                                                                                                                                                                                         >  1:1011320/1‑239 (MQ=255)                                     taatGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCCCGC‑CGGCCAAGTCA‑ATCCAAAGGCTGGGTAatcatc                                                                                                                                                                                                  >  1:49478/1‑238 (MQ=255)                                       atGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGACGTTC‑CTGCAACGTCTGGGTAATCATcctc                                                                                                                                                                                               >  1:393055/1‑237 (MQ=255)                                          gcgcgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑CAA‑ACCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCcgcgc                                                                                                                                                                                             >  1:508672/1‑238 (MQ=255)                                          gcgcgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCCCCGCCCCGA‑CGGCGAAAGTC‑ATGCAAAGTCTGGGGAatcactgcccct                                                                                                                                                                                            >  1:1234658/1‑231 (MQ=255)                                          gcgcgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCCGGGGAATCATcccccc                                                                                                                                                                                             >  1:916769/1‑234 (MQ=255)                                          gcgcgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCCCCGCGA‑CGGCCAAGTTC‑ATGCAATGTCTGGGGAATCATCcgccc                                                                                                                                                                                             >  1:1234906/1‑236 (MQ=255)                                             cgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCg                                                                                                                                                                                         >  1:1191024/1‑239 (MQ=255)                                             cgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCCCCCCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGGAATCATcccgccgccg                                                                                                                                                                                         >  1:58263/1‑231 (MQ=255)                                             cgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCCCCCCGA‑CGGCCAAGTTC‑ATTCAATTCCTGGGCAAACATCCTCGccccc                                                                                                                                                                                         >  1:1032586/1‑238 (MQ=255)                                                cccggggAGGAAACCTCCCCCTTCAGGCTATTAGTTTTCCCACCCGGTACGTGGATGGTAAAACCGTCAAGGTTGTCCAGAAGAAAATTAAGGTCAAAGGCGAAGGCCTTATCCATTATCCCGCCTCCAGAA‑CCATGGCCGCCCCAACCCCACAGGCAGCACAGGCGGAAACTAAACC‑AAA‑CCCTCCGCCGCGA‑CGGAGAAGGTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTc                                                                                                                                                                                      <  2:1311082/233‑1 (MQ=255)                                                            aaCTCCGCCTTCGTGGTTTTATTTTTCCAACCGTGGACGTCGATGGTGCAACCGGAAATGTTGCCAAGGCGAACATTAAGGGTAAACGCTGATGCCATTTCCATTATCCCGCCTCAAAAA‑CCATTCCCGCGCCAAGCCCACGGGCAGCACAGGGGGTAACTAACCC‑AAA‑CCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGa                                                                                                                                                                          <  2:145231/239‑1 (MQ=255)                                                            aaCTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGGAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCCATGTCTGGGTAATCATCCGCGCGCCGGCCCCCGCGAGGGGa                                                                                                                                                                          >  1:1079242/1‑239 (MQ=255)                                                                    cTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCCCCGACA‑CGGCGAAGTTC‑ATGCAATGTCTTGGTAATCCTCCGCGCGCCGGTGGCCGCGAAGGGAAGGCCggg                                                                                                                                                                  >  1:136438/1‑237 (MQ=255)                                                                              ttatttttCCCCCCCGGTAGGTGGAGGGTGATAACGGAAATGTTGCCCAAACAAACATAAAGGGTAAACGCTGATTTCATTCCCATTTTCCCGCTCCCAAAC‑CCATCGCCGCGCCAGGCCCCCCGGCAGCACAGGGGGTAACTAAACC‑AAA‑CCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAAGGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGaa                                                                                                                                                        <  2:1234906/235‑1 (MQ=255)                                                                              ttCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑CAA‑TCCCCCGCCCCGC‑CCGCGAAGTTC‑ATGCAATGTCTGGGGAAACAACCACGCACCGGGCCCCCCGCAGGGATGCCCGAAAGCAAACGaa                                                                                                                                                        >  2:359196/1‑239 (MQ=255)                                                                              ttCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCCATTCCTGGGTAATCATCCGCCCCCCGGGCGCCCCGCAGGGGTGCCCGTAAGCAATCGa                                                                                                                                                         >  1:622849/1‑238 (MQ=255)                                                                              ttCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCCGGGTAATCATCCGCCCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGa                                                                                                                                                         >  1:906123/1‑238 (MQ=255)                                                                                      ctcACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGGAATAATCCGCGCGCCGGTCGCCCCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAg                                                                                                                                                >  1:611350/1‑239 (MQ=255)                                                                                               tACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGCCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CCGCGCAGTTC‑ATGCAATGGCTGGGGAATCATCCGCCCCCCGGTGGCCCGGAAGGGATGCCCCGAAGCAATCAAAACCCCacacagcgaaaa                                                                                                                                       >  2:469418/1‑228 (MQ=255)                                                                                                         gTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑CAA‑ACCCCCGCCCCCC‑CGGCCAAGTTCAATCCAATG‑CTGGGTCAACATCCCCCCGCCCGTCGCCGACAAGGGGAGGCCGTTGACAATCGGACCCCCAAGCCCGTT‑AAAATTGCTAtag                                                                                                                            >  2:56082/1‑237 (MQ=255)                                                                                                         gTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAAATTGCTATc                                                                                                                             >  1:316049/1‑238 (MQ=255)                                                                                                            aTAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGGTC‑ATGCAATGTCTGGGGAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCCCCAAGCACGGT‑AAATTTGCCATCGTcc                                                                                                                         >  1:635899/1‑239 (MQ=255)                                                                                                                 ggCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCCCCCCGA‑CCGCGCAGTTC‑ATGCAATGTCTGGGTAATCACACCCGCGCCGGCCCCCCCGAAGGGGTGCCCCGAACGAATCCACCCCCCCACCACCTTAAAATTTGCAAtgcccccct                                                                                                                      >  2:1113788/1‑230 (MQ=255)                                                                                                                     attttggcaagaagaACATTAAGGGTAACCGCTAATGTCATTCCCATTATCCCCCCTCCAGAA‑CCATGGCCGCGCCAAGCCCCCCGGCAGCACAGGGGGAAACTAACCC‑AAA‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCa                                                                                                                <  1:1193263/230‑1 (MQ=255)                                                                                                                        ttgtCCAGACGAACATTAAGGGTAACCGCTGATGTCATTTCCATTATCCCGCCTCCAGAA‑CCATTGCCGCGCCAACCCCACCGGCAGCACGGGCGGTAACTAACCC‑AAA‑CCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGtt                                                                                                                                          <  1:1149763/211‑1 (MQ=255)                                                                                                                           tCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑CAATTTGCTATCGTTCCCTTCGCCAGTGGCa                                                                                                          >  1:198879/1‑239 (MQ=255)                                                                                                                           cccAGACGAACATTAGGGGTAAACGGTGATGCTATTTCCATTATTCCGCCCCCAGAC‑CCATTGCCGGGCCAAGCCCCCCGGCAGCACAGGCGGTAACTAAACA‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCa                                                                                                          <  1:1113788/238‑1 (MQ=255)                                                                                                                               gcgtaaCATTAAGGGTAAGCGCTGATGTCATTCCAATTATTCCGCCTCCAAAA‑CCATTGCCGCCCAAAGCCCACCGGCAGCCCAGGCGGTAACTAACCC‑AAA‑CCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCAtgtg                                                                                                      <  1:925794/235‑1 (MQ=255)                                                                                                                                 ccacatttAAGGTAAAACGCTGATGTCATTCCCATTATCCCGCCCCCAGAA‑CCATGGGCGCCCCAAGCCCCCCGGCAGCACGGGCGGTAAATAAACC‑AAA‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCa                                                                                                    <  1:469418/233‑1 (MQ=255)                                                                                                                                     cATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAACCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGCAGGGAAGCCCGTAAGCAATACAACCCCCCACCACGTT‑CACTTTGCTATCCTCCACCTCCCCAGAGACATGTGccccc                                                                                                 >  2:126085/1‑234 (MQ=255)                                                                                                                                           gggTAACCGCTGATGTCATTCCCATTATCCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCCCCGGCAGCAAAGGCGGTAACTAAACC‑AAA‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACg                                                                                                  <  1:545805/231‑1 (MQ=255)                                                                                                                                            ggTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGGAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCCCCAGTGCCCTTGGCCCGcccccgg                                                                                           >  1:1217151/1‑234 (MQ=255)                                                                                                                                                aaCGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAAATTGCTATCGTCCATTTCCCCAGTGGCATTGGCACCCCCCCGTGGTTc                                                                                      >  1:818134/1‑238 (MQ=255)                                                                                                                                                aaCGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGGTCCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCACGCCCGTT‑AAATTTTCTCTCCTCCCCTTCCCCAGGGGCATGTGAACCCCCCAACGCTTc                                                                                      >  2:1263787/1‑238 (MQ=255)                                                                                                                                                      gATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCCCGCCCCACTGGTTCCCGCGCa                                                                               >  1:1050240/1‑239 (MQ=255)                                                                                                                                                      gATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGAATGTGGACGCCCCCGAGCTTCAcacgcc                                                                               >  1:14074/1‑234 (MQ=255)                                                                                                                                                      gATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAAACGCCAAGCACGTT‑AAATTTGCTATCGTCCCATTCGCCAGTGGCCTTTGCACGCCCCAAGGCTTCCCGAGcc                                                                               >  1:313803/1‑238 (MQ=255)                                                                                                                                                            aTTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCATTGGCATGGGCAAGCCCCCGGGCTTCAccggacaaac                                                                           >  1:190322/1‑228 (MQ=255)                                                                                                                                                            aTTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGGTCACGCCCCAGTGCCTCACGAGCCAAACGt                                                                         >  1:1066584/1‑239 (MQ=255)                                                                                                                                                             tttCCATTATTCCGCCTCCAGAA‑CCATTGCCGCCCCAAGCCCCCCGGCAGCACGGGCGTTAACTAAACC‑AAA‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGtt                                                                        <  1:723094/239‑1 (MQ=255)                                                                                                                                                                  attagcccGCCCCCAGAC‑CCATTGCCGCGCCAACCCCGCCGGCAGCACGGGCGTAAACAAAACC‑AAA‑CCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTa                                                                   <  1:126085/233‑1 (MQ=255)                                                                                                                                                                     tttCCCCCCCCACAACCCATCGCCGCCCCGACCCCCCCGGCCGCCCAGG‑GGTAAAAAAACC‑AAACCCCCCCGCGGGGA‑CGGCGAAGTTC‑ATGCAATGTGTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTAcc                                                                 <  2:1011320/238‑1 (MQ=255)                                                                                                                                                                           cccccAAAA‑CCATTGCCGGCCCAAGCCCCCCGGCAGCACGGGCGTTAACAAAACC‑AAA‑CCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAAc                                                          <  1:1170456/236‑1 (MQ=255)                                                                                                                                                                           ccTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACt                                                                                                                      >  1:728825/1‑179 (MQ=255)                                                                                                                                                                                        tgcccGCCCCAAGCCCCCCGGCAGCCGAGGTGTTAACAAACCC‑AAC‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGcc                                              <  1:1164112/236‑1 (MQ=255)                                                                                                                                                                                                ccagccccACCGGCAGCACAGGGGGTACCAAAACC‑AAA‑CCCGCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTg                                      <  1:343340/234‑1 (MQ=255)                                                                                                                                                                                                 caccccccccGGCTACACGGGGGATAACACAACC‑AAC‑CCCCCCGCCGCGG‑CGGCGAAGTTC‑CTGCAATGCTTGGGTAACCCACGCCGCGCGGGTCGCCACGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGa                                     <  2:611350/235‑1 (MQ=255)                                                                                                                                                                                                      ccACCGGCAGCACGGGCGTTACCAAACCC‑AAC‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGc                                   <  1:359196/236‑1 (MQ=255)                                                                                                                                                                                                         ccGGCAACAAGGGCGGAAACTAACCC‑AAA‑CCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCaaa                             <  1:1263787/239‑1 (MQ=255)                                                                                                                                                                                                             cagccaaGGCGGTAACTAACCC‑AAC‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGc                          <  1:1271102/232‑1 (MQ=255)                                                                                                                                                                                                             caccacgGGCGGAAACTAAACC‑AAA‑TCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGc                          <  1:342092/235‑1 (MQ=255)                                                                                                                                                                                                                        tgtaCTACACC‑AAC‑CCCCCCCCCGCGG‑GGGCGAAGTGC‑ACGCATGGTCTGGATAACCCTCCGCGCGCGGGTCGACGTGAAGGGATGCCCTTAAGAGAGCGCACCGCCAAGCACGTT‑AAATTTGCTATCGACCACTTCGCCAGGGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCatat               <  2:704015/235‑1 (MQ=255)                                                                                                                                                                                                                        gTAACTAAACC‑AAA‑TCCCCCGCCGCGA‑CGGCGAAGTTA‑ACGCAATGTCTGGGAAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCatat               <  1:79202/238‑1 (MQ=255)                                                                                                                                                                                                                          aaCTAACCC‑AAA‑TCCCCCGCCGGGC‑AGCAGAAGTTC‑ATGCAAGGCCTGGGTAAACATCCGCGCGCCGGCCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCATTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCg             <  1:613212/238‑1 (MQ=255)                                                                                                                                                                                                                           acaaaccc‑AAA‑CCCACCGCCGGGA‑CGGCAAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCGAt           <  1:817890/233‑1 (MQ=255)                                                                                                                                                                                                                           acaaaacc‑aaa‑ccCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGCCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCGAt           <  1:490493/236‑1 (MQ=255)                                                                                                                                                                                                                           aacacaCC‑ACA‑CCCCCCCCCGCGG‑AGGGGAAGTTC‑ATGCCATGGCTGGGTAATCCTCCGCGGGCCGGTCGCCGGGAAGGGATGCCCGTAACCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCGAt           <  2:1069821/234‑1 (MQ=255)                                                                                                                                                                                                                            ccaaacc‑aaaccCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCGAt           <  1:722406/237‑1 (MQ=255)                                                                                                                                                                                                                            acaaaCCAAAA‑CCCCCCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCGAt           <  1:1241311/237‑1 (MQ=255)                                                                                                                                                                                                                               aaCC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGa                                                                                                                                                         <  2:1021423/94‑1 (MQ=255)                                                                                                                                                                                                                               aaCC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGa                                                                                                                                                         >  1:1021423/1‑94 (MQ=255)                                                                                                                                                                                                                               aaCC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCACCGCCTGAGCTGCAAAAGCTTCGTGCATATCGATCAAt       >  1:807608/1‑239 (MQ=255)                                                                                                                                                                                                                               aaCC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCCATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCACAACCTTCACCACCCACCAACTCAAACTTCCCCGGCGCCCGAACTCCAAAGACCTCGGGCAAATcaaaaaaa       >  2:631572/1‑232 (MQ=255)                                                                                                                                                                                                                                     aa‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACCCCCCAGTGCTTCACGGGCAAAACGGTCACTACCCAGCAACTGAATATTCGCCCGCGTCTGAGCGGCAAAGGCTTCGTGCCTATCGATCAATTTCAg  >  1:1069294/1‑239 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                CGCCTCGCAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAA‑CCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACC‑AAA‑TCCACCGCCGCGA‑CGGCGAAGTTC‑ATGCAATGTCTGGGTAATCATCCGCGCGCCGGTCGCCGCGAAGGGATGCCCGTAAGCAATCGAACCGCCAAGCACGTT‑AAATTTGCTATCGTCCACTTCGCCAGTGGCATGTGCACGCCCCAGTGCTTCACGAGCAAAACGTTCACTACCCAGCAACTGAATATTCGCCAGCGTCTGAGCTGCAAAGGCTTCGTGCATATCGATCAATGTCAG  >  NC_000913/2458997‑2459461

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