Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,823,288 C→T G161D (GGC→GAC)  recA ← DNA recombination/repair protein RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,2880CT94.6% 235.0 / 2.9 74G161D (GGC→GAC) recADNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (34/36);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

GCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATC  >  NC_000913/2823047‑2823530
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
gCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGc                                                                                                                                                                                                                                           >  1:359022/1‑251 (MQ=255)
 cTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGcc                                                                                                                                                                                                                                          >  1:310165/1‑251 (MQ=255)
  tACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCt                                                                                                                                                                                                                                         >  2:18758/1‑251 (MQ=255)
         cacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAttt                                                                                                                                                                                                                                  <  1:404688/251‑1 (MQ=255)
            gTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCg                                                                                                                                                                                                                               >  1:307696/1‑251 (MQ=255)
                         ttCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTc                                                                                                                                                                                                                    <  1:356396/249‑1 (MQ=255)
                              cgcgCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgccg                                                                                                                                                                                                             >  2:223009/1‑251 (MQ=255)
                                   cGATACGACGGATGTCGAGACAAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACg                                                                                                                                                                                                        <  1:124012/251‑1 (MQ=255)
                                    gctACGACGGATGTCGAGACGAACCGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCCATTTTCATACGGATCTGGTTGCTGAAGATCAGCAGCGTGTTGGCCTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCACGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACgg                                                                                                                                                                                                       >  2:278815/3‑251 (MQ=255)
                                       acgacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAg                                                                                                                                                                                                     <  1:131115/250‑1 (MQ=255)
                                           cgGATGTCGTGCCGAACAGAGGCGTAGAAGTTCCGCGCGTTCCCACCGGTAGGGGTTTCCGGGTTACCGAACATCACCCCAATTTTCATACGGATCTGGTTGATGAAGATCAGCGGCGTGTTGGCCTGCTGCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAac                                                                                                                                                                                                <  2:251599/251‑1 (MQ=255)
                                            ggATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCTTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAacg                                                                                                                                                                                               <  1:339596/251‑1 (MQ=255)
                                            ggATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAacg                                                                                                                                                                                               >  1:473117/1‑251 (MQ=255)
                                              aTGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCTGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAacgac                                                                                                                                                                                             >  1:153557/1‑251 (MQ=255)
                                                  cGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGAt                                                                                                                                                                                          >  1:522309/1‑250 (MQ=255)
                                                    agaCGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACCCCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAAc                                                                                                                                                                                       <  2:338232/251‑1 (MQ=255)
                                                           cAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAATATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACt                                                                                                                                                                                <  1:347679/251‑1 (MQ=255)
                                                            agagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTg                                                                                                                                                                               >  1:223011/1‑251 (MQ=255)
                                                              agGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTgcg                                                                                                                                                                             <  1:324782/251‑1 (MQ=255)
                                                                      aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAAcc                                                                                                                                                                     >  1:378535/1‑250 (MQ=255)
                                                                         ttCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACgcgc                                                                                                                                                                  >  1:521931/1‑251 (MQ=255)
                                                                               gcgTTCCCCCCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACt                                                                                                                                                                                <  2:262326/231‑1 (MQ=255)
                                                                               gcgTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACt                                                                                                                                                                                >  1:262326/1‑231 (MQ=255)
                                                                               gcgTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGc                                                                                                                                                            >  1:556047/1‑251 (MQ=255)
                                                                                cgTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCg                                                                                                                                                           >  2:393883/1‑251 (MQ=255)
                                                                                      caccGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTcac                                                                                                                                                       >  2:9350/1‑249 (MQ=255)
                                                                                         cGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTCTCATACGGATCTGGTTGATGAAGAGCAGCAGCGTGTTGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAtt                                                                                                                                                  <  1:49348/251‑1 (MQ=255)
                                                                                            tAGTGGTTCCCGGTTTCACGAACTTCACACCAATTTTCATACGGATCTGGTTGATGAAGTTCAGAATCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGGGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTc                                                                                                                                                <  1:76526/250‑1 (MQ=255)
                                                                                                ggTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGTTCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTg                                                                                                                                           <  2:522309/251‑1 (MQ=255)
                                                                                                ggTTTCCGGGTGACCAAACTTCCCACCATTTTTCATACGGATCTGGTTGATGAAAACCAGCACCGTGTTGGATGGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTg                                                                                                                                           <  2:83535/251‑1 (MQ=255)
                                                                                                  tttCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAGCGCGCCAGGGCGTCACAGATTTCCAGTGcc                                                                                                                                         >  1:462680/1‑251 (MQ=255)
                                                                                                   ttCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGcc                                                                                                                                         <  1:523806/250‑1 (MQ=255)
                                                                                                    tCCGGGTTCACGACCATCACCCCAATTTTCATACGGATCGGGTTGATGAAGAACAGCATGGTGTTGGACTTCTCCAGGTTACCCGCCAGCTTCCGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCt                                                                                                                                        <  2:45607/250‑1 (MQ=255)
                                                                                                      cggggtaccgaacattgctcataTTTTCACACGGATCTGGGTGATGAATGTCAGCAGCGTGGGGGTCTCCTGCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATGCGTGCCGCAAGTCCCATGTGAGAGTCGTCGATTTCGCCTACGATTTCCGCTTTCGGCGTCATTGCCGCCGCGGAGTCAACGACGATAACGTCTGCTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                     <  2:565193/229‑1 (MQ=255)
                                                                                                      cGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                     >  1:81964/1‑251 (MQ=255)
                                                                                                       gggTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGccc                                                                                                                                    >  2:74733/1‑249 (MQ=255)
                                                                                                       gggTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTc                                                                                                                                    >  1:466541/1‑251 (MQ=255)
                                                                                                                aaCATCACACCAATTTTCATACGGATCGGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTc                                                                                                                           <  1:2672/251‑1 (MQ=255)
                                                                                                                  cATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGCACGCGCCAGGGCGTCACAGTTTTCCAGTGCCTGCTCGCCGGAGTCCg                                                                                                                         >  2:152937/1‑251 (MQ=255)
                                                                                                                        acaaaTTTCCATACGGATTTGGTTGATGAAGATCAGTAGCGGTTTGGCCTGTTTAAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTggg                                                                                                                   <  2:579302/248‑1 (MQ=255)
                                                                                                                         ccAATTTGCATACGGATCTGGTTGAGGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCGTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGa                                                                                                                  <  1:581380/251‑1 (MQ=255)
                                                                                                                           aaTTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGGGTCCGGCTGGGAGc                                                                                                                >  1:275426/1‑251 (MQ=255)
                                                                                                                             ttttCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGGGcac                                                                                                              >  2:418185/1‑251 (MQ=255)
                                                                                                                                  gctcggtgcTGGTTGATGAAGATAGGCGTCGTGTGGGACTGCTTCGGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCAcagcag                                                                                                         <  2:303561/243‑1 (MQ=255)
                                                                                                                                  aTACGGATCTGGTTGATGAAGATCAGCAGAGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCAcagc                                                                                                           >  1:149459/1‑249 (MQ=255)
                                                                                                                                     cGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTCCGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAATCGTCGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCCCAGCAGGtt                                                                                                      >  1:374778/1‑251 (MQ=255)
                                                                                                                                       gtccTGGTTGATGAAGAGCGGCAGCGTGTTGCACTGCGTCGGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGt                                                                                                    <  2:310165/248‑1 (MQ=255)
                                                                                                                                                gatgaAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCCGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCCGGTTGTCGATTTCgc                                                                                           >  2:506539/1‑250 (MQ=255)
                                                                                                                                                 atgaAGTTCAGCAGCGTGTTGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGAc                                                                                          <  2:484488/251‑1 (MQ=255)
                                                                                                                                                    aaGATAAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTAGCCGGTGTCCGGCTGGGGGCACAGCAGGTTTTCGATATCGACGcc                                                                                       >  2:479769/1‑251 (MQ=255)
                                                                                                                                                      gATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCa                                                                                     <  1:240264/251‑1 (MQ=255)
                                                                                                                                                               cGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTg                                                                            >  1:161942/1‑251 (MQ=255)
                                                                                                                                                               cGTGTTGGACTGCTGCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTg                                                                            <  1:230058/251‑1 (MQ=255)
                                                                                                                                                                          gCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTg                                                                   <  1:242937/249‑1 (MQ=255)
                                                                                                                                                                             tCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTcc                                                              >  1:405592/1‑251 (MQ=255)
                                                                                                                                                                             tCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTc                                                               <  1:393883/250‑1 (MQ=255)
                                                                                                                                                                                 gttccccGCCAGCTTTCGCATCGCCTGGCTCATCATACGTTCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAGTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgcg                                                          <  1:64529/247‑1 (MQ=255)
                                                                                                                                                                                 gTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgcg                                                          >  2:484191/1‑251 (MQ=255)
                                                                                                                                                                                     cccGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCgtg                                                      <  1:152937/251‑1 (MQ=255)
                                                                                                                                                                                          cAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTc                                                   >  1:464964/1‑249 (MQ=255)
                                                                                                                                                                                                  gCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCCTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATa                                         <  1:217050/251‑1 (MQ=255)
                                                                                                                                                                                                    aTCGCCTGGGTCATCATACGGGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCg                                            <  2:314570/246‑1 (MQ=255)
                                                                                                                                                                                                    aTCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCg                                            >  1:314570/1‑246 (MQ=255)
                                                                                                                                                                                                       gCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGc                                    >  1:464829/1‑251 (MQ=255)
                                                                                                                                                                                                        ccTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGcc                                   >  2:68463/1‑250 (MQ=255)
                                                                                                                                                                                                           ggCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAg                                <  1:566657/251‑1 (MQ=255)
                                                                                                                                                                                                            gCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAgg                               <  1:114637/251‑1 (MQ=255)
                                                                                                                                                                                                                         gCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcg                  >  1:374042/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTCCCTTCAcgc                 >  2:78222/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctg               <  2:473117/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cccATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACgctgcgctgc         <  1:223009/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGctgcgctgcg        <  1:483213/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gagaGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGAt   <  1:18758/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gagaGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGAGGCGATc  >  2:542197/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
GCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATC  >  NC_000913/2823047‑2823530

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