Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,776,2630TG43.1% 96.6 / 15.7 58V226G (GTG→GGG) ydiFputative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base T (5/28);  new base G (25/0);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 8.35e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTTG  >  NC_000913/1776116‑1776412
                                                                                                                                                   |                                                                                                                                                         
aaGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCa                                                                                                                                                        <  1:651217/151‑1 (MQ=255)
  gCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTTGATGATGCAGGGGccg                                                                                                                                                       >  2:333112/1‑148 (MQ=255)
    gATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGaa                                                                                                                                                    >  2:831644/1‑150 (MQ=255)
        gcgcCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATg                                                                                                                                                <  2:681070/150‑1 (MQ=255)
             aGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGTTGATTCAGGGGCAGAAAATgggt                                                                                                                                             >  2:588597/1‑146 (MQ=255)
                    gCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGcc                                                                                                                                    <  2:697840/150‑1 (MQ=255)
                    gCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCa                                                                                                                                   <  2:957012/151‑1 (MQ=255)
                     ccTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCAc                                                                                                                                  <  2:553774/151‑1 (MQ=255)
                         cATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCt                                                                                                                               <  1:395608/150‑1 (MQ=255)
                          aTTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTTGATGGTGCAGGGGGAGAAAATGGGTaaaaaa                                                                                                                                       >  1:717467/1‑141 (MQ=255)
                              gcgcTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGTGATGATGCAGGGGGAGAAAAAGGGTTAAAAAACCACGCCGCCTcc                                                                                                                         >  1:264694/1‑151 (MQ=255)
                                gcTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTa                                                                                                                       <  1:204894/151‑1 (MQ=255)
                                  tACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGTTGATGCAGGGGGAGAAAAAGGGTAAAAAAAACACGCCGGATCCCaa                                                                                                                      >  1:516859/1‑150 (MQ=255)
                                            gACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGGTGCAGGGGCAGAAAAAGGGTAAGAAAACCACGCCGCATACCAAAACTTTCCGt                                                                                                            >  2:685919/1‑150 (MQ=255)
                                             aCAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGAAAATGGGTAAGAAAACCACGCTGCATCCTAAATCTGTcc                                                                                                              >  2:749566/1‑147 (MQ=255)
                                                gTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTcc                                                                                                      <  1:264981/151‑1 (MQ=255)
                                                 tGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGCAGGGGCAGGAAATGGGTAAGAAAACCACGCTGCATccgatccg                                                                                                                   >  1:906588/1‑132 (MQ=255)
                                                 tGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATc                                                                                                                          <  2:906603/131‑1 (MQ=255)
                                                     ggCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATTATTGAGGGGCCGGAAATGGGTAAAAAAAACACGCTGCCTCCTAAAATTGTTCGTGTTTCggg                                                                                                   >  1:734427/1‑149 (MQ=255)
                                                        tACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTGGATGAAGGAAGGGCAGGAAAAGGGTTAGAAAAGCACGGCGCCTCCCAAATCTTTCCCGATTTCggggttt                                                                                               >  1:747203/1‑146 (MQ=255)
                                                           gCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCCGAAAATGGGTAAGAAAGACACGCTGCATCCTTAATCTGTTCGTATTCCGGGGTATCtg                                                                                            >  1:478235/1‑150 (MQ=255)
                                                           gCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAAAGGGTAAGAAAACCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCt                                                                                             >  2:741793/1‑149 (MQ=255)
                                                            ccACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGCAGCAGGGGCAGGAAATGGTTAAGAAAGCCACCCTGCATCCTAAAACTTTCCGGATTCCggggtttc                                                                                              >  2:763633/1‑142 (MQ=255)
                                                                  ttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGata                                                                                     <  2:608188/150‑1 (MQ=255)
                                                                   tGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTTATGATGGAGGGGCAGAAAAAGGTTAAGAAAGACACGCTGCATCCTAAAACTTTCCCTGTTTCGGGTTATTtggggg                                                                                        >  1:547730/1‑146 (MQ=255)
                                                                   tGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGGTTAGAAAACCACGCCGCCTCCTAAATCTGTCCGTCTTTCGGGGTTTCTGGTGGGTAtt                                                                                   >  1:889565/1‑151 (MQ=255)
                                                                    gAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATGCAGGGGCCGAAAAAGGGTAAGAAAGCCACGCTGCATCCTAAATCTTTCCGTATTCCCGGGTTTCTGGGGGGTAtt                                                                                   >  2:236526/1‑150 (MQ=255)
                                                                    gAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATtg                                                                                  <  1:17921/151‑1 (MQ=255)
                                                                     aaGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGCATTTTGATGATTCAGGTGCAGAAAATGGGTAAGAAAACCACGCTTCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATAtttt                                                                                 >  2:223170/1‑149 (MQ=255)
                                                                       gATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGGTTGAGGGGCAGAAAATGGTTAAAAAAGACACGCCGCATCCTAAAACTTTTCCTATTTCGGGGTTTTTTGTGGGTAtttt                                                                                 >  1:507886/1‑147 (MQ=255)
                                                                       gATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTgtgg                                                                               <  2:626572/151‑1 (MQ=255)
                                                                               gATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGa                                                                        <  2:736515/150‑1 (MQ=255)
                                                                                  gTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCg                                                                   <  1:763616/151‑1 (MQ=255)
                                                                                        cGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCa                                                              <  1:801055/150‑1 (MQ=255)
                                                                                         gACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTTTTGTGATGATGGAGGTGCCGAAAAAGGGTAAAAAAGCCACGCTGCATCCTTAATCTTTTCGTATTTCGGGGTATCTGGTGGGTATTTTGGTGGTCGATCCGGCTCaa                                                             >  1:314100/1‑150 (MQ=255)
                                                                                         gACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTGGATGATGGAGGGGGCGAAAAAGGGTAAGAAAACCACGCTGCATCCTAAATCTTTCCGTATTTCGGGGTTTCTTGTGGGTATTGTGGGGGGCGCTTCGGATCaa                                                             >  1:905240/1‑150 (MQ=255)
                                                                                           cGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAAAGGTTAAGAAAACCACGCTGCATCCTAAAACTGTCCGGATTTCGGGGTTTCTGGGGGGTATTGTGGTGGGCGATCCGGATCAAAcc                                                          >  2:660873/1‑151 (MQ=255)
                                                                                              aTTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCaa                                                       <  1:870829/151‑1 (MQ=255)
                                                                                                tGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTTTTGTGATGATGCAGGGGCAGAAAAAGGGTAAGAAAGCCACGCCGCATCCTAAATTTGTCCGTATTCCGGGTTATCTGGTGGGTATTGTGGGGGTCGGTCCCGGTCAAAcccca                                                       >  2:18691/1‑147 (MQ=255)
                                                                                                 ggTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGTGATGATGGAGGGGCAGGAAAAGGGTTAAAAAAGCACGCCGCATCCTAAAACTTTTCGCATTTCCGGGTTTTTGGGGGATATTGTGGGGGGGGGTCCCGGATCAAACCCa                                                        >  1:659794/1‑148 (MQ=255)
                                                                                                   ttATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGGTAAGAAAACCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGGTATTGTGGTGGTCGATCCGGATCAAACCCAACTGt                                                   >  2:833588/1‑150 (MQ=255)
                                                                                                         cccAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCgg                                            <  2:814076/151‑1 (MQ=255)
                                                                                                         cccAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCgg                                            <  2:240440/151‑1 (MQ=255)
                                                                                                                 gTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGAAAAAGGGTAAAAAAAACACGCTGCATTCTAAATCTTTTCGTATTTCGGGGTATCTGGTGGGTATTTTGGTGGGCGATCCGGGTCAAACCCCACTTTTTGGCGGGGGAcccg                                   >  1:943837/1‑149 (MQ=255)
                                                                                                                   gCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGAAAAAGGTTAAAAAAGCCACGCTGCATCCTAAATCTGTCCGTATTTCGGGGTATCTGGTGGGTATTGTGGTGGTCGATCCGGGTCAAACCCCACTGTATGGCGGGGGACCgg                                   >  1:103832/1‑149 (MQ=255)
                                                                                                                   gCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGt                                  <  2:951529/150‑1 (MQ=255)
                                                                                                                        aTAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCg                            <  2:406423/151‑1 (MQ=255)
                                                                                                                        aTAACGGCGGTATTGTGATGATGCAGGGGGAGAAAAAGGTTAAGAAAGCCACCCTGCCTCCTAAATCTGTCCGTATTTCGGGGTTTTTGGTGGGTATTGTTGTGGTCGATCCGGATCAAACCCCACCGTATGGGGGGGCACCCGGTTAAc                              >  1:818267/1‑150 (MQ=255)
                                                                                                                          aaCGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCt                          <  2:478245/151‑1 (MQ=255)
                                                                                                                           aCGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCtt                         <  1:714123/151‑1 (MQ=255)
                                                                                                                               cggTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTAtt                     <  2:30965/151‑1 (MQ=255)
                                                                                                                               cggTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTAtt                     <  1:833572/151‑1 (MQ=255)
                                                                                                                                ggTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTAttt                    >  2:215788/1‑151 (MQ=255)
                                                                                                                                        gatgatGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGAc            <  2:473209/151‑1 (MQ=255)
                                                                                                                                           gatgCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACtt          <  2:516870/150‑1 (MQ=255)
                                                                                                                                            atgCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTc         >  2:821560/1‑150 (MQ=255)
                                                                                                                                                  gTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTtg  <  2:620834/151‑1 (MQ=255)
                                                                                                                                                  gTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTtg  <  1:18692/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
AAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTTG  >  NC_000913/1776116‑1776412

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