Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,314,261 C→T 6.5% P81L (CCG→CTG)  ompW → outer membrane protein W

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,314,2610CT 6.5% 189.9 / 3.0 61P81L (CCG→CTG) ompWouter membrane protein W
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base T (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GTTCTGCAACCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCAGGTATTAACTACACCAC  >  NC_000913/1314114‑1314405
                                                                                                                                                   |                                                                                                                                                 
gTTCTGCAACCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGcc                                                                                                                                                   <  6:236561/148‑1 (MQ=255)
gTTCTGCAACCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGcc                                                                                                                                                   <  2:243245/148‑1 (MQ=255)
     gCAACCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTcc                                                                                                                                              >  7:15094/1‑148 (MQ=255)
     gCAACCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTcc                                                                                                                                              >  3:226070/1‑148 (MQ=255)
      cAACCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGtt                                                                                                                                                <  7:119359/145‑1 (MQ=255)
        aCCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCg                                                                                                                                             >  1:297989/1‑146 (MQ=255)
        aCCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCg                                                                                                                                             <  2:297989/146‑1 (MQ=255)
             aCGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATaa                                                                                                                                       >  5:427866/1‑147 (MQ=255)
                  cAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCg                                                                                                                                  <  4:144838/147‑1 (MQ=255)
                     cAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGCATTACTGGCAGCGACGCTGTTCCGCCATAAAATCGGCAc                                                                                                                              >  2:84394/1‑148 (MQ=255)
                     cAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCTGTTCCGCCATAAAATCGGCAc                                                                                                                              >  2:90398/1‑148 (MQ=255)
                         aGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCgg                                                                                                                          <  2:34839/148‑1 (MQ=255)
                              cTGGGGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGa                                                                                                                      <  5:72048/147‑1 (MQ=255)
                                ggtggtACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCg                                                                                                                   >  5:3781/1‑148 (MQ=255)
                                ggtggtACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCATTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCg                                                                                                                   >  7:421823/1‑148 (MQ=255)
                                     tACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGat                                                                                                              <  8:421823/148‑1 (MQ=255)
                                       cGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGatat                                                                                                            <  6:63904/148‑1 (MQ=255)
                                         ttAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTg                                                                                                          <  2:221489/148‑1 (MQ=255)
                                           aGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCTACAACATTGGAGTGGAATTACTGGCAGCGACGCCGTTCAGACATAAAATCGGCACACGGGAGACCGGCGATATTGc                                                                                                         >  4:176442/1‑147 (MQ=255)
                                                 tCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGc                                                                                                                                >  6:313419/1‑118 (MQ=255)
                                                 tCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGc                                                                                                                                <  5:313419/118‑1 (MQ=255)
                                                 tCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGc                                                                                                                                <  7:286407/118‑1 (MQ=255)
                                                 tCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGc                                                                                                                                >  8:286407/1‑118 (MQ=255)
                                                 tCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGtt                                                                                                  <  6:427866/148‑1 (MQ=255)
                                                  cTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGtt                                                                                                  >  3:393375/1‑147 (MQ=255)
                                                    gggTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGat                                                                                                              <  1:122998/133‑1 (MQ=255)
                                                    gggTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGat                                                                                                              >  2:122998/1‑133 (MQ=255)
                                                     ggTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTtcat                                                                                               <  5:257314/147‑1 (MQ=255)
                                                        ggATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCTGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCatcatc                                                                                           <  1:90398/148‑1 (MQ=255)
                                                         gATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGGGGAATTACTGGCAGCGACGCTGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCatcatc                                                                                           <  1:84394/147‑1 (MQ=255)
                                                            tCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGat                                                                                                              >  2:89828/1‑125 (MQ=255)
                                                            tCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGat                                                                                                              <  1:89828/125‑1 (MQ=255)
                                                                cGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTAATCATCTGccacc                                                                                    >  3:329769/1‑147 (MQ=255)
                                                                       aaTAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGa                                                                            >  8:143016/1‑148 (MQ=255)
                                                                       aaTAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGa                                                                            <  8:15094/148‑1 (MQ=255)
                                                                       aaTAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGGGGAATTACTGGCAGCGACGCCGTTCCGACATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGa                                                                            <  4:226070/148‑1 (MQ=255)
                                                                                cAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGt                                                                  >  8:301271/1‑148 (MQ=255)
                                                                                           tACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCa                                                                                                                               >  7:363251/1‑77 (MQ=255)
                                                                                           tACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCa                                                                                                                               <  8:363251/77‑1 (MQ=255)
                                                                                           tACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGt                                                       >  8:73690/1‑148 (MQ=255)
                                                                                            aCGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTgg                                                        <  7:143016/146‑1 (MQ=255)
                                                                                            aCGTTTACTTATATGGCGACCGACAACATTGGGGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTgg                                                        <  3:176442/146‑1 (MQ=255)
                                                                                             cGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTg                                                      <  7:73690/147‑1 (MQ=255)
                                                                                               tttACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTg                                                      <  4:367467/145‑1 (MQ=255)
                                                                                               tttACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTg                                                      >  3:367467/1‑145 (MQ=255)
                                                                                                ttACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGa                                                                                                                      >  1:220336/1‑81 (MQ=255)
                                                                                                ttACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGa                                                                                                                      <  2:220336/81‑1 (MQ=255)
                                                                                                     tataTGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGccacc                                                                                    >  1:181250/1‑110 (MQ=255)
                                                                                                     tataTGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGccacc                                                                                    <  2:181250/110‑1 (MQ=255)
                                                                                                      atatGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCcagc                                              <  7:154535/146‑1 (MQ=255)
                                                                                                                  acaacaTTCGTCTCGAAATACTGGCAGCCACGCCCATCCGCCATAAAAACGGCACCCGGGCGACCGGCGATATTGAAACCCTTCATCATCTGGCACCCACACTGATGCGCAACTGGTATTGAGCACATCCCAGGAGCAAAttccttcc                                >  6:353039/1‑148 (MQ=255)
                                                                                                                   caacaTTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCt                               >  2:224748/1‑148 (MQ=255)
                                                                                                                   caacaTTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCt                               >  3:320197/1‑148 (MQ=255)
                                                                                                                         tGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGtt                         >  5:77116/1‑148 (MQ=255)
                                                                                                                           gtgtGGAATTACTGGCAGCGACGCCGGTCCGCCATAAAATCGGCACCCGGGCGACCGTCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTg                        <  4:77480/147‑1 (MQ=255)
                                                                                                                               ggAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCACCCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTAGTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTGgggg                     >  7:418574/1‑146 (MQ=255)
                                                                                                                                gAATAACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCa                   <  5:212429/147‑1 (MQ=255)
                                                                                                                                        ggCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCAGGTATTaa           >  8:280028/1‑147 (MQ=255)
                                                                                                                                        ggCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCAGGTATTAAc          >  8:143688/1‑148 (MQ=255)
                                                                                                                                             cGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCAGGTATTAACTacac     >  7:38681/1‑148 (MQ=255)
                                                                                                                                               aCGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCAGGTATTAACTAcacca   >  3:398037/1‑148 (MQ=255)
                                                                                                                                               aCGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCAGGTATTAACTAcacca   >  7:128514/1‑148 (MQ=255)
                                                                                                                                                cGCCGTTCCGCCATAAAATTGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAATTCCGTCCTTACGTTGGGGCAGGTATTAACTAcaccac  >  6:332622/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTTCTGCAACCGTACGTCCAACAGAAGGTGCTGGTGGTACGTTAGGAAGTCTGGGTGGATTCAGCGTGACCAATAACACGCAACTGGGCCTTACGTTTACTTATATGGCGACCGACAACATTGGTGTGGAATTACTGGCAGCGACGCCGTTCCGCCATAAAATCGGCACCCGGGCGACCGGCGATATTGCAACCGTTCATCATCTGCCACCAACACTGATGGCGCAGTGGTATTTTGGTGATGCCAGCAGCAAATTCCGTCCTTACGTTGGGGCAGGTATTAACTACACCAC  >  NC_000913/1314114‑1314405

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