Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,002,938 (A)8→7 intergenic (‑137/+48) aspC ← / ← ompF aspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,002,9310A.97.2% 833.9 / 16.5 177intergenic (‑130/+55)aspC/ompFaspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base . (81/91);  total (85/92)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 209 total aligned reads displayed.

ATTTACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGA  >  REL606/1002856‑1003002
                                                                           |                                                                       
aTTTACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGtata                                                                         <  2:3904204/76‑1 (MQ=255)
aTTTACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGtata                                                                         >  2:3172738/1‑76 (MQ=255)
 tttACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaa                                                                        <  1:2303672/76‑1 (MQ=255)
  ttACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaa                                                                       <  2:2899904/76‑1 (MQ=255)
  ttACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaa                                                                       <  2:2628428/76‑1 (MQ=255)
   tACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGATGGTATaaaa                                                                      >  1:2352847/1‑76 (MQ=255)
   tACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaaa                                                                      >  2:4899154/1‑76 (MQ=255)
   tACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaaa                                                                      >  2:516460/1‑76 (MQ=255)
     cAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaaaaa                                                                    <  2:5261885/76‑1 (MQ=255)
       aaaTTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaaaaaac                                                                  >  1:6046746/1‑75 (MQ=255)
       aaaTTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaaaaaac                                                                  <  2:1368031/76‑2 (MQ=255)
        aaTTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATaaaaaaaca                                                                 <  2:3422844/76‑3 (MQ=255)
             cAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGAcc                                                           <  1:4645098/76‑1 (MQ=255)
             cAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGAcc                                                           >  2:2405710/1‑76 (MQ=255)
              aaCTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCa                                                          >  1:1094201/1‑76 (MQ=255)
               aCTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCa                                                          >  2:158950/1‑76 (MQ=255)
               aCTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCaa                                                         <  1:1889447/76‑1 (MQ=255)
                cTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCaaa                                                        <  1:566605/76‑1 (MQ=255)
                cTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCaaa                                                        >  1:929982/1‑76 (MQ=255)
                 tATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAg                                                       >  1:3885137/1‑76 (MQ=255)
                 tATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAg                                                       >  1:3370147/1‑76 (MQ=255)
                  aTTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGt                                                      >  1:6227477/1‑76 (MQ=255)
                   ttACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAACACAGGACCAAAGTc                                                     >  2:4448028/1‑76 (MQ=255)
                    tACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTcc                                                    >  2:856504/1‑76 (MQ=255)
                    tACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTcc                                                    >  2:3518179/1‑76 (MQ=255)
                     aCTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCt                                                   <  1:6275126/76‑1 (MQ=255)
                     aCTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCt                                                   >  1:4180124/1‑76 (MQ=255)
                     aCTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCt                                                   <  2:4057953/76‑1 (MQ=255)
                      cTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTg                                                  <  2:589405/76‑1 (MQ=255)
                      cTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTg                                                  >  2:2743363/1‑76 (MQ=255)
                       tgatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGt                                                 >  1:5026765/1‑76 (MQ=255)
                       tgatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGt                                                 >  2:3218444/1‑76 (MQ=255)
                       tgatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGt                                                 <  2:2877691/76‑1 (MQ=255)
                        gatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtt                                                <  2:386153/76‑1 (MQ=255)
                        gatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtt                                                >  1:1790143/1‑76 (MQ=255)
                        gatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtt                                                >  2:3352742/1‑76 (MQ=255)
                         atgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGttt                                               <  2:4818038/76‑1 (MQ=255)
                         atgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGttt                                               >  1:1710797/1‑76 (MQ=255)
                          tgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtttt                                              <  2:1976162/76‑1 (MQ=255)
                          tgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtttt                                              <  2:2896565/76‑1 (MQ=255)
                          tgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtttt                                              >  2:5985083/1‑76 (MQ=255)
                           gaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGttttt                                             <  2:4026909/76‑1 (MQ=255)
                           gaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGttttt                                             <  1:3441664/76‑1 (MQ=255)
                           gaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGttttt                                             <  2:6088278/76‑1 (MQ=255)
                           gaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGttttt                                             >  1:5619501/1‑76 (MQ=255)
                            aaaaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtttttt                                            <  2:1787926/76‑1 (MQ=255)
                            aaaaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtttttt                                            <  2:2053935/76‑1 (MQ=255)
                            aaaaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGtttttt                                            <  2:3288696/76‑1 (MQ=255)
                               aCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCgg                                         <  2:2352847/76‑1 (MQ=255)
                               aCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCgg                                         >  1:5822004/1‑76 (MQ=255)
                                 tCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCa                                       >  1:499680/1‑76 (MQ=255)
                                 tCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCa                                       >  2:767397/1‑76 (MQ=255)
                                  cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAt                                      >  2:4675519/1‑76 (MQ=255)
                                  cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAt                                      >  2:3361365/1‑76 (MQ=255)
                                    ggCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAttt                                    <  2:5668627/76‑1 (MQ=255)
                                     ggtGTTTTTTCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTa                                   <  2:4705770/74‑1 (MQ=255)
                                      cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTaa                                  >  2:507075/1‑76 (MQ=255)
                                      cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTaa                                  >  1:1319503/1‑76 (MQ=255)
                                      cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTaa                                  <  2:5234436/76‑1 (MQ=255)
                                       tGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                 <  2:4411066/76‑1 (MQ=255)
                                       tGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                 <  2:5011827/76‑1 (MQ=255)
                                       tGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                 <  1:5380490/76‑1 (MQ=255)
                                       tGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                 >  2:5319951/1‑76 (MQ=255)
                                        gTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                <  1:1492484/76‑1 (MQ=255)
                                          ttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                              <  2:416239/76‑1 (MQ=255)
                                          ttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                              <  2:579507/76‑1 (MQ=255)
                                           tttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAag                             <  2:1260794/76‑1 (MQ=255)
                                           tttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAag                             >  1:6120004/1‑76 (MQ=255)
                                           tttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAag                             >  2:5999982/1‑76 (MQ=255)
                                           tttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAag                             <  1:2434764/76‑1 (MQ=255)
                                            ttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAaga                            <  2:1597352/76‑1 (MQ=255)
                                            ttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAaga                            >  1:609508/1‑76 (MQ=255)
                                               cAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                         >  2:898618/1‑76 (MQ=255)
                                                   aCGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGc                     >  2:635522/1‑76 (MQ=255)
                                                    cGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                    >  1:1376395/1‑76 (MQ=255)
                                                    cGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                    >  2:2227323/1‑76 (MQ=255)
                                                      tGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAt                  <  2:18889/76‑1 (MQ=255)
                                                       gagaTTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAtt                 <  2:4503237/76‑1 (MQ=255)
                                                       gagaTTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAtt                 <  1:1210847/76‑1 (MQ=255)
                                                       gagaTTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAtt                 >  1:5670033/1‑76 (MQ=255)
                                                       gagaTTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAtt                 <  1:5962293/76‑1 (MQ=255)
                                                        agaTTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                <  1:2158029/76‑1 (MQ=255)
                                                         gaTTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg               >  1:5157931/1‑76 (MQ=255)
                                                         gaTTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg               <  1:4412347/76‑1 (MQ=255)
                                                           ttGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGaa             <  1:3750343/76‑1 (MQ=255)
                                                           gtGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGaa             <  2:3715855/75‑1 (MQ=255)
                                                             gCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt           >  2:1165976/1‑76 (MQ=255)
                                                              ctctGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTg          >  1:4659401/1‑76 (MQ=255)
                                                              ctctGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTg          >  2:572525/1‑76 (MQ=255)
                                                                 tGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTa       <  1:3951135/76‑1 (MQ=255)
                                                                  ggAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa      <  2:4277599/76‑1 (MQ=255)
                                                                  ggAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa      >  2:1754285/1‑76 (MQ=255)
                                                                  ggAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa      >  1:2628588/1‑76 (MQ=255)
                                                                  ggAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa      >  1:3070376/1‑76 (MQ=255)
                                                                  ggAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa      >  2:112190/1‑76 (MQ=255)
                                                                   gAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa     <  1:5297424/76‑1 (MQ=255)
                                                                   gAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa     >  2:2959890/1‑76 (MQ=255)
                                                                   gAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa     <  2:1695028/76‑1 (MQ=255)
                                                                     aGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt           <  1:4096166/68‑1 (MQ=255)
                                                                      ggTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGa  >  1:3819623/1‑76 (MQ=255)
                                                                           |                                                                       
ATTTACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGA  >  REL606/1002856‑1003002

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