Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,026,252 2 bp→GC coding (110‑111/1164 nt) fadA ← thiolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,026,2520CG99.2% 331.7 / ‑6.0 126T37S (ACG→AGCfadAthiolase I
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (59/66);  total (60/66)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
*NC_0009134,026,2530GC99.2% 326.7 / NA 129T37S (ACG→AGC) fadAthiolase I
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base C (58/69);  minor base A (0/1);  total (59/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCA  >  NC_000913/4026181‑4026322
                                                                       ||                                                                     
gTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCCGCGCCGGGGTACGCCCCCGCAGGGtgc                                                                    >  2:731168/1‑74 (MQ=255)
  ctgctgCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGc                                                                  >  1:3130748/1‑76 (MQ=255)
   tgctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCa                                                                 <  1:5369958/76‑1 (MQ=255)
    ggtgcCCACAACCCCAGTAAATGTCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAt                                                                <  2:4966665/74‑1 (MQ=255)
    gctgcACACACCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAt                                                                <  1:1323156/76‑1 (MQ=255)
    gctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTGCGCAt                                                                >  2:4656917/1‑76 (MQ=255)
     ctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAtt                                                               <  1:2039934/76‑1 (MQ=255)
     cggcACACACCCCCAGGAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAtt                                                               <  2:2640678/74‑1 (MQ=255)
       ccccacaCCCCCAGTAAATACCGTCGGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaa                                                             <  2:653041/73‑1 (MQ=255)
        cacacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCCGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaaa                                                            >  2:1847275/1‑76 (MQ=255)
         ccacaCCCCCAGAAAATATCGGCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAAt                                                           <  2:5714109/75‑1 (MQ=255)
         acacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCCGCAGGCTACGCATTAAAt                                                           >  2:2523242/1‑76 (MQ=255)
         acacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAAt                                                           >  1:2797388/1‑76 (MQ=255)
          cacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATg                                                          >  1:3476249/1‑76 (MQ=255)
            caACCCCAGTAAATATCGTCGAGGGCCGCCGCTGCCAGCGCCGGGGTACGCGCCAGCAGGCCACGCATTCAATGac                                                        >  2:5687047/1‑75 (MQ=255)
             accccgAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTACATGAGc                                                       <  2:4462284/74‑1 (MQ=255)
             acccccAGTAAATATCGTCGGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGc                                                       <  1:1384545/74‑1 (MQ=255)
              aCCCCAGTAAATACCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                      <  1:4807949/76‑1 (MQ=255)
               ccccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCgg                                                     <  2:2542718/76‑1 (MQ=255)
                cccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGga                                                    <  2:5176650/76‑1 (MQ=255)
                 ccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgag                                                   >  1:411637/1‑76 (MQ=255)
                  cAGTAACTATCGGCTGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                  <  2:4438167/76‑1 (MQ=255)
                  cAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                  <  1:4853526/76‑1 (MQ=255)
                  cAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                  >  1:3827009/1‑76 (MQ=255)
                  cAGTAAATACCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                  <  2:3504358/76‑1 (MQ=255)
                   aGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagag                                                 >  1:740792/1‑76 (MQ=255)
                    gaaaaTAGCGTCGGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagaga                                                <  2:3524671/74‑1 (MQ=255)
                    gTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagaga                                                >  2:4649198/1‑76 (MQ=255)
                     tcaaTATCGTCGGGGGCCGCCTCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGatat                                               <  2:4045803/74‑1 (MQ=255)
                     tAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                               <  1:2041969/76‑1 (MQ=255)
                      aaaTATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtc                                              <  1:5169541/76‑1 (MQ=255)
                      aaaTATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtc                                              <  1:5578216/76‑1 (MQ=255)
                      aaaTATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtc                                              >  2:2269274/1‑76 (MQ=255)
                        atgtcggcgGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtctt                                            <  2:253045/73‑1 (MQ=255)
                        atatCGTGGAGGGCCGCCGCTTCCCGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtctt                                            <  2:542494/76‑1 (MQ=255)
                        atatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtctt                                            >  1:1257104/1‑76 (MQ=255)
                         tatCGCCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttc                                           <  1:3257192/76‑1 (MQ=255)
                          atCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttct                                          >  1:5363047/1‑76 (MQ=255)
                           tcgtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTg                                         <  1:821841/76‑1 (MQ=255)
                            cggcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgc                                        <  2:4032123/73‑1 (MQ=255)
                             gtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCCGCAGGCTACGCGTTAAATGAGCGGAGAGCTCCTCCgcc                                       >  2:385122/1‑75 (MQ=255)
                             gtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgca                                       >  1:4822639/1‑76 (MQ=255)
                              tcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGCTCTTCTgcac                                      >  2:5411519/1‑76 (MQ=255)
                               cgggggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                     <  2:345791/73‑1 (MQ=255)
                               cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                     <  2:4248331/76‑1 (MQ=255)
                               cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                     <  1:1231289/76‑1 (MQ=255)
                                gggggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                    <  2:4107531/74‑1 (MQ=255)
                                gacggCCGCCGCTTCCATCGCCGGGTTACGCGCCAGCAGGCTACGCATCAACCGAGCGGAGCGCACCTCTgtacgc                                    >  2:1039282/4‑76 (MQ=255)
                                gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                    <  2:2931493/76‑1 (MQ=255)
                                gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                    >  2:743084/1‑76 (MQ=255)
                                gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                    >  1:2272263/1‑76 (MQ=255)
                                gAGGGCCGCCGCTTCCAGCGCCAGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                    >  1:1911156/1‑76 (MQ=255)
                                    gccgccGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGt                                <  2:219889/76‑1 (MQ=255)
                                    gccgccGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGt                                >  1:1547983/1‑76 (MQ=255)
                                    gccgccGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGt                                <  1:3682618/76‑1 (MQ=255)
                                     ccgccgCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGtt                               <  2:2880332/76‑1 (MQ=255)
                                       gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTAc                             <  2:1346186/76‑1 (MQ=255)
                                        ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACg                            >  1:4221906/1‑76 (MQ=255)
                                        ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACg                            <  2:1883438/76‑1 (MQ=255)
                                        ccgcTGCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACg                            <  2:5701819/76‑1 (MQ=255)
                                         cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGa                           >  2:1013545/1‑76 (MQ=255)
                                          gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAGATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                          >  1:2158426/1‑76 (MQ=255)
                                          gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                          >  1:447073/1‑76 (MQ=255)
                                          gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                          >  1:4540753/1‑76 (MQ=255)
                                          gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGCTACGaa                          >  1:1341801/1‑76 (MQ=255)
                                          gcTTCCAGAGCCGGGTTACGCGAAAGCAGGCGACGCATTAAATGAGCGGAGAGATCTGCTCCACGCACGTCACGaa                          >  2:5390035/1‑76 (MQ=255)
                                            ttCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaaaa                        >  1:272425/1‑76 (MQ=255)
                                              ccAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGc                      <  2:4117178/76‑1 (MQ=255)
                                              ccAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGc                      <  2:3060861/76‑1 (MQ=255)
                                                aGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAc                    <  2:1745474/76‑1 (MQ=255)
                                                 gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAcc                   >  2:4940102/1‑76 (MQ=255)
                                                 gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAcc                   <  1:180413/76‑1 (MQ=255)
                                                 gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAcc                   <  1:4592726/76‑1 (MQ=255)
                                                 gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAcc                   >  1:3325880/1‑76 (MQ=255)
                                                 gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAcc                   >  1:3785358/1‑76 (MQ=255)
                                                  cgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCg                  <  1:940522/76‑1 (MQ=255)
                                                  cgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCg                  <  2:5410116/76‑1 (MQ=255)
                                                      gggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt              >  1:3433401/1‑76 (MQ=255)
                                                      gggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt              >  2:2889005/1‑76 (MQ=255)
                                                      gggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTGCTGCACGCACGTTAAGAAAAGCAACGCCCt              >  1:5751933/1‑76 (MQ=255)
                                                      gggTTACGCGCCAGCAGGCTACGCATTAAATGAGAGGAGAGATCTGCTGCACGCACGTTACGAAAAGCACCGCCCt              >  2:5173733/1‑76 (MQ=255)
                                                      gggTTACGCGCCAGCAGGCTAAGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTGACGAAAAGCACCGCCCt              >  2:4015339/1‑76 (MQ=255)
                                                        gTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc            <  1:1846121/76‑1 (MQ=255)
                                                        gTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc            >  2:2983201/1‑76 (MQ=255)
                                                         ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCGCCGCCCTTCg           >  2:2930715/1‑76 (MQ=255)
                                                         ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCg           >  1:3232397/1‑76 (MQ=255)
                                                         ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCg           <  2:4448612/76‑1 (MQ=255)
                                                          tACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGCAAAGCACCGCCCTTCGa          >  1:4322196/1‑76 (MQ=255)
                                                          tACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGa          <  1:1828883/76‑1 (MQ=255)
                                                           aCGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGaa         >  2:4771006/1‑76 (MQ=255)
                                                           aCGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGATCTTCGaa         >  1:2931724/1‑76 (MQ=255)
                                                           aCGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATATTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGaa         >  1:1520295/1‑76 (MQ=255)
                                                             gcgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACg       <  2:5034228/76‑1 (MQ=255)
                                                             gcgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACATTACGAAAAGCACCGCCCTTCGAACg       <  2:5230571/76‑1 (MQ=255)
                                                              cgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTCCGACCga      >  2:3362028/1‑75 (MQ=255)
                                                                ccAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGcc    <  1:1984462/76‑1 (MQ=255)
                                                                  agcagtATACGCGTGAAATGAGCGGAGAAGTCTTCTGCACGCACGTTACGAAAAGCGCCGCCCGTCGAACGGCCCa  <  1:296215/76‑1 (MQ=255)
                                                                  agcaggCTACGCATTAAATGAGCTGAGCGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCa  <  1:42564/76‑1 (MQ=255)
                                                                  agcaggCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCa  <  1:279145/76‑1 (MQ=255)
                                                                  agcaggCTACGCATGAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCa  <  2:4737081/76‑1 (MQ=255)
                                                                       ||                                                                     
GTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCA  >  NC_000913/4026181‑4026322

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