Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 428,188 G→T R181L (CGC→CTC)  secD → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913428,1880GT92.5% 255.8 / 6.7 80R181L (CGC→CTC) secDSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (29/45);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGG  >  NC_000913/428089‑428283
                                                                                                   |                                                                                               
aaCAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCcg                                                                                                 >  2:159344/1‑100 (MQ=255)
  cAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCctcg                                                                                               <  2:1534535/100‑1 (MQ=255)
   aaaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGa                                                                                              >  2:577749/1‑100 (MQ=255)
    aaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGAt                                                                                             <  1:160544/100‑1 (MQ=255)
    aaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGAt                                                                                             >  2:901841/1‑100 (MQ=255)
     aaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                            >  2:14394/1‑100 (MQ=255)
     aaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                            <  1:841042/100‑1 (MQ=255)
     aaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                            <  1:881157/100‑1 (MQ=255)
          cGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTaaa                                                                                       >  1:1318328/1‑100 (MQ=255)
            aTAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGc                                                                                     <  2:589140/100‑1 (MQ=255)
               gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                  >  2:219442/1‑100 (MQ=255)
               gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                  >  2:247903/1‑100 (MQ=255)
               gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                  >  2:169150/1‑100 (MQ=255)
               gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                  >  2:162586/1‑100 (MQ=255)
                 cTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtg                                                                                <  2:531321/100‑1 (MQ=255)
                  tACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGtga                                                                               <  2:1373485/100‑1 (MQ=255)
                  gaCGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCCCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGt                                                                                 <  2:1460664/97‑1 (MQ=255)
                    cGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatg                                                                             <  1:14394/100‑1 (MQ=255)
                    cGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatg                                                                             <  2:663914/100‑1 (MQ=255)
                    cGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatg                                                                             >  1:881461/1‑100 (MQ=255)
                     gCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatga                                                                            >  1:573670/1‑100 (MQ=255)
                          gACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCa                                                                       >  1:502859/1‑100 (MQ=255)
                            ccTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAtt                                                                     <  2:208139/100‑1 (MQ=255)
                               gcgcGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCg                                                                  <  2:557782/100‑1 (MQ=255)
                               gcgcGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCg                                                                  <  2:729286/100‑1 (MQ=255)
                                cgctAAAAGGGTATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCg                                                                  <  1:1576584/99‑1 (MQ=255)
                                cgcgAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGt                                                                 >  2:612732/1‑100 (MQ=255)
                                     aaaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTg                                                            <  1:159344/100‑1 (MQ=255)
                                       aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAg                                                          <  1:673108/100‑1 (MQ=255)
                                       aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCGGAg                                                          <  2:114047/100‑1 (MQ=255)
                                       aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCACTGCGTATCTGAg                                                          <  1:1040974/100‑1 (MQ=255)
                                        gggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGTGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGc                                                         <  2:80894/100‑1 (MQ=255)
                                        gggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGc                                                         >  1:658291/1‑100 (MQ=255)
                                         ggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCa                                                        <  1:1146532/100‑1 (MQ=255)
                                          gCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCaa                                                       >  2:1371857/1‑100 (MQ=255)
                                          gCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCaa                                                       >  2:51496/1‑100 (MQ=255)
                                           cATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAg                                                      <  1:1574290/100‑1 (MQ=255)
                                            aTCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgc                                                     >  2:1170442/1‑100 (MQ=255)
                                             tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                    <  1:313145/100‑1 (MQ=255)
                                             tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                    <  1:1447991/100‑1 (MQ=255)
                                             tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                    <  1:304307/100‑1 (MQ=255)
                                                cGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCa                                                 <  2:634738/100‑1 (MQ=255)
                                                cGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCa                                                 <  1:794420/100‑1 (MQ=255)
                                                   ataCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATcc                                              <  2:140773/100‑1 (MQ=255)
                                                      ccACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGa                                           >  1:1334764/1‑100 (MQ=255)
                                                           gTTCGTAAAGAAAACTACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTgg                                      <  1:1553211/100‑1 (MQ=255)
                                                            ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                     <  2:498976/100‑1 (MQ=255)
                                                            ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                     >  1:1121063/1‑100 (MQ=255)
                                                             tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                    <  2:765507/100‑1 (MQ=255)
                                                             tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                    >  1:1059270/1‑100 (MQ=255)
                                                              cGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAACCATTGAGTATCTGAGCAAGCGCCATCCGGACCTGGTGa                                   >  2:1109408/1‑100 (MQ=255)
                                                              cGTAAAGAAAACAACTACGGACTAAGCATCACTTTGCTCGATGCTAAAGATCGTGATGAAGCCATTGAGTATATGAGCAAGCGCCATCCGGACCTGGTGa                                   >  2:1076360/1‑100 (MQ=255)
                                                               ggaaaGAAAACAAGTACGGCCTGAGGATCACTTTCCGCGTTGCTAAAGCGCGTGATGAAGCCATTGCGTCTCTGAGCAAGCGCCATCCGGACCTGGTg                                    <  2:1261307/96‑1 (MQ=255)
                                                               gTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAt                                  <  2:281155/100‑1 (MQ=255)
                                                                tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGATGAAGCGCCATCCGGCCCTCGTgct                                  >  2:1504205/1‑97 (MQ=255)
                                                                tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                                 <  1:429839/100‑1 (MQ=255)
                                                                tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                                 <  2:1487525/100‑1 (MQ=255)
                                                                 aaaGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTa                                >  2:640503/1‑100 (MQ=255)
                                                                  aaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTag                               <  2:1556763/100‑1 (MQ=255)
                                                                    gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                             >  1:1429626/1‑100 (MQ=255)
                                                                    gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                             >  2:1425506/1‑100 (MQ=255)
                                                                     aaaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcag                            <  1:162586/100‑1 (MQ=255)
                                                                        acaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCa                         <  1:1202539/100‑1 (MQ=255)
                                                                          aacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAgg                       <  1:583518/100‑1 (MQ=255)
                                                                          aacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAgg                       <  1:153943/100‑1 (MQ=255)
                                                                          aacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAgg                       <  1:612734/100‑1 (MQ=255)
                                                                            cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                     <  1:1170447/100‑1 (MQ=255)
                                                                            cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                     <  2:1483763/100‑1 (MQ=255)
                                                                              aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                   <  2:985431/100‑1 (MQ=255)
                                                                              aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                   <  2:834156/100‑1 (MQ=255)
                                                                              aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                   <  1:219442/100‑1 (MQ=255)
                                                                              aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                   >  1:1286220/1‑100 (MQ=255)
                                                                                ggCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagca                 <  1:956154/100‑1 (MQ=255)
                                                                                 gtcTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGGCATTGCGTATCTGAGCAAGCGCCATCCGGACTTGGTGATTAGCAGCCAGGgcag                   <  2:54498/95‑1 (MQ=255)
                                                                                 gCCTGAGCATCGCTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCTATCCGGACCTGGTGATTAGCAGCCAGGGCCGCaa                >  2:829065/1‑100 (MQ=255)
                                                                                   cTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAAcc              <  2:1555269/100‑1 (MQ=255)
                                                                                        cATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtg         <  1:1492674/100‑1 (MQ=255)
                                                                                         aTCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgc        >  1:925576/1‑100 (MQ=255)
                                                                                           cACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgt      >  2:181028/1‑100 (MQ=255)
                                                                                               ttCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCgg  <  1:1132428/100‑1 (MQ=255)
                                                                                               ttCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCgg  >  1:670321/1‑100 (MQ=255)
                                                                                                tCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtgcg   >  2:1534394/1‑98 (MQ=255)
                                                                                                   |                                                                                               
AACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGG  >  NC_000913/428089‑428283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑