Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,985,084 A→G V437A (GTA→GCA)  ECB_03714 ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,985,0840AG98.8% 695.1 / ‑3.1 170V437A (GTA→GCA) ECB_03714conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (55/113);  total (55/115)
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.60e-01

Only 100 of 203 total aligned reads displayed.

TTATTCATGTCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTACCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGTCG  >  REL606/3985012‑3985156
                                                                        |                                                                         
ttATTCATGTCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCa                                                                        <  2:3460812/76‑1 (MQ=255)
ttATTCATGTCGATGTGTGAAACTTGGTCCAGCAATTCGATTGCAGCGGTAGTGAACGCGCCATAACCAAGTGCCa                                                                        <  2:5052790/76‑1 (MQ=255)
 tATTCATGTCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAg                                                                       <  2:5941885/76‑1 (MQ=255)
  aTTCATGTCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAgg                                                                      <  2:4102295/76‑1 (MQ=255)
     cATGTCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAAACAAGTGGCAGGTCg                                                                   >  1:4029097/1‑76 (MQ=255)
        gTCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGGAGTGGGAGGGCCCTAAACACGGGCCAGGGCGaaa                                                                >  1:4946696/1‑75 (MQ=255)
         tCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGGGTAGTGGGAGGGCCCCAAACAAGTGCCAGGGCGGATa                                                              >  1:6147279/1‑76 (MQ=255)
         tCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATa                                                              <  1:3032444/76‑1 (MQ=255)
          cGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATaa                                                             <  1:1175673/76‑1 (MQ=255)
          cGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATaa                                                             <  1:265658/76‑1 (MQ=255)
          cGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGGGGAAGGGCCCACAACAAATGCCAGGGCGAATaa                                                             >  2:2283232/1‑76 (MQ=255)
           gATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAac                                                            <  2:328238/76‑1 (MQ=255)
           gATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTACCAGGTCGAATaat                                                            <  2:1444006/76‑2 (MQ=255)
            aTGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAaca                                                           <  1:626917/76‑1 (MQ=255)
            aTGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAaca                                                           <  1:5837311/76‑1 (MQ=255)
             tgtgtgAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAacac                                                          <  2:747379/76‑1 (MQ=255)
             tgtgtgAAAGTTGGTCCAGTAATTCGCTGGCAGGCGGAGTGGGAGGGCCCTAACCAAGTGCCCGGGCGAAAAacac                                                          >  2:4342198/1‑76 (MQ=255)
              gtgtgAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACacg                                                         <  1:5622529/76‑1 (MQ=255)
              gtgtgAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGGGGAAGGGCCATAAACAAGGTCCCGGGCGAACTAacaa                                                          >  2:4580423/1‑75 (MQ=255)
                gtgAAAGTTTGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACacgac                                                       <  1:4990478/76‑1 (MQ=255)
                gtgAAAGTTGGTCCAGTAATTCGCTGGCAGGCGGAGTGGGAGCGCCCTAACCAAATGCCAGGGCGGATAACaagac                                                       >  1:1854283/1‑76 (MQ=255)
                 tgAAAGTTGGTCCAGTAATTCGCTGGCAGGGGTAGTGGGAGCGCCATAAACCAGTGCCCGGGCGAATAACacgccg                                                      >  2:4551672/1‑76 (MQ=255)
                 tgAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACacgacg                                                      <  1:1075677/76‑1 (MQ=255)
                 tgAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACacgacg                                                      <  2:1625960/76‑1 (MQ=255)
                  gAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACgg                                                     <  2:4299778/76‑1 (MQ=255)
                     aGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGGGCCCTAAACAAGGGCCAGGGCGGATAAAAAGACGGCAg                                                  >  1:4091392/1‑76 (MQ=255)
                      gTTGGTCCAGTAATTCGCTGGCAGGGGGAGTGGAAGGGCCATAAACAAGTGCCACGTCGAAAAACACAACGGcggg                                                 >  2:5222280/1‑73 (MQ=255)
                      gTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCCTAAACAAGAGCCAGGGCGAATAACAAGACGGCAgg                                                 >  1:3740339/1‑75 (MQ=255)
                       ttGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAgtg                                                <  2:5759125/76‑1 (MQ=255)
                        tGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAgtgt                                               <  2:134542/76‑1 (MQ=255)
                        tGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAgtgt                                               <  1:1057267/76‑1 (MQ=255)
                           tCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGtgtgtg                                            <  1:3515683/76‑1 (MQ=255)
                            ccAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGa                                           <  2:141044/76‑1 (MQ=255)
                            ccAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGa                                           <  1:4311379/76‑1 (MQ=255)
                             cAGTAATTCGCTGGCAGGCGTAGTGGAAGGGCCATAAACAAAGGCCCGGGCGGATAACACGACGGCCGggtgtggg                                          >  2:4486189/1‑74 (MQ=255)
                             cAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAg                                          <  1:248960/76‑1 (MQ=255)
                               gTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTc                                        <  1:1417151/76‑1 (MQ=255)
                               gTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTc                                        <  1:1496546/76‑1 (MQ=255)
                                  aTTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGtt                                     <  2:4240231/76‑1 (MQ=255)
                                  aTTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGtt                                     <  2:3989784/76‑1 (MQ=255)
                                   ttCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTg                                    <  1:4567026/76‑1 (MQ=255)
                                    tCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGa                                   <  2:1045920/76‑1 (MQ=255)
                                      gCTGGTAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGt                                 <  2:5539736/76‑1 (MQ=255)
                                      gCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGt                                 <  1:3636851/76‑1 (MQ=255)
                                         ggCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCg                              <  2:338494/76‑1 (MQ=255)
                                         ggCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCg                              <  2:2871592/76‑1 (MQ=255)
                                         ggCAGGAGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCg                              <  2:94244/76‑1 (MQ=255)
                                          gCAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGt                             <  2:2991222/76‑1 (MQ=255)
                                           cAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGtt                            <  1:12265/76‑1 (MQ=255)
                                           cAGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGtt                            <  1:6036356/76‑1 (MQ=255)
                                            aGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTg                           <  1:3598235/76‑1 (MQ=255)
                                            aGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTg                           <  1:5810897/76‑1 (MQ=255)
                                            aGGCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTg                           <  1:5555650/76‑1 (MQ=255)
                                              gCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAg                         <  1:907874/76‑1 (MQ=255)
                                              gCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAg                         >  1:5043950/1‑76 (MQ=255)
                                              gCGTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAg                         <  2:3936936/76‑1 (MQ=255)
                                                gTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCa                       <  1:5244110/76‑1 (MQ=255)
                                                gTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCa                       >  1:2460810/1‑76 (MQ=255)
                                                gTAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCa                       >  2:1711099/1‑76 (MQ=255)
                                                 tAGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAg                      <  2:5926701/76‑1 (MQ=255)
                                                  aGTGGAAGCGCCCTAACCAAGTGCCAGGTCGAATAACACGACGGCAGAGGGTGAGTCGTTGAGTTCGTTTAGCAGa                     >  2:2067913/1‑76 (MQ=255)
                                                  aGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGa                     >  1:5519469/1‑76 (MQ=255)
                                                  aGTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGa                     <  1:2599568/76‑1 (MQ=255)
                                                   gTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGAc                    <  2:5657473/76‑1 (MQ=255)
                                                   gTGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGAc                    <  2:1153007/76‑1 (MQ=255)
                                                    tGGAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACg                   >  1:2520962/1‑76 (MQ=255)
                                                     ggAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGt                  >  2:5667738/1‑76 (MQ=255)
                                                     ggAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGt                  >  1:4093779/1‑76 (MQ=255)
                                                      gAAGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGtt                 <  2:6155877/76‑1 (MQ=255)
                                                       aaGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTg                >  1:1136403/1‑76 (MQ=255)
                                                       aaGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTg                <  1:1821196/76‑1 (MQ=255)
                                                       aaGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTg                <  1:4864829/76‑1 (MQ=255)
                                                       aaGCGCCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTg                >  1:2922270/1‑76 (MQ=255)
                                                          cgcCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGAtt             >  2:2557660/1‑76 (MQ=255)
                                                          cgcCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGAtt             <  1:3183257/76‑1 (MQ=255)
                                                          cgcCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGAtt             <  1:1755885/76‑1 (MQ=255)
                                                           gcCATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTg            <  2:2692966/76‑1 (MQ=255)
                                                            ccATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgc           <  1:517151/76‑1 (MQ=255)
                                                            ccATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgc           >  2:2788451/1‑76 (MQ=255)
                                                            ccATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgc           >  2:1934711/1‑76 (MQ=255)
                                                             cATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgcg          >  1:388504/1‑76 (MQ=255)
                                                             cATAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgcg          <  1:2211641/76‑1 (MQ=255)
                                                              aTAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgcgc         <  2:4064436/76‑1 (MQ=255)
                                                              aTAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgcgc         <  1:3976104/76‑1 (MQ=255)
                                                              aTAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgcgc         >  1:2132123/1‑76 (MQ=255)
                                                              aTAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgcgc         >  2:3853835/1‑76 (MQ=255)
                                                              aTAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTgcgc         >  2:5181595/1‑76 (MQ=255)
                                                               tAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCa        >  1:2334170/1‑76 (MQ=255)
                                                               tAACCAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCa        >  1:5700584/1‑76 (MQ=255)
                                                                  ccAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCagag     >  2:5469131/1‑76 (MQ=255)
                                                                  ccAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCagag     <  2:5492025/76‑1 (MQ=255)
                                                                  ccAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCagag     >  2:2445554/1‑76 (MQ=255)
                                                                  ccAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCagag     <  1:5364231/76‑1 (MQ=255)
                                                                   cAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGt    >  2:4853376/1‑76 (MQ=255)
                                                                   cAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGt    >  2:3233076/1‑76 (MQ=255)
                                                                   cAAGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGt    <  2:6282677/76‑1 (MQ=255)
                                                                    aaGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGTc   <  1:3117860/76‑1 (MQ=255)
                                                                    aaGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGTc   >  1:1646511/1‑76 (MQ=255)
                                                                     aGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGTCg  >  1:2087170/1‑76 (MQ=255)
                                                                     aGTGCCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGTCg  >  1:5956289/1‑76 (MQ=255)
                                                                        |                                                                         
TTATTCATGTCGATGTGTGAAAGTTGGTCCAGTAATTCGCTGGCAGGCGTAGTGGAAGCGCCATAACCAAGTACCAGGTCGAATAACACGACGGCAGTGTGTGAGTCGTTGAGTTCGTTGAGCAGACGTTGATTGCGCAGAGTCG  >  REL606/3985012‑3985156

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