Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,026,350 A→G 8.7% E627G (GAA→GGA)  helD → DNA helicase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,026,3500AG 8.7% 130.2 / 2.2 57E627G (GAA→GGA) helDDNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (31/18);  new base G (3/2);  total (35/22)
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.82e-01

GCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAA  >  NC_000913/1026211‑1026492
                                                                                                                                           |                                                                                                                                              
gCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCCGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGCCGAGTATGGAaga                                                                                                                                             >  1:177535/1‑143 (MQ=255)
   acacGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCTGCTGCGGCGCGGGAGTCTATTATGGAAGATgc                                                                                                                                          >  1:297361/1‑143 (MQ=255)
    cacGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAgcggcg                                                                                                                                         >  1:97391/1‑143 (MQ=255)
     acGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGgcgc                                                                                                                                        <  1:363869/143‑1 (MQ=255)
     acGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGgcgc                                                                                                                                        >  1:481666/1‑143 (MQ=255)
         gcgCCGAGGTTGCTAAGCGCCTTTATTACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATTGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACt                                                                                                                                    <  2:69626/141‑1 (MQ=255)
           gCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGc                                                                                                                                  >  1:55219/1‑143 (MQ=255)
                 ggtGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGGTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCgg                                                                                                                            <  2:313822/141‑1 (MQ=255)
                 ggtGCAAATCGACGTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCAGCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCTGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCgg                                                                                                                            <  2:31822/141‑1 (MQ=255)
                      aaaTCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGCAGAGATGGTTTTCCGGCTGCGGCGCGGGAGACGATGATGGCAGAGGCGCTACTGCCACCGGTGGAg                                                                                                                       >  2:97429/1‑143 (MQ=255)
                      aaaTCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGGGTCGATGATGGGAGAGGCGCTACTGCCACCGGGTGAg                                                                                                                       >  2:583881/1‑143 (MQ=255)
                      aaaTCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGCAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGGAGAGGCGCTACTGCCACCGGTTGAg                                                                                                                       >  2:587071/1‑143 (MQ=255)
                               ttATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCAGCTGCCACCGGTTGAGGATTTCCCg                                                                                                              >  2:432418/1‑143 (MQ=255)
                                   gACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGAGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACg                                                                                                          >  2:829276/1‑143 (MQ=255)
                                           aTGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGAGTTCCCGGACGCTgaagaa                                                                                                  >  2:876318/1‑143 (MQ=255)
                                              ccAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGCAGTGCTGGTTTTCCGGCTGCGGCGCGGGGGTCGATTATGGGAGAGGCGACACTGCGCGCGGTTGAGGCTTTCTCGGACGCTgcagcgccg                                                                                               >  2:298120/1‑138 (MQ=255)
                                                          aaCAGGCGGATTTCGTCCTCAGCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGCTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTAc                                                                                   <  2:657778/143‑1 (MQ=255)
                                                             aGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGCCGCTGAAGCAAGGCGGTTAATCTACGTg                                                                                >  1:546881/1‑143 (MQ=255)
                                                             aGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTg                                                                                >  1:453437/1‑143 (MQ=255)
                                                             aGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTg                                                                                >  1:410495/1‑143 (MQ=255)
                                                                cGGATTGCGTTAGCAGCGTGGGCTTGCAGGCGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGgcg                                                                             <  2:425239/143‑1 (MQ=255)
                                                                 ggATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGgcgc                                                                            <  2:55219/143‑1 (MQ=255)
                                                                    ttACGTCATGATCGTTGGCTGGCAGGAGGGAAGTGATGGTGTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGa                                                                         <  2:477759/143‑1 (MQ=255)
                                                                      aCGTGATCATCGTGGGCTTGCAGGAGGGAAGGGGTGGTGTGCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGAcc                                                                       <  2:279831/143‑1 (MQ=255)
                                                                       cttcatcaTCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGAGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATCTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGAccc                                                                      <  2:315362/141‑1 (MQ=255)
                                                                           atcatcGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGc                                                                  >  1:487251/1‑143 (MQ=255)
                                                                             catcGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGTCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCAc                                                                >  2:227327/1‑143 (MQ=255)
                                                                               tcGTTGGCTTGCAGGAGGGAAGGGCTGGTTTCCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGc                                                              <  2:481666/143‑1 (MQ=255)
                                                                                   tGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGAGTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGAGCACGCCATc                                                          >  2:703437/1‑143 (MQ=255)
                                                                                    ggCTTGCGGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCg                                                         <  1:854739/143‑1 (MQ=255)
                                                                                    ggCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATGATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCg                                                         <  2:519807/143‑1 (MQ=255)
                                                                                     gCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGGCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGCAGAACAGCCGCTAATGTTCTTGGCGCTGACCCGGGCACGCCATCgg                                                        <  1:870235/143‑1 (MQ=255)
                                                                                     gCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGGCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCCCTGACCCGGGCACGCCATCgg                                                        >  1:69785/1‑143 (MQ=255)
                                                                                       ttGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGCGTCGAGTATGGAAGAGGCGCTACTGCCCCCGGCTGAGGATTTCCCGGACGCTGAAGACCTGCGGTTAATGTACGAGGCGCTGCCCCGGCCCCGCCACCCGGt                                                      >  2:365558/1‑143 (MQ=255)
                                                                                       ttGCAGGAGGGAAGTGATGAATTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACAGCCACCGGTCGAGGATTGCCCGGACGCTGAAGAACGGAGGTAAACGTACGTGGGGCTGACCCGGGGACGCCATCGGGt                                                      >  2:100133/1‑143 (MQ=255)
                                                                                             gaggGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATGATGGAAGAGGCGCTACTGCCACCGGTGGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTGCGTGGCGCTGGCCCGGGCACGCCATCGGGGATggtc                                                >  2:71045/1‑141 (MQ=255)
                                                                                                ggAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAAAGAACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGcgca                                             >  2:293892/1‑140 (MQ=255)
                                                                                                   aGTGATGGTTTTGCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGtt                                          <  2:33776/143‑1 (MQ=255)
                                                                                                    gTGATGGTTTTCCGGCCGCGGGGCGGCCGTCGACTACGGCAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACCCTGCAGAAAGGCCGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGttt                                         <  1:298120/143‑1 (MQ=255)
                                                                                                      gATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGATGATGTTTTTCTGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTaa                                       >  1:419100/1‑143 (MQ=255)
                                                                                                       aTGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGGTGAGGATTTCCCGGACGGTGAAGAACGGCGGCTAATGTACGTGGCGCCGACCCGGGCACGCCGTAGGGGATGGGCAgtgtgtcac                                      >  2:254043/1‑138 (MQ=255)
                                                                                                           ttttCCGGGTGCGGCGCGGGGGTCGATTATGGGAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAag                                  <  1:583881/143‑1 (MQ=255)
                                                                                                           ttttCCGGCTGCGGCGCGGGAGTCGATTATGGGAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAag                                  <  1:587071/143‑1 (MQ=255)
                                                                                                             ttCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTGGTGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAagag                                >  1:464457/1‑143 (MQ=255)
                                                                                                                        gcgcGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGCTAATGTACCTGGCGCTGAACCGGGCACGCCATCGGGGATGGGCACTGGTTAACACAGAGAATcccgcccc                     >  2:617525/1‑138 (MQ=255)
                                                                                                                        gcgcGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGATACTGTGCGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTGACAAAGAGAATCCctctgc                     >  2:547452/1‑141 (MQ=255)
                                                                                                                            gggAGTCGATTATGGCAGAGGCGCGACTGCCGCCGGTTGGGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCttt                 <  2:699327/143‑1 (MQ=255)
                                                                                                                             ggAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGACTCCCTCTCCCTCtg                >  1:101303/1‑143 (MQ=255)
                                                                                                                               aGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTtgtg              <  1:829276/143‑1 (MQ=255)
                                                                                                                                   gATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTATGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGACTCCCTCTCCCTTTGTGGaaa          >  2:623620/1‑143 (MQ=255)
                                                                                                                                   gATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAaga                                 >  2:294556/1‑120 (MQ=255)
                                                                                                                                   gATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAaga                                 <  1:294556/120‑1 (MQ=255)
                                                                                                                                    aTTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGacgt         >  1:508492/1‑140 (MQ=255)
                                                                                                                                    aTTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAAt         <  1:876318/143‑1 (MQ=255)
                                                                                                                                      tATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATAc       >  2:168158/1‑143 (MQ=255)
                                                                                                                                       aTGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACt      >  1:518170/1‑143 (MQ=255)
                                                                                                                                        tGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTg     >  2:184511/1‑143 (MQ=255)
                                                                                                                                           aaGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAGCTACTGaag  >  1:277057/1‑142 (MQ=255)
                                                                                                                                           aaGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATATTGaaa  <  2:410495/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATCATCGTTGGCTTGCAGGAGGGAAGTGATGGTTTTCCGGCTGCGGCGCGGGAGTCGATTATGGAAGAGGCGCTACTGCCACCGGTTGAGGATTTCCCGGACGCTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAA  >  NC_000913/1026211‑1026492

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: