Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,044,3570AC35.1% 301.1 / 69.1 148V284G (GTG→GGG) rfbDdTDP‑4‑dehydrorhamnose reductase subunit, NAD(P)‑binding, of dTDP‑L‑rhamnose synthase
Reads supporting (aligned to +/- strand):  ref base A (86/10);  new base C (0/52);  total (86/62)
Fisher's exact test for biased strand distribution p-value = 3.36e-30
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 166 total aligned reads displayed.

CATCATGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATGGGGTCGA  >  REL606/2044282‑2044431
                                                                              |                                                                          
catcatGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCAcaccca                                                                            <  2:1641595/76‑1 (MQ=255)
   catGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCAcacccacca                                                                         <  2:724875/76‑2 (MQ=255)
   catGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCt                                                                         <  2:272295/76‑1 (MQ=255)
   catGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCt                                                                         >  1:2480703/1‑76 (MQ=255)
   catGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCt                                                                         >  1:2407766/1‑76 (MQ=255)
    atGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTg                                                                        <  2:1838519/76‑1 (MQ=255)
    atGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTg                                                                        >  1:2148948/1‑76 (MQ=255)
     tGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGc                                                                       >  2:1071083/1‑76 (MQ=255)
      gAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGcc                                                                      >  1:4119901/1‑76 (MQ=255)
      gAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGcc                                                                      <  1:216026/76‑1 (MQ=255)
        gCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAg                                                                    <  2:1370737/76‑1 (MQ=255)
         cGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGt                                                                   >  1:726116/1‑76 (MQ=255)
           agaTGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCa                                                                 >  2:1197935/1‑76 (MQ=255)
              tGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGc                                                              >  2:4249719/1‑76 (MQ=255)
               taaaaaaTTTTTTAAATTGCCCGTAGTCGAAAAAATTTGTTTGGCCTTTCTTTCCCCCCCCCCCTGCCAGTCAGGc                                                              <  2:3175544/74‑1 (MQ=255)
               taaaaaaTTTGTTTAAATTGCCGTAGTCGAAAATAATTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGc                                                              <  2:2410759/74‑1 (MQ=255)
                cAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCaa                                                            <  1:2833756/76‑1 (MQ=255)
                cAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCaa                                                            >  2:4141197/1‑76 (MQ=255)
                cAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCaa                                                            >  1:1507792/1‑76 (MQ=255)
                cAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCaa                                                            >  1:3230895/1‑76 (MQ=255)
                  aaatTTTGTAAAATGGCCGTAGTCGAAAATAATTCGTTGAGCTTTCGTTTCCCCCCCCCCTGCCAGTCAGGCAAGa                                                          <  2:3019017/76‑1 (MQ=255)
                  aaaaTTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGa                                                          >  1:4190642/1‑76 (MQ=255)
                   aaaTTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGAc                                                         >  2:3492967/1‑76 (MQ=255)
                    aaTTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACa                                                        >  1:1481914/1‑76 (MQ=255)
                     aTTTTTTAAATTGCCGTAGTCTAAAATAATTCGTGGAGCATCCGTTTCACACCCCCCTGCCAGTCAGGCAAGACaa                                                       <  2:1085421/76‑1 (MQ=255)
                     aTTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACaa                                                       >  1:3955518/1‑76 (MQ=255)
                      tttGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAg                                                      <  2:2663279/76‑1 (MQ=255)
                      tttGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAg                                                      >  1:2335519/1‑76 (MQ=255)
                      tttGTTAAATGGCCGTAGTCAAAAATAATTCGTTGAGCATTCGTTCCACCCCCCCCTGCCAGTCAGGCAAGACAAg                                                      <  1:1459639/76‑1 (MQ=255)
                      tttGTAAAATTGCCGTAGTGGTAAATATTTCGTTGAGCTTTCGTTCCCCCCCCCCCTGCCAGTCAGGCAAGACAAg                                                      <  1:1389543/76‑1 (MQ=255)
                       ttGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAgc                                                     >  2:2311837/1‑76 (MQ=255)
                        tGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAgcg                                                    >  2:712552/1‑76 (MQ=255)
                         ttttAAGTTGCCGTAGTCGGAAATAATCCTTTGGACATTCGTTGCACCCCCCCCTGCCAGTCAGGCAAGACAAgcg                                                    <  2:449367/74‑1 (MQ=255)
                         gTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAgcgc                                                   >  1:1284105/1‑76 (MQ=255)
                           ttAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCa                                                  <  1:1219909/76‑1 (MQ=255)
                           ttAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCa                                                  >  1:288194/1‑76 (MQ=255)
                             aatttGCCGTAGTGTAAATTATTTCTTGGACCATTCGTTCCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCaaa                                                <  1:4173681/73‑1 (MQ=255)
                             aaaTGGCGGTAGTCGAAAATAATTCGTGGAGCTTTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCaaa                                                <  2:343442/76‑1 (MQ=255)
                              aaTTGCCGTAGTCGAAAATATTTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAg                                               <  2:2589196/76‑1 (MQ=255)
                               ttggCCGTAGTCGTAAATATTTCGTTGAGCTTTCTTTTCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGt                                              <  2:2823954/73‑1 (MQ=255)
                               aTTGCCGTAGTCTAAAATAATTCGTTGAGCATTCGTTCCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGt                                              <  1:1589338/76‑1 (MQ=255)
                               aTTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGt                                              >  1:2857091/1‑76 (MQ=255)
                               aTTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGt                                              >  2:2113243/1‑76 (MQ=255)
                                 tGCTGTATTCGTAAATAATTGTTTGACCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTc                                            <  2:4647925/76‑1 (MQ=255)
                                 tGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTc                                            <  2:423269/76‑1 (MQ=255)
                                  gCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTct                                           >  2:2063696/1‑76 (MQ=255)
                                  gCCGTAGTCGAAAATAATTGTTGGAGCATTCTTTCCACACCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTct                                           <  2:3639078/76‑1 (MQ=255)
                                   cccGTAGTCGAAAATAATTCTTTGAGCTTTCTTTCCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTct                                           <  2:4023764/75‑1 (MQ=255)
                                       tAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTctgctg                                       >  2:1403290/1‑76 (MQ=255)
                                        aGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGa                                      >  2:1973466/1‑76 (MQ=255)
                                        aGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGa                                      >  1:3841988/1‑76 (MQ=255)
                                        aGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGa                                      >  2:291480/1‑76 (MQ=255)
                                         gTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGaa                                     >  2:2255594/1‑76 (MQ=255)
                                          tCGAAAATATTTCGTGGAGCATTCTTTTCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGaaa                                    <  2:531565/76‑1 (MQ=255)
                                            gTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACt                                  >  1:3386730/1‑76 (MQ=255)
                                            gTAAATAATTCGTGGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACt                                  <  2:1192195/76‑1 (MQ=255)
                                            atAAATAATTCGTTGAGCTTTCTTTCCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACt                                  <  1:298400/75‑1 (MQ=255)
                                             tAAATATTTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACtt                                 <  2:650679/76‑1 (MQ=255)
                                             tAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACtt                                 >  2:3199481/1‑76 (MQ=255)
                                             tAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACtt                                 >  2:3779397/1‑76 (MQ=255)
                                              aaaTAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACttt                                >  2:2498/1‑76 (MQ=255)
                                               aataatTCGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACtttt                               <  2:1280005/76‑1 (MQ=255)
                                               aataatTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACtttt                               >  2:3751360/1‑76 (MQ=255)
                                               aataatTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACtttt                               >  1:832013/1‑76 (MQ=255)
                                               aataatTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACtttt                               >  1:4792828/1‑76 (MQ=255)
                                                atattTCGTGGAGCATCCGTTCCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTc                              <  1:4608004/76‑1 (MQ=255)
                                                ataatTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTc                              >  2:1045191/1‑76 (MQ=255)
                                                    atTCGTTGAGCTTTCTTTCCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGt                           <  1:3884050/76‑1 (MQ=255)
                                                    atTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGt                           >  1:12692/1‑76 (MQ=255)
                                                     tgcGTTGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTa                          <  1:2541799/74‑1 (MQ=255)
                                                     tgcGTGGAGCATTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTa                          <  2:4901470/74‑1 (MQ=255)
                                                       gtttGACCTTTCGTTCCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTAtt                        <  2:2611841/74‑1 (MQ=255)
                                                       cGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTAtt                        >  1:4895411/1‑76 (MQ=255)
                                                       cGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTAtt                        >  2:3612339/1‑76 (MQ=255)
                                                       cGTTGACCTTTCGTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTAtt                        <  2:3538992/76‑1 (MQ=255)
                                                        gtcgAGCATTCTTTTCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTg                       <  1:1217110/73‑1 (MQ=255)
                                                          tGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAg                     >  1:2589818/1‑76 (MQ=255)
                                                          tGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAg                     >  2:4280832/1‑76 (MQ=255)
                                                          tGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAg                     >  2:2340364/1‑76 (MQ=255)
                                                           gAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAgg                    >  1:5056402/1‑76 (MQ=255)
                                                              cATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCga                 >  1:3711494/1‑76 (MQ=255)
                                                               aTTCGTTCCACCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCgag                <  1:983174/76‑1 (MQ=255)
                                                                ttCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCgaga               >  2:1507294/1‑76 (MQ=255)
                                                                ttCGTTCCCCCCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCgaga               <  1:1710827/76‑1 (MQ=255)
                                                                 tCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGagag              >  1:3280319/1‑76 (MQ=255)
                                                                    tttCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTa           >  2:4794290/1‑76 (MQ=255)
                                                                    tttCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTa           >  2:2284294/1‑76 (MQ=255)
                                                                     ttCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTAt          >  1:2031180/1‑76 (MQ=255)
                                                                     ttCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTAt          >  1:4428624/1‑76 (MQ=255)
                                                                     tccaccCCCCCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTAt          <  1:4661262/74‑1 (MQ=255)
                                                                      tCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATg         >  2:4725752/1‑76 (MQ=255)
                                                                      tCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATg         >  2:2193655/1‑76 (MQ=255)
                                                                      tCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATg         >  1:4053653/1‑76 (MQ=255)
                                                                      tCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATg         >  2:84448/1‑76 (MQ=255)
                                                                       cacaCCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATgg        >  1:2767164/1‑76 (MQ=255)
                                                                         ccccccccTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATgggg      <  2:493853/74‑1 (MQ=255)
                                                                         ccccccccTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATgggg      <  2:5030358/74‑1 (MQ=255)
                                                                         cacccaccTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATgggg      >  2:1950144/1‑76 (MQ=255)
                                                                          cccccccTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATGGGGt     <  2:1634458/75‑1 (MQ=255)
                                                                             caccTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATGGGGTCGa  >  1:4982790/1‑76 (MQ=255)
                                                                              |                                                                          
CATCATGAGCGAGATGCAAAAATTTGTTAAATTGCCGTAGTCGTAAATAATTCGTTGAGCATTCGTTTCACACCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAACTTTTCTGTATTGAGGCGAGAGTTATGGGGTCGA  >  REL606/2044282‑2044431

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