Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,854 T→G 6.7% intergenic (‑81/+37) dcuA ← / ← aspA C4‑dicarboxylate transporter DcuA/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,8540TG 6.7% 151.7 / 4.2 60intergenic (‑81/+37)dcuA/aspAC4‑dicarboxylate transporter DcuA/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base G (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-02

CCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACC  >  NC_000913/4366713‑4366992
                                                                                                                                             |                                                                                                                                          
cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTGCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                           <  1:909761/143‑1 (MQ=255)
 ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACACGACAAAAGGCAAGTCATCt                                                                                                                                          >  1:612138/1‑143 (MQ=255)
 ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCt                                                                                                                                          >  2:446532/1‑143 (MQ=255)
  cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTg                                                                                                                                         >  1:144445/1‑143 (MQ=255)
  cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTCGTTTTTTAACATGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAACGGCACGTCATCTg                                                                                                                                         >  1:808194/1‑143 (MQ=255)
   ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGCTCTACAACTAGCCTATTAGCCTTCCTTGTCTTTTGACAAGTTGAGATTAGATTGTTATTTTTAAGTTACTACTAACACTAAAACAGACACGTCATCAGa                                                                                                                                        >  2:13955/1‑143 (MQ=255)
      aTCGGGCGACCAGGAAGATCGCCAGCAAAACTATGATGGATTCTACAACGAGCATATTAGCCTTCCTTGTGTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAAAAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                     <  2:297547/143‑1 (MQ=255)
            cgcCCAAGAAGCTCGCCAGCAAAACTCTGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAACTTGTTAGTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCttt                                                                                                                               <  1:696999/143‑1 (MQ=255)
             gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCtttt                                                                                                                              <  1:594319/143‑1 (MQ=255)
               ccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCCTATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCtttttt                                                                                                                            <  2:554926/143‑1 (MQ=255)
                cAAGGAGATCGCCTGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCttttttc                                                                                                                           <  2:315706/143‑2 (MQ=255)
                cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGGAAAAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                           <  2:672377/143‑1 (MQ=255)
                 aagaagATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                          >  1:66067/1‑143 (MQ=255)
                  agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACACGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTattatt                                                                                                                         >  2:206005/1‑143 (MQ=255)
                  agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                         >  1:618685/1‑143 (MQ=255)
                  agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCTTTTTTCtt                                                                                                                         <  1:666874/143‑1 (MQ=255)
                          gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATCTTAGATTGTTATTTTTAAGTTACTGCTCCCAAGAAAAAAGGCACGTCATGTGACGGGCCGTTTTTTTTTGtactac                                                                                                                 >  2:825324/1‑143 (MQ=255)
                          gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtactac                                                                                                                 >  1:9646/1‑143 (MQ=255)
                                         atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGCAAAAAGGCACGTCATCTGACCTCGCTTTTTTATTGGTAATCCCGTGTACGCTTACTg                                                                                                  >  2:167493/1‑143 (MQ=255)
                                         atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTg                                                                                                  >  1:382272/1‑143 (MQ=255)
                                              ttCTACAACTAGCATATTAGCCTTCCTAGTTTTTTAACAAGATGATATTAGATCGTTATTTTTAAGTTACTCATCACAAGAAAAGAGGCACGACATCTGACGTGCCTCTTTTATTTGTACTACACTGTACGATTACTATTCGc                                                                                             >  2:150360/1‑143 (MQ=255)
                                               tCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                            >  1:1028291/1‑143 (MQ=255)
                                                cTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACACTGTACGATTACTGTTCGCtt                                                                                           >  2:547586/1‑143 (MQ=255)
                                                  aCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtc                                                                                         <  1:927938/143‑1 (MQ=255)
                                                   cAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                        <  2:302253/143‑1 (MQ=255)
                                                    aaCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                       <  1:574032/143‑1 (MQ=255)
                                                    aaCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                       >  1:26279/1‑143 (MQ=255)
                                                    aaCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                       >  1:563479/1‑143 (MQ=255)
                                                      cgaGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                     <  2:704305/141‑1 (MQ=255)
                                                      cgaGCATATTAGCCTGTCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                     <  1:1054583/141‑1 (MQ=255)
                                                      cTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAAACTATACGATTACTGTTCGCTTtcatca                                                                                     <  2:396364/143‑1 (MQ=255)
                                                       tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                    >  1:738334/1‑143 (MQ=255)
                                                          cATATTAGCCTTCCTTGTTGTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtat                                                                                 <  2:630911/143‑1 (MQ=255)
                                                                     tccctGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGCACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCt                                                                      <  1:192579/143‑1 (MQ=255)
                                                                         tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTATACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGt                                                                  >  1:78960/1‑143 (MQ=255)
                                                                                tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTGTAGCGTTTTGCTTTGTAAGCCgg                                                           >  1:305412/1‑143 (MQ=255)
                                                                                    aaGCTGATCTTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTCGACACGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGc                                                       <  1:696700/143‑1 (MQ=255)
                                                                                      gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                     >  1:1027176/1‑143 (MQ=255)
                                                                                        tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCa                                                   <  1:446532/143‑1 (MQ=255)
                                                                                               aGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGGGCATCAGATTCTg                                            >  1:464340/1‑143 (MQ=255)
                                                                                                gATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCTTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                           <  2:239150/143‑1 (MQ=255)
                                                                                                     ttATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGga                                      >  1:46786/1‑143 (MQ=255)
                                                                                                         ttttAAGTTTCTGCTCACAAGAAAAAAGGCACGTCATCTGCCGTGCCTTTTTTATTTGTCCTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACTGAGaaa                                  <  2:547737/143‑1 (MQ=255)
                                                                                                          tttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTAGCTTTCATCAGTATAGCGTTTTGCTTTGTACCCCGGGTGCATGAGATGCTGTACGGAGaaaa                                 >  2:774920/1‑143 (MQ=255)
                                                                                                            tAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGCAAata                               >  2:1050863/1‑143 (MQ=255)
                                                                                                                          caAGAAAAACGGCACGTCATCTGACGTGCCTTTTGTATTTGTACTGCCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCGTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                 <  2:57482/143‑1 (MQ=255)
                                                                                                                           aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATTTACCGCTATACGGAGAAAATATCGTCAGGTTCCGCt                >  2:228728/1‑143 (MQ=255)
                                                                                                                              cgaaaaGGCACCGCCTCTGACGTGCCTTTTTTATTTGTACTAACCTGGACGTGTACTGTTCGCTTTCATCATTATTGCGTTTTTCTCTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa             <  1:206005/141‑1 (MQ=255)
                                                                                                                              aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa             <  1:928593/143‑1 (MQ=255)
                                                                                                                              aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGATCGGTTTCATCAGTATAGCATTTTGCTTTGTCAGCCGGGTGCATCAGATTATGTACGGAGAAAATGTAGTCAAGTACCgagtaa             >  2:635340/1‑138 (MQ=255)
                                                                                                                                 aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTACAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc          <  2:563479/143‑1 (MQ=255)
                                                                                                                                   aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCaa        <  2:923000/143‑1 (MQ=255)
                                                                                                                                   aGGCACGCCATCTGACGTGCCTTTTTTATTTGTGCTACCCGGTGCGATTACCGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGCTTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCaa        <  2:178910/143‑1 (MQ=255)
                                                                                                                                    ggCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc       >  1:1029899/1‑143 (MQ=255)
                                                                                                                                    ggCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGCGAAAATATCGTCAAGTTCCGCTTCAGTCAAc       >  1:648446/1‑143 (MQ=255)
                                                                                                                                     gCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAGCAGATTCTGTACGGAGAAAATATCGTCAAGTACCGCCTCAGTCACCa      >  2:947509/1‑143 (MQ=255)
                                                                                                                                     gCACGCCATCTTACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGCGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACa      <  1:434923/143‑1 (MQ=255)
                                                                                                                                        cGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAc   <  2:46786/143‑1 (MQ=255)
                                                                                                                                         gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc  <  1:614710/143‑1 (MQ=255)
                                                                                                                                         gTCATCTGACGTGCCTTTTATATTAGTACAACCTTGTACGATTATTATTCGCTTTCAACAGTATAGAGTTATGCTTTGTAAGCCGGGTGCATCAGAATCTGTACGGAGAAAATATCGTCAAGTTACGCATCAGTCAACAAAcc  >  2:133814/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                          
CCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACC  >  NC_000913/4366713‑4366992

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