Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,583,050 G→A Q446* (CAA→TAA)  hsdR ← endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,583,0500GA98.0% 424.9 / NA 148Q446* (CAA→TAA) hsdRendonuclease R
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (84/60);  minor base C (1/2);  total (85/63)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 150 total aligned reads displayed.

CATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGT  >  NC_000913/4582978‑4583125
                                                                        |                                                                           
cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACgg                                                                          <  2:1740251/76‑1 (MQ=255)
 aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                         <  1:2158454/76‑1 (MQ=255)
 aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                         >  2:5660576/1‑76 (MQ=255)
  tCTTCCTGATAATAACGCAGGCCCAGCGCGCTCATGCCAGGGGTATCGGCAAACCACTGATTCTGTTTTTACGGtt                                                                        >  1:1208702/1‑76 (MQ=255)
   cTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTc                                                                       <  1:2689156/76‑1 (MQ=255)
    ttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCg                                                                      <  1:1652216/76‑1 (MQ=255)
     tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGc                                                                     >  1:4617628/1‑76 (MQ=255)
      ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCt                                                                    <  2:3221007/76‑1 (MQ=255)
        tGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                  >  2:2155321/1‑76 (MQ=255)
         gATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcc                                                                 >  2:3674284/1‑76 (MQ=255)
           taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGa                                                               <  2:653624/76‑1 (MQ=255)
           taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCCGCCGa                                                               >  2:714604/1‑76 (MQ=255)
             ataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGc                                                             <  1:2753176/76‑1 (MQ=255)
              taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCa                                                            <  1:1759614/76‑1 (MQ=255)
              taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCa                                                            >  2:2371302/1‑76 (MQ=255)
              taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCa                                                            >  1:3637527/1‑76 (MQ=255)
              taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCa                                                            <  2:3958425/76‑1 (MQ=255)
               aaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAt                                                           >  2:2968554/1‑76 (MQ=255)
               aaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAt                                                           <  1:442724/76‑1 (MQ=255)
               aaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCCGCCGGGCAt                                                           >  2:53576/1‑76 (MQ=255)
               aaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGGTCGCTGCCGAGCAt                                                           >  1:90038/1‑76 (MQ=255)
                 cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                         >  1:5562933/1‑76 (MQ=255)
                 cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                         >  1:4225636/1‑76 (MQ=255)
                  gCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTc                                                        >  2:4139437/1‑76 (MQ=255)
                  gCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTc                                                        >  2:4717932/1‑76 (MQ=255)
                  gCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTc                                                        >  2:5223530/1‑76 (MQ=255)
                   cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTcc                                                       >  2:1330104/1‑76 (MQ=255)
                    aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCca                                                      >  2:732872/1‑76 (MQ=255)
                    aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCca                                                      >  1:4615850/1‑76 (MQ=255)
                     ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcag                                                     >  1:3788433/1‑76 (MQ=255)
                     ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcag                                                     <  2:1735393/76‑1 (MQ=255)
                       cccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagca                                                   <  1:2708050/76‑1 (MQ=255)
                        ccAGCTCTCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                  <  1:1239708/76‑1 (MQ=255)
                        ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCGCTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                  <  2:4113622/76‑1 (MQ=255)
                        ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                  <  1:5269823/76‑1 (MQ=255)
                        ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                  <  1:5092949/76‑1 (MQ=255)
                         cAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                 >  2:3950883/1‑76 (MQ=255)
                          aGCTCGCTCATGCAAGGGTTATCGGCAAACCACTGATTCTGTTTTTCAGGTGCGCTGCCGAGCATTTCCagcagcc                                                >  2:1924124/1‑75 (MQ=255)
                           gctggctcATGCCAGGGTGATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                               <  2:786467/76‑1 (MQ=255)
                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCCGctc                                               >  2:1411966/1‑76 (MQ=255)
                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                               >  1:4997956/1‑76 (MQ=255)
                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                               <  2:1599124/76‑1 (MQ=255)
                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                               <  2:5327182/76‑1 (MQ=255)
                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                               >  2:1273064/1‑76 (MQ=255)
                            ctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctct                                              <  2:2131527/76‑1 (MQ=255)
                              cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                            <  2:1684526/76‑1 (MQ=255)
                               gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCCGCCGAGCATTTCCAGCAGCTCTTcc                                           >  2:3012603/1‑76 (MQ=255)
                                ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCg                                          >  1:2185049/1‑76 (MQ=255)
                                 tcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCgg                                         >  1:4788158/1‑76 (MQ=255)
                                  cATGCCAGGGTTATCGGCAAACCACTGATTCTGCTTTTACGGCTCGCTGCCGAGCATTTTCAGCAGCTCTTCCggg                                        >  2:3969614/1‑76 (MQ=255)
                                   aTGCCAGGGTTATCGGCAAANNACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCATCTCTTCCGGGc                                       >  1:5771138/1‑76 (MQ=255)
                                   aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                       <  1:2968554/76‑1 (MQ=255)
                                   aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                       >  2:3469815/1‑76 (MQ=255)
                                   aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                       <  1:3970511/76‑1 (MQ=255)
                                   aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                       >  1:2519116/1‑76 (MQ=255)
                                    tGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCg                                      >  2:4637155/1‑76 (MQ=255)
                                      gcAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                    <  2:5182941/75‑1 (MQ=255)
                                      ccggggTTACCGACAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                    <  1:787285/73‑1 (MQ=255)
                                      ccAGGGTTATCGGCAAACCGATGATTCTGCTTTTACGGACAGCTGACGCGCACTTCCAGCAGCTCCTCCGTGCggg                                    >  1:3079482/1‑75 (MQ=255)
                                      ccAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                    >  2:3036038/1‑76 (MQ=255)
                                      ccAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                    >  1:3076853/1‑76 (MQ=255)
                                       cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGtc                                   >  2:4782171/1‑75 (MQ=255)
                                       cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGCTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTg                                   >  2:416403/1‑76 (MQ=255)
                                         gggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGcc                                 >  2:3998288/1‑76 (MQ=255)
                                         gggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGAGGTGcc                                 <  1:4890222/76‑1 (MQ=255)
                                            ttATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACt                              <  2:4241781/76‑1 (MQ=255)
                                             tATCGGCAAACCGCTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTc                             <  2:4123401/76‑1 (MQ=255)
                                             tATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTc                             <  2:1175580/76‑1 (MQ=255)
                                             tATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTc                             >  1:3630023/1‑76 (MQ=255)
                                             tATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTc                             <  1:188365/76‑1 (MQ=255)
                                               tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCCGCTCTTACTGGCGTGGGCACTctg                           >  1:3194396/1‑74 (MQ=255)
                                               tCGGCAAACCACTGATTCTGGTGTTACGGTTCGCTGCCGAGCATTGCCAGCAGCTCTTCCGCGCGCTGCCCCTcag                           >  2:3492343/1‑74 (MQ=255)
                                                cGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCGGCTCTTCCGGGCGGTGCCACTCggg                          <  1:1411966/76‑1 (MQ=255)
                                                cGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                          >  2:2675755/1‑76 (MQ=255)
                                                 ggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGt                         <  1:2396494/76‑1 (MQ=255)
                                                  gCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTa                        >  1:5049907/1‑76 (MQ=255)
                                                  gCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTa                        >  1:2380578/1‑76 (MQ=255)
                                                   cAACCCACTGATTCTGTTTTGACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                       <  1:3751782/76‑1 (MQ=255)
                                                   cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                       <  1:3409510/76‑1 (MQ=255)
                                                   cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                       <  1:4219523/76‑1 (MQ=255)
                                                   cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                       <  1:4093278/76‑1 (MQ=255)
                                                    aaaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAg                      >  2:4530981/1‑76 (MQ=255)
                                                     aaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAgg                     <  2:3158552/76‑1 (MQ=255)
                                                      aCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                    >  1:4882011/1‑76 (MQ=255)
                                                      aCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCGCTTCCGGGCGGTGCAACGCGGGTAAGGc                    >  1:4422600/1‑76 (MQ=255)
                                                        cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                  >  1:5786756/1‑76 (MQ=255)
                                                        cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                  <  2:1390624/76‑1 (MQ=255)
                                                          cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                >  1:2842807/1‑76 (MQ=255)
                                                          cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                >  2:5629732/1‑76 (MQ=255)
                                                          cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                >  2:4913982/1‑76 (MQ=255)
                                                            ggttCTGTTTTTCCGGTTCGCTGCCGAGCAGTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGa              <  2:4707945/74‑1 (MQ=255)
                                                             aTTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGAc             >  2:1864341/1‑76 (MQ=255)
                                                             aTTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGAc             <  1:4701015/76‑1 (MQ=255)
                                                             aTTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGAc             <  2:4615849/76‑1 (MQ=255)
                                                               tCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACat           >  1:4969822/1‑76 (MQ=255)
                                                                  gTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGTCACTCGGGTAAGGCTTTCGACATAtt        <  1:1460364/76‑1 (MQ=255)
                                                                  gTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAtt        >  2:114817/1‑76 (MQ=255)
                                                                  gTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAtt        <  1:1865084/76‑1 (MQ=255)
                                                                   tttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTg       <  2:2888483/76‑1 (MQ=255)
                                                                        gCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGt  <  2:372229/76‑1 (MQ=255)
                                                                        |                                                                           
CATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGT  >  NC_000913/4582978‑4583125

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: