Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 762,870 A→T 11.5% intergenic (+131/‑144) sucB → / → sucC dihydrolipoyltranssuccinase/succinyl‑CoA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913762,8700AT11.5% 185.0 / 15.3 61intergenic (+131/‑144)sucB/sucCdihydrolipoyltranssuccinase/succinyl‑CoA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base T (2/5);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

GACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATGGACAGAACACAT  >  NC_000913/762731‑763015
                                                                                                                                           |                                                                                                                                                 
gACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATGCGGCGCTTGcgccttat                                                                                                                                          <  2:188924/149‑5 (MQ=38)
gACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATGCGGCGCTTGcgccttat                                                                                                                                          <  3:262394/149‑5 (MQ=38)
gACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCAt                                                                                                                                          >  4:321536/1‑149 (MQ=255)
 aCGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATc                                                                                                                                         <  3:148266/149‑1 (MQ=255)
 aCGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATc                                                                                                                                         <  1:465807/149‑1 (MQ=255)
     gtagtaGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATGCGGCGCTTGCGCCTTATCCGAc                                                                                                                                     <  2:360676/149‑1 (MQ=38)
           gtTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATGCGGCGCTTGCGCCTTATCCGAcct                                                                                                                                   <  4:142851/145‑3 (MQ=38)
           gtTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATGCGGCGCTTGCGCCTTATCCGAcct                                                                                                                                   >  3:142851/1‑143 (MQ=38)
            tttAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAgcgc                                                                                                                                               <  1:219003/132‑1 (MQ=255)
            tttAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAgcgc                                                                                                                                               >  2:219003/1‑132 (MQ=255)
            tttAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGAt                                                                                                                              >  1:58448/1‑149 (MQ=255)
                aGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCt                                                                                                                          <  3:64441/149‑1 (MQ=255)
                  tttCACCTGTACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCGAGCGCCGCATCCGACAAGCGATGCCTgc                                                                                                                        >  4:409281/1‑148 (MQ=255)
                   ttCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGAt                                                                                                                       <  2:29725/149‑1 (MQ=255)
                        cTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATcc                                                                                                                                        >  1:469841/1‑127 (MQ=255)
                        cTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATcc                                                                                                                                        <  2:469841/127‑1 (MQ=255)
                          gCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGt                                                                                                                >  2:75534/1‑149 (MQ=255)
                               gTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAAc                                                                                                           <  2:387793/149‑1 (MQ=255)
                               gTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAAc                                                                                                           <  1:26839/149‑1 (MQ=255)
                                 aGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACgt                                                                                                         <  2:387148/149‑1 (MQ=255)
                                  gACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACgtg                                                                                                        >  2:68900/1‑149 (MQ=255)
                                   aCCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACgtgt                                                                                                       <  4:427678/149‑1 (MQ=255)
                                   aCCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACgtgt                                                                                                       <  3:169320/149‑1 (MQ=255)
                                   aCCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACgtgt                                                                                                       >  3:132249/1‑149 (MQ=255)
                                     cGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCt                                                                                                     <  4:352116/149‑1 (MQ=255)
                                           aGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAg                                                                                               <  3:321536/149‑1 (MQ=255)
                                           aGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAg                                                                                               <  4:132249/149‑1 (MQ=255)
                                             gCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGc                                                                                             >  2:188456/1‑149 (MQ=255)
                                               aTTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACaa                                                                                                                                   >  3:58782/1‑109 (MQ=255)
                                               aTTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACaa                                                                                                                                   <  4:58782/109‑1 (MQ=255)
                                                           tcCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGa                                                                               <  2:444556/149‑1 (MQ=255)
                                                           tcCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGa                                                                               >  4:262398/1‑149 (MQ=255)
                                                            ccGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGAc                                                                              <  4:221343/149‑1 (MQ=255)
                                                             cGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGCCAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACa                                                                             >  3:452425/1‑149 (MQ=255)
                                                              ggCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACaa                                                                            <  1:75534/149‑1 (MQ=255)
                                                                   ttGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGccc                                                                       <  2:58448/149‑1 (MQ=255)
                                                                    tGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGc                                                                                             <  2:69288/126‑1 (MQ=255)
                                                                    tGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGc                                                                                             >  1:69288/1‑126 (MQ=255)
                                                                           tGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGa                                                               <  3:262398/149‑1 (MQ=255)
                                                                                 gACGCTGACGCGTCTTATCGGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGaa                                                         >  4:186053/1‑149 (MQ=255)
                                                                                              cTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTAc                                            <  1:188456/149‑1 (MQ=255)
                                                                                               ttATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACg                                           >  3:305771/1‑149 (MQ=255)
                                                                                                  tCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGtt                                        >  1:314888/1‑149 (MQ=255)
                                                                                                   cAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGttt                                       <  1:68900/149‑1 (MQ=255)
                                                                                                       ccTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTaaaa                                   <  2:372364/149‑1 (MQ=255)
                                                                                                              gACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACgtgt                                                                                                       >  1:83338/1‑74 (MQ=255)
                                                                                                              gACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACgtgt                                                                                                       <  2:83338/74‑1 (MQ=255)
                                                                                                                       tGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCTACAATGCCCGGACGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGaa                   <  4:452425/149‑1 (MQ=255)
                                                                                                                          ttgTCGGATGCGGCGCTTGCGCCTTATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTaaaa                                   >  1:2713/3‑130 (MQ=38)
                                                                                                                          ttgTCGGATGCGGCGCTTGCGCCTTATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTaaaa                                   <  2:2713/128‑1 (MQ=38)
                                                                                                                            gTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATg              >  4:87595/1‑149 (MQ=255)
                                                                                                                            gTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATg              >  4:102995/1‑149 (MQ=255)
                                                                                                                               ggATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATGGAc           <  4:305771/149‑1 (MQ=255)
                                                                                                                                gATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTc                                                   >  2:156530/1‑108 (MQ=255)
                                                                                                                                gATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTc                                                   <  1:156530/108‑1 (MQ=255)
                                                                                                                                gATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATGGACa          >  1:190950/1‑149 (MQ=255)
                                                                                                                                  tAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGaa                   <  2:111419/138‑1 (MQ=255)
                                                                                                                                  tAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGaa                   >  1:111419/1‑138 (MQ=255)
                                                                                                                                   aaGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATGGACAGaa       <  3:186053/149‑1 (MQ=255)
                                                                                                                                      gcgcAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATGGACAGAacac    >  3:277901/1‑149 (MQ=255)
                                                                                                                                        gcAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATGGACAGAACACAt  >  3:62197/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCACCAATGTAGGTCGGATAAGGCGCAAGCGCCGCATCCGACAAGCGATGCCTGATGTGACGTTTAACGTGTCTTATCAGGCCTACGGGTGACCGACAATGCCCGGAAGCGATACGAAATATTCGGTCTACGGTTTAAAAGATAACGATTACTGAAGGATGGACAGAACACAT  >  NC_000913/762731‑763015

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