Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 789,580 C→A E134* (GAA→TAA)  galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,5800CA98.7% 252.2 / NA 76E134* (GAA→TAA) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/41);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTG  >  NC_000913/789338‑789818
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
tcatttACTTTTGTTGATGACGCCAGCCACTCCTTTGGGCATGGAGACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTTCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCATTACGGTTCCGTCCGCGTCTTACAGTGaa                                                                                                                                                                                                                                          <  2:124260/246‑1 (MQ=255)
gAAGGTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCGTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGAACCATGATCCCGCAGTTACAGCCTACAACCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAAGACGGTTCAGACCGCGACGTTCAGTGAAg                                                                                                                                                                                                                                         <  2:383477/251‑1 (MQ=255)
  aGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCg                                                                                                                                                                                                                                       >  2:328869/1‑251 (MQ=255)
           tgatgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACtt                                                                                                                                                                                                                              >  1:150886/1‑251 (MQ=255)
             atgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCCTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTa                                                                                                                                                                                                                             >  1:394424/1‑250 (MQ=255)
              tgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGCACTTAAc                                                                                                                                                                                                                           >  1:136518/1‑251 (MQ=255)
                      acAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCAcc                                                                                                                                                                                                                   <  1:259503/251‑1 (MQ=255)
                         gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑g                                                                                                                                                                                                               >  2:203453/1‑251 (MQ=255)
                         gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑g                                                                                                                                                                                                               >  1:173873/1‑251 (MQ=255)
                         gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGAGCCCGCAGTTACAGCCTTCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCCGCTGCAATACGGTTCCGACCGCGACTTACAGTTAAGCGGAAGAACTTAACCCGGCACCCT‑g                                                                                                                                                                                                               <  2:458368/251‑1 (MQ=255)
                             cacCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGc                                                                                                                                                                                                           >  1:300016/1‑251 (MQ=255)
                              acCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGca                                                                                                                                                                                                          <  1:49149/251‑1 (MQ=255)
                                cTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCACGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGcac                                                                                                                                                                                                         >  1:358364/1‑250 (MQ=255)
                                    gggCATGGAAGCTGCTTTGGTCCCCAGTGAGCTGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAATCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAAGACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCTCCCT‑GCGGCACATTGc                                                                                                                                                                                                    <  2:217573/251‑1 (MQ=255)
                                          ggAAACTGCTGTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTGTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGa                                                                                                                                                                                              <  2:320431/251‑1 (MQ=255)
                                                     tGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTc                                                                                                                                                                                    <  1:10078/250‑1 (MQ=255)
                                                      ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTccc                                                                                                                                                                                  >  1:131355/1‑250 (MQ=255)
                                                      ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTcc                                                                                                                                                                                   <  2:491944/250‑1 (MQ=255)
                                                       gTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCAc                                                                                                                                                                                 <  2:150886/251‑1 (MQ=255)
                                                       gTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCAAGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCAc                                                                                                                                                                                 <  1:258022/251‑1 (MQ=255)
                                                         ccccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTGACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCAGCAATCCGGTTCCGACCGCGACTTACAGTGAAGCGGAAGATCTTAACCCGGCCCCCT‑GCGGCACATTGCCGCTGATCACCATGTCCAcg                                                                                                                                                                                <  1:313225/250‑1 (MQ=255)
                                                             aGTGAGCGGCAATCGAGCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCTCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACgccgcc                                                                                                                                                                           <  1:264674/251‑1 (MQ=255)
                                                               tGAGCGGCAATCGCTCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACACT‑GCGGCACATTGCCGCTGATCACCATGTCCACGccgccg                                                                                                                                                                          >  1:83766/1‑250 (MQ=255)
                                                               tGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGccgccg                                                                                                                                                                          >  1:118979/1‑250 (MQ=255)
                                                                gcgCGGCAATCTATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTATGCGCCGTCCAGCGGCAGATTATAAAGCCGCTGCAATACCGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTATCCCGTCACCCT‑GCGGCACATTGCCGCTGATCCCCATGTCCACGCCGCCGaa                                                                                                                                                                        <  1:547740/249‑1 (MQ=255)
                                                                       aaTCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtgttg                                                                                                                                                                 <  2:129565/251‑1 (MQ=255)
                                                                                        tGATCTTTCTTGCCAAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATg                                                                                                                                                >  1:124405/1‑251 (MQ=255)
                                                                                          aTCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTACGCGCGATTTGTGCGCCGTCAAGCGGCAGATGATAAAGATGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGAGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATgct                                                                                                                                              >  1:357106/1‑249 (MQ=255)
                                                                                              ttcttGCCGAGCGCGGAAATTAGCTTATCCATGATTCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAGTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑TCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGGTGTTACGCAGTTGCAGAGGTTTca                                                                                                                                           <  2:252570/250‑1 (MQ=255)
                                                                                                  tGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGCTGTTTCAccacg                                                                                                                                      >  2:169114/1‑251 (MQ=255)
                                                                                                   gCCGAGCGCGGAAATTAGCTGAGCCATGTACCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGAACGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCGGCTGTAATACGGTTCCGACCGCGACTTACAGGGAAGCGGAAGAACTTTACCTGGC‑CCCTCGCGTCACATTGCCGCTGATCACCATGTCCACGCCGTCGAAGCTGATGTTACTCAGTTGCAgtcgttgacccacg                                                                                                                                      <  2:490407/250‑14 (MQ=255)
                                                                                                    ccGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcc                                                                                                                                    >  1:280898/1‑251 (MQ=255)
                                                                                                            ggAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCg                                                                                                                                                                                                  <  2:18552/181‑1 (MQ=255)
                                                                                                            ggAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCg                                                                                                                                                                                                  >  1:18552/1‑181 (MQ=255)
                                                                                                                   aGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGc                                                                                                                     <  1:106324/251‑1 (MQ=255)
                                                                                                                           cATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTACTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑CCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGAt                                                                                                             >  1:328188/1‑251 (MQ=255)
                                                                                                                              gATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGt                                                                                                          <  1:333651/251‑1 (MQ=255)
                                                                                                                                 cccGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTTATAGtttt                                                                                                       <  2:171526/251‑1 (MQ=255)
                                                                                                                                 cccGCAGTTACAGCCTACAAACTGGTTGTCTGCTTCCTGACCGTTAAGCGCGATTTGGGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATCCGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTTATAGtttt                                                                                                       <  2:131355/251‑1 (MQ=255)
                                                                                                                                  ccGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTc                                                                                                      <  1:116728/251‑1 (MQ=255)
                                                                                                                                    gCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGACTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACg                                                                                                                                                            <  1:537565/195‑1 (MQ=255)
                                                                                                                                    gCAGTAACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGCTTTGTGCGCCGTCCAGCGGCAGATGATAACGCTGATGCAATACGGTTCCGAACGCGACTTACAGTGGCGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGCAGCTGTTGTTACg                                                                                                                                                            >  2:537565/1‑195 (MQ=255)
                                                                                                                                     cAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtg                                                                                                   <  2:17809/251‑1 (MQ=255)
                                                                                                                                      aGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGCGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCCCATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTCCGCAGTTGCAGATGTTTCACACCGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtgg                                                                                                  >  1:340873/1‑250 (MQ=255)
                                                                                                                                         taccgCCTACAAACTGGTTTTCTGCTTACTGACCGTTAAGCGCGATTTGTGCGACGTCCAGCGGCAGATGATAAAGCTGCTGCCATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtg                                                                                                                                                                    >  2:370403/5‑182 (MQ=255)
                                                                                                                                         tACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTGGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGGTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtg                                                                                                                                                                    <  1:370403/182‑1 (MQ=255)
                                                                                                                                          aCAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCg                                                                                               >  2:182073/1‑250 (MQ=255)
                                                                                                                                            aGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCCCCACGCCACGAACGAAGTTAGCCCATTGATAGTTTTCCTGTGCGAc                                                                                             >  1:401354/1‑250 (MQ=255)
                                                                                                                                                 aCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATgg                                                                                        <  2:173873/250‑1 (MQ=255)
                                                                                                                                                  cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTTCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCATTACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCCGCACCCT‑GCGGCACATTGCCGCTGATCGCCATGTCCACGCCGCCGAAGCGGTTGTTACGCAGTTGCATATGTTTCACCACGCCACGAACGTAGTTAGCCCATTTATAGGTTTCATGTGCGACAATGGgc                                                                                      <  2:389492/251‑1 (MQ=255)
                                                                                                                                                        ggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTACCCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCAt                                                                                  <  1:380261/249‑1 (MQ=255)
                                                                                                                                                          ttttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCg                                                                                >  2:455924/1‑249 (MQ=255)
                                                                                                                                                           tttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAgg                                                                             >  1:41701/1‑251 (MQ=255)
                                                                                                                                                           tttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAg                                                                              <  2:118979/250‑1 (MQ=255)
                                                                                                                                                           tttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGGCCAGCGGCAGATGAGAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATCGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAg                                                                              <  2:485408/250‑1 (MQ=255)
                                                                                                                                                                    ccTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGCAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcg                                                                    >  2:551860/1‑251 (MQ=255)
                                                                                                                                                                    ccTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcg                                                                    <  2:153627/251‑1 (MQ=255)
                                                                                                                                                                       gACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtcg                                                                 >  1:340765/1‑251 (MQ=255)
                                                                                                                                                                        aCCGTTAAGCGCGATTTGGGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGCCTTACAGTGAAGCGGAAGCACTTAACCCTGCACCCC‑GCTGCACATTCCCGCTGATCCCCATGCCCACGCCGACGAAGCTGTTGTTACTCAGATGTAAATGTTTCACCACGCTACGAACGTAGTTAGCCCAGTGATAGTGTTCATGTGAGACAATAGGCGCATGTAGCGAAAACGCGTCGa                                                                <  2:545158/251‑1 (MQ=255)
                                                                                                                                                                           gTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGACGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCt                                                             >  2:64141/1‑251 (MQ=255)
                                                                                                                                                                                         gtgCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCg                                               >  1:69057/1‑251 (MQ=255)
                                                                                                                                                                                             gcCGTGCAGGGGCAGAAGATAAAGCGGCTGCAAGACGGATCCGCCCGCGACTTACAGAGAAGCGGAAGAACTTAACCCGGCTCCCT‑GCGGCAAATTGCCGCTGATCACCGTGTCCACGCCGCCGAAGCTGTTGTTACGATGCTGCAGATGTTTCACCACGCCACGAACGTAGTTGGTACATTGATAGTTTTCATGTGCGACCACGGGCGCATAGAGGGAAAACTGGTCAAGCTGATTTTCATAATCgggtg                                           >  2:347140/1‑248 (MQ=255)
                                                                                                                                                                                                gTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGcc                                         >  2:119668/1‑250 (MQ=255)
                                                                                                                                                                                                gTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGcc                                         >  2:107097/1‑250 (MQ=255)
                                                                                                                                                                                                  ccAGCGGCAGATGATAAAGCTGCTGCAATACGGTCCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATc                                      <  1:512304/251‑1 (MQ=255)
                                                                                                                                                                                                   cAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAAc                                                                                                                              <  1:378198/162‑1 (MQ=255)
                                                                                                                                                                                                   cAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAAc                                                                                                                              >  2:378198/1‑162 (MQ=255)
                                                                                                                                                                                                    aGCGGCAGATGATAAAGCTGCTTCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCAGAATCGGCTGCCATCAc                                    <  2:519866/251‑1 (MQ=255)
                                                                                                                                                                                                      cGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcgc                                  >  1:204428/1‑251 (MQ=255)
                                                                                                                                                                                                         cAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGaa                               <  1:513763/251‑1 (MQ=255)
                                                                                                                                                                                                               atAAAGCTGCTGCAATCCGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTAc                         <  1:491284/251‑1 (MQ=255)
                                                                                                                                                                                                                tAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACtttccg                        >  2:151327/1‑248 (MQ=255)
                                                                                                                                                                                                                 aaaGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACg                        <  2:492687/250‑1 (MQ=255)
                                                                                                                                                                                                                          gcAATACGGTTCCGCCCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATcgcg               <  2:437330/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTTCCGACCGCGCCTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGcac       <  2:3825/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCG‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTGTACGGTCAGCGCGTGGTGCACAACTg  <  2:358364/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
GAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCT‑GCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTG  >  NC_000913/789338‑789818

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