Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,336	A→G	E201E (GAA→GAG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336	0	A	G	94.5%	206.6 / 0.6	74	E201E (GAA→GAG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (32/37);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGC  >  NC_000913/4094‑4583                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        atTGCGGGCGATAAGCCAGTGACCATTCTGCCCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAgg                                                                                                                                                                                                                                                 >  1:368800/1‑250 (MQ=255)    gCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGgc                                                                                                                                                                                                                                              <  1:98832/251‑1 (MQ=255)       ggcgATAAGCCAGTGACCATTCTGACCGCAACCTCAGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGGGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGGGgcgcc                                                                                                                                                                                                                                           >  2:286116/4‑250 (MQ=255)          gATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCtggg                                                                                                                                                                                                                                        >  2:325572/1‑248 (MQ=255)                     gacaaTTCTGACCGCGACCTCCGGGGATAACGGAGTGGCAGTGGCTTATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCGCTAGCCACGAGGCAAAAACAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTca                                                                                                                                                                                                                              <  2:547929/246‑2 (MQ=255)                     gACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTcag                                                                                                                                                                                                                             <  1:19743/251‑3 (MQ=255)                     gACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTca                                                                                                                                                                                                                              <  1:491467/250‑2 (MQ=255)                         aTTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTaa                                                                                                                                                                                                                         <  1:259999/251‑1 (MQ=255)                              gACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTTAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGa                                                                                                                                                                                                                    <  1:467670/251‑1 (MQ=255)                                 cgcgACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTa                                                                                                                                                                                                                 <  1:190235/251‑1 (MQ=255)                                           ggTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCg                                                                                                                                                                                                       >  1:316306/1‑251 (MQ=255)                                            gtgcgaCCGGAGAGACAGTGGCTAATGCTTTCTACGGTTTACCGAATGGGAAAGAGGTTATCCTCTATACACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGTCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCgg                                                                                                                                                                                                      <  2:120943/246‑2 (MQ=255)                                            gTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGACATATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCTGTTGATGATGA‑AGGGCTGACAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGt                                                                                                                                                                                                      >  2:364255/1‑251 (MQ=255)                                                    ggAGCCGCAGGGGCTCATGCTTTCTACCGTTTACCGAAGGTGAAAGTGGTTATCCTGTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGGACATTGGGCGGCAATATCGAAACTTTGGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGTTGAGGCAGGCGTTTGATGATGA‑AGAGCTTAAAGTGTCGCTGGGGTTAATCACAGCTAACTCGCTTAACATGAGTCGTTTGCTGgc                                                                                                                                                                                              <  2:655131/251‑1 (MQ=255)                                                        cGGCAGAGGCTCATGCTATCTACGGTTTACCGAATGTGCAAGTGGTTATACTCTATCCACGAGGCAAAATCAGTCCCCTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAAATGTTGCCATAGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCAGAAAGGGGCGCTGTGGTTAAACTCAACTAACTCGATTAACATCAGCTGTTTGCTGgcgac                                                                                                                                                                                           >  2:560991/1‑248 (MQ=255)                                                         ggCAGTGGCTCATTATTTCTCCGGATTACCGAATGTGAAAGGGGTTATCCTCTATCCAAGCGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCAAAACTGTTGCAATAGACGGCGATTTAGATGACTGTCAGGCGCTGGAGAAGCAGGAGTTTAATGATGA‑CGAGCTGAAAGGGGCGCTGGGGGTAAAATCCGCTAACTCGTTTAACAACAGCAGTTTGCGGGCGCAgg                                                                                                                                                                                         >  2:468910/1‑250 (MQ=255)                                                            aGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAttt                                                                                                                                                                                      <  1:325572/251‑1 (MQ=255)                                                                            tACGGTTTACCGCATGTGAACGAGGTTATCCTCTATCCACGCGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGCTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGCGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGCTTTGCTCCTACTTTGAAGc                                                                                                                                                                      >  2:503454/1‑251 (MQ=255)                                                                             aCGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAg                                                                                                                                                                       <  1:262686/249‑1 (MQ=255)                                                                             aCGGTTTACAAAAAGTGAAAGTGTTTAGCCTGTATCCACGAGGCAAATTCAATCCACTGCAAGCAATACTGTTCTGTACATTGGGCGGCAATATCGAACGTGATGCCACCGACGTCGATTTCGATGCCTGTCAGGCGCTGGTGGTGCAGGCGTTTGGTGATGA‑AGAACGAAAAGTGGCGCTAGGGTTCAACTCGGATAACTCGATTAAAATCAGCCTATCGCTGGCGCAGATGTGCTAAGACTTTGAAGcc                                                                                                                                                                     >  2:268457/1‑250 (MQ=255)                                                                                gTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                  >  1:4511/1‑251 (MQ=255)                                                                                 tttACCGAATGTGAAAGTGGTTATCCGCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                  <  2:456319/250‑1 (MQ=255)                                                                                       gAATGTGAAAGTGGTAATACTCTATCCACGAGGCAAAATCAGTCCAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGCAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGAAGGCGTTTGCTGATGA‑AGACCTGAAAGAGGCGCTAGGGTTAAACACGGCTAAATCGATTAAAAACAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTgcgc                                                                                                                                                              >  2:593325/1‑248 (MQ=255)                                                                                          tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCCGCTAACTCGATTACCATCAGCCGTTTGCTGGCGCGGATTTGCTACTACTTTGAAGcag                                                                                                                                                                    >  2:418725/1‑237 (MQ=255)                                                                                          tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGcag                                                                                                                                                                    <  1:418725/239‑3 (MQ=255)                                                                                                                ccACGAGGCAAAATCATGCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCt                                                                                                                                  <  1:114863/251‑1 (MQ=255)                                                                                                                        cAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAACGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtct                                                                                                                          >  2:279846/1‑251 (MQ=255)                                                                                                                        cAAAATCAGGCCCCTGAAAGAAAAACTGGTCTGTACATTGGGAGGCAATATCGAAACTGTTGCCACCGACGGCGATGTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGTCACTAGGGTTAAAATCGGCTAACTCGATGAACATCAGCCGTTTGCTGGCGCAGATTTGCTCCTACTTGGAAGCTGTTGCGCGGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGtct                                                                                                                          <  2:51029/251‑1 (MQ=255)                                                                                                                            aTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                      <  1:537250/251‑1 (MQ=255)                                                                                                                               aGTCCACTGCAAGAAAAACTGTTCTGGACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGACAGTGGCGCTGGGGTTAAACTCAGCTCACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCATTTGCGCAGCGGCCGCAGAAAGCGCGCACCCACCTGGTTGtcgcgtgtcc                                                                                                                   >  2:581461/1‑243 (MQ=255)                                                                                                                                        cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                           >  1:539105/1‑250 (MQ=255)                                                                                                                                        cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaaa                                                                                                          <  1:286116/251‑1 (MQ=255)                                                                                                                                         aaGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAAATCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAAc                                                                                                         >  1:379782/1‑251 (MQ=255)                                                                                                                                          aGAACAACAGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCAACGGCGCTTTCGATGTCTGACAAGCCCTGGTGAAGCAGACGCTTGATGATGC‑GGAGCTGAAAGGGGCGGTGGGGGTAAACTCAGCTAACTCGATTACCATCAGCGGTGTGCTGGCGCAGGTTGGATACTACTTTGAAGCAGTTGCGTCGCTCCCGCAGGAAGCGCACAACCAGCTGGATGTACCGGTGcccacaggaaact                                                                                                        >  2:13155/1‑240 (MQ=255)                                                                                                                                          aGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACt                                                                                                        >  2:536611/1‑251 (MQ=255)                                                                                                                                                 cTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGCCGGCGATTTCGATGCCTGTAAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCgg                                                                                                 >  2:34376/1‑250 (MQ=255)                                                                                                                                                         tACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGTGTTGAcgg                                                                                         >  2:372779/1‑251 (MQ=255)                                                                                                                                                          aCATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                           >  1:137070/1‑216 (MQ=255)                                                                                                                                                          aCATTGGGCGGCAATATCGAAACTGTGGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                           <  2:137070/216‑1 (MQ=255)                                                                                                                                                           cATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGAcggcg                                                                                       >  1:397026/1‑251 (MQ=255)                                                                                                                                                              tGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                      >  1:487639/1‑217 (MQ=255)                                                                                                                                                               gggCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGt                                                                                    <  1:313348/250‑1 (MQ=255)                                                                                                                                                               gggCGGCAATAACGAACCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGTAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTc                                                                                   <  2:122206/251‑1 (MQ=255)                                                                                                                                                                   ggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGACTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCACACGGAAACTTCGGCGATTTGCCGGCgggggtgtt                                                                               >  2:319302/1‑245 (MQ=255)                                                                                                                                                                     cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCCGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGCAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGt                                                                                                                                                                   >  1:645176/1‑165 (MQ=255)                                                                                                                                                                     cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATAAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGtacgctgg                                                                             >  1:380436/1‑244 (MQ=255)                                                                                                                                                                        tatCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTggggg                                                                          >  1:393309/1‑248 (MQ=255)                                                                                                                                                                                  gTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCg                                                                 <  1:104744/250‑1 (MQ=255)                                                                                                                                                                                     gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCt                                                             <  1:536611/251‑1 (MQ=255)                                                                                                                                                                                      ccatggccggagagttcggggCCTGTCAGGCGCTGTTGAAGCGTGCGTTTGAGGA‑GAGAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGAGTAACTTCAGCCGTTGGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCTCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACTCGGTCTg                                                            <  2:588811/231‑1 (MQ=255)                                                                                                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGGAGGAGTTTGATGATGA‑AGAGCTGAAAGTGGCGATGGGGTTAAACTCAGCTAACTCGATTAAAATCAGCCGTTTGCTGGCGCAGAAGTGCTACTACTTTAAAGAAGTTGCGAAGATGCAGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAACATTCGGCGATTTGACGGCGGGActgtgggcgggagggacaggggtg                                                            >  2:659677/1‑230 (MQ=255)                                                                                                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCCGGtgggctgggggaggaccgtggatgg                                                            >  2:514265/1‑227 (MQ=255)                                                                                                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                              >  1:452215/1‑201 (MQ=255)                                                                                                                                                                                      ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                              <  2:452215/201‑1 (MQ=255)                                                                                                                                                                                        aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGcc                                                          <  1:245566/251‑1 (MQ=255)                                                                                                                                                                                        aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCGGCTGGCggcgacaacccgcgcggc                                                          >  2:312113/1‑234 (MQ=255)                                                                                                                                                                                           gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                  >  1:146957/1‑144 (MQ=255)                                                                                                                                                                                           gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                  <  2:146957/144‑1 (MQ=255)                                                                                                                                                                                            aCGGCGATTTCGATTCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTg                                                      <  1:330595/251‑1 (MQ=255)                                                                                                                                                                                             gggCGATTTGAATGCCTGTCAGGCGTGGGTTAAGCAGGCGTTGGAGGAGGA‑AGAGCGTAAAGTGGCGCTGGGTTTAAACTCAGCTAACGCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGa                                                     <  2:314301/250‑1 (MQ=255)                                                                                                                                                                                             cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                     >  1:401063/1‑187 (MQ=255)                                                                                                                                                                                             cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTGGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                     <  2:401063/187‑1 (MQ=255)                                                                                                                                                                                              ggCGATTTCGATGCCTGTCAGGCGCTGGCGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGTGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAAACAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGaa                                                    <  1:259989/251‑1 (MQ=255)                                                                                                                                                                                                cGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGATGCTGGCggaggtccaccgtctgccggggaga                                                   >  1:232509/1‑226 (MQ=255)                                                                                                                                                                                                  aTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGt                                                <  1:617040/251‑1 (MQ=255)                                                                                                                                                                                                     tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGttt                                              <  1:431307/250‑1 (MQ=255)                                                                                                                                                                                                           ccTGGCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACGCAGCTAAGCCGATGAACATCAGCCGGTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCTATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATtgct                                       <  2:397026/251‑1 (MQ=255)                                                                                                                                                                                                                cAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGAc                                  <  2:316306/251‑1 (MQ=255)                                                                                                                                                                                                                          ggagaccggCGTTTGATGATGA‑AGAGCTGAAAGGGGCGCTGGGGTAAAACTAAGCTAACTCGATTAACAACAGCGGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAGCGTGAAc                        <  2:4511/244‑1 (MQ=255)                                                                                                                                                                                                                                aGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGTTTAAACTCAGCTAACTCGCTTAACATCGGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTa                                            <  2:542452/225‑1 (MQ=255)                                                                                                                                                                                                                                 ggCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCg                 <  1:268960/251‑1 (MQ=255)                                                                                                                                                                                                                                       tgatgatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTTTTTCTGCGACCAACGTGAACGATACCGTGcccc           >  1:193019/1‑249 (MQ=255)                                                                                                                                                                                                                                          ggatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                                                  <  2:19485/176‑1 (MQ=255)                                                                                                                                                                                                                                           gatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACCCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCGGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGttt       <  2:193019/251‑1 (MQ=255)                                                                                                                                                                                                                                                 agaaCTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGc  <  1:312113/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        ATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGC  >  NC_000913/4094‑4583

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