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

CGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCT  >  NC_000913/1776119‑1776410
                                                                                                                                                |                                                                                                                                                     
cGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGaa                                                                                                                                                <  1:93735/151‑1 (MQ=255)
cGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAg                                                                                                                                                   >  2:1026627/1‑149 (MQ=255)
 gATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGGAGGGGccg                                                                                                                                                   >  2:1165697/1‑146 (MQ=255)
       gcCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCCGAAAATgggt                                                                                                                                         >  2:788429/1‑149 (MQ=255)
        ccAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTa                                                                                                                                        <  2:70410/151‑1 (MQ=255)
         cAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGGTGCAGGGGCAGAAAAAgggt                                                                                                                                         >  1:236382/1‑147 (MQ=255)
                 gCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCa                                                                                                                               <  1:682249/151‑1 (MQ=255)
                     tCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATTATGCAGGGGCAGAAAAAGGGTAAGAAAACCAc                                                                                                                              >  1:566741/1‑148 (MQ=255)
                     tCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACg                                                                                                                             <  1:579507/149‑1 (MQ=255)
                       aTTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGc                                                                                                                         <  2:604776/151‑1 (MQ=255)
                          cgcgCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCAt                                                                                                                       <  2:461174/150‑1 (MQ=255)
                          cgcgCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAAAGGTTAAGAAAGCCACGCTGCAt                                                                                                                       >  2:423383/1‑150 (MQ=255)
                            cgcTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCt                                                                                                                    <  1:218366/151‑1 (MQ=255)
                                accaccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTaa                                                                                                                  <  2:160073/149‑1 (MQ=255)
                                 ccaccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATc                                                                                                               <  1:138156/151‑1 (MQ=255)
                                   accTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGTTGATGCAGGGGCAGGAAAAGGGTAAGAAAAACACGCTGCCTCCTAAAtat                                                                                                              >  1:912280/1‑148 (MQ=255)
                                     cTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGGAGGGGCAGAAAATGGTTAAGAAAGACACGCTGCATCCTAAATCTGt                                                                                                            >  2:168748/1‑150 (MQ=255)
                                     cTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGCAGGGGCAGAAAATGGGTAAGAAAACCACCCTGCCTCCTAAATCTGt                                                                                                            >  1:161255/1‑150 (MQ=255)
                                       gCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTcc                                                                                                          <  2:912286/150‑1 (MQ=255)
                                           cAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTAtt                                                                                                     <  1:467558/151‑1 (MQ=255)
                                                aaGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCggg                                                                                                <  1:315082/151‑1 (MQ=255)
                                                  ggCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGtt                                                                                              <  1:349436/151‑1 (MQ=255)
                                                      aCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTAt                                                                                            <  2:353694/149‑1 (MQ=255)
                                                      aCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATc                                                                                           <  2:1129774/150‑1 (MQ=255)
                                                             ttttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGGTGGAGGGGGAGAAAAAGGTTAAAAAAAGCACGCTGCATCCTAAATCTTTCCGTATTCCGGGGTATCttgtggg                                                                                    >  1:250827/1‑149 (MQ=255)
                                                              tttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGat                                                                                   <  1:1140388/150‑1 (MQ=255)
                                                                tGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGCTGGTGAAAGGGCAGAAAATAGTTAAAAAAACCACACGGCATCCCAAATCTGTTCGTATTCCGGGTTTTCttgcggggatt                                                                                >  2:217776/1‑141 (MQ=255)
                                                                  aaGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAAGGGGAGAAAAAGGTTAAAAAAGACACGCCTCATCCTCAAAATGTCCCGATTCCGGGGTATTTGGGGGGTAtt                                                                                >  2:900336/1‑149 (MQ=255)
                                                                     aTGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATGGAGGGGCAGAAAAAGGGTAAAAAAGACACGCGGCATTCCAAATATGTTCGTATTCCCGGGTATTTGGTGGATAAtggtg                                                                            >  1:430682/1‑147 (MQ=255)
                                                                       gAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGAAAATGGGTAAAAAAACCACCCCGCCTCCTAAATTTGTCCGTATTTCGGGGTTTCTGGGGGGTATTTtggtgg                                                                         >  2:243478/1‑151 (MQ=255)
                                                                       gAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTgtggtg                                                                          <  2:258458/150‑1 (MQ=255)
                                                                          gTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGGGGTTGTGGAGGTGCAGCAAAAGGGTTAAAAAGCCCCGCCGCATCCTAAATCTGTTCCTGTTTCGGGTTATCTGGTGGGTAtt                                                                                >  1:134317/1‑141 (MQ=255)
                                                                                 aTCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCgg                                                                <  1:42059/150‑1 (MQ=255)
                                                                                       aCGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGGTGCAGGGGCAGAAAAAGGGTAAGAAAACCACCCTGCATCCTAAATCTTTTCGTATTCCGGGGTTTTTTGTGGGTATTGTGGTGGGCGATcccggtccaa                                                          >  1:1039252/1‑142 (MQ=255)
                                                                                        cGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAAcc                                                        <  1:1048458/151‑1 (MQ=255)
                                                                                         gCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGAAAAAGGGTAAGAAAACCACGCTGCATCCTAAAACTGTTCGTATTCCGGGGTATCTGGGGGGTATTGTGGGGGGCGATCCGGGTCAAAacc                                                       >  2:1114096/1‑151 (MQ=255)
                                                                                            ttGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGGTGCAGGGGCAGAAAATGGGTGAGAAAACCACGCTGCCTCCTAAATCTGTCCGTTTTTCCGGGTATCTGGGGGGAATTGTGGTGGGCGCTCCGGGTCAAACCCaa                                                     >  1:288243/1‑150 (MQ=255)
                                                                                                  aTTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAACCACGCTGCATCCTAAATCTGTCCGTATTCCGGGGTATCTTGTGGGTATTGTGGTGGTCGATCCGGGTCAAACCCAACTGTa                                                >  2:1042185/1‑149 (MQ=255)
                                                                                                    tGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTGAATCTGTCCGTATTCCGGGTTATCTGTTGTATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATgg                                             <  2:153486/150‑1 (MQ=255)
                                                                                                     gCCCAGGCGGTGCACAATAACGGCGGCATTGTGATGATGCAGGGGCCGAAAAAGGGTAAGAAAACCACGCCGCCTCCTAAATCTTTTCGTATTCCGGGTTTTTTGGGGGGTATTTTGGGGGGTGATCCCGGTCAAACCCAACTTTATgggg                                           >  1:255955/1‑149 (MQ=255)
                                                                                                      cccAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGGTTAGAAAACCCCCCCGCCTCCTAAATCTGTCCGTATTTCGGGGTTTCTGGGGGGTATTTTGGGGGTCGATCCCGCTCAAAACCAACAGTATTGCgg                                          >  1:1172064/1‑151 (MQ=255)
                                                                                                      cccAGGCGGTGCACAATAACGGCGGAATTTTGATGATGCAGGGGGAGAAAAAGGTTAAGAAAACCACGCCGCCTCCTAAATCTGTCCGTATTTCGGGGTTTCTTGTGGGTATTTTGGGGGGCGATCCGGATCAAACCCCACTGTATGGCgg                                          >  2:1135037/1‑151 (MQ=255)
                                                                                                          ggCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCCGAAAATGGGTAAGAAAACCACCCCGCATCCTAAATCTGTCCGGATTTCGGGGTATCTGGTGGGTATTGTGGTGGGCGATCCGGGTCAAACCCCAATGTATGgggggg                                        >  1:257557/1‑147 (MQ=255)
                                                                                                            cGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCAcc                                   <  2:105385/151‑1 (MQ=255)
                                                                                                                gCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGt                                <  1:271084/150‑1 (MQ=255)
                                                                                                                 cacaATAACGGCGGTATTGTGGTGATGCAGGGGCAGAAAAAGGGTAAAAAAACCCCGCTGCATCCTAAATCTGGCCGGATTCCGGGGTATTTGGGGGGGATTTTGGGGGTCGATCCCGGGCAAACCCAAATGTATGGCGGTGCAcccggt                               >  1:314226/1‑146 (MQ=255)
                                                                                                                   caATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAAAGGTTAAGAAAACCCCCCCGCATCCTAAAACTGTCCGTATTCCGGGGTATCTGGGGGGTATTGTGGTGGTCGATCCGGATCAAACCCCACCTTATGGCGGGGCACCCGGTTa                              >  1:337245/1‑150 (MQ=255)
                                                                                                                        aCGGCGGTATTGTGATGATGCAAGGGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGGATTTCGGGTTATCTGGGGGGTATTGTGGTGGGCGATCCGGATCAAACCCAACTGTATGGCGGGGCCCCGGTTAAccccct                       >  1:431588/1‑147 (MQ=255)
                                                                                                                          ggcggTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCttt                      >  2:952140/1‑150 (MQ=255)
                                                                                                                               tATTGAGATGATGAAGGAGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGAGGCCGATCCGGATCAAACCCAAATGTATGGCGGCGGAACGGTTAACCGCTTTATTTCt                >  1:900217/1‑151 (MQ=255)
                                                                                                                                aTTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTg               <  1:452944/151‑1 (MQ=255)
                                                                                                                                  tgtgATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTgg              <  1:1152553/150‑1 (MQ=255)
                                                                                                                                      atgatgCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACt         <  2:1194944/151‑1 (MQ=255)
                                                                                                                                      atgatgCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACt         >  1:1069857/1‑151 (MQ=255)
                                                                                                                                         atgCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTc       >  2:434657/1‑150 (MQ=255)
                                                                                                                                           gCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCAcc    >  1:346932/1‑151 (MQ=255)
                                                                                                                                            cAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCAccc   >  1:514486/1‑151 (MQ=255)
                                                                                                                                             aGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCt  <  1:578679/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCT  >  NC_000913/1776119‑1776410

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