Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,713,568 T→C W225R (TGG→CGG)  srmB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,713,5680TC100.0% 263.4 / NA 73W225R (TGG→CGG) srmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/35);  total (38/35)

ACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAA  >  NC_000913/2713445‑2713691
                                                                                                                           |                                                                                                                           
aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTAcc                                                                                                                    >  2:1522995/1‑133 (MQ=255)
aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTAcc                                                                                                                    <  1:1522995/133‑1 (MQ=255)
aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                              >  1:957332/1‑139 (MQ=255)
aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                              >  2:560475/1‑139 (MQ=255)
aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                              >  2:200576/1‑139 (MQ=255)
aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGAATCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                              >  2:1360787/1‑139 (MQ=255)
     gCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatgat                                                                                                         >  1:930775/1‑139 (MQ=255)
      ctctTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                        <  1:560475/139‑1 (MQ=255)
                cAACGCTGGAAGGCTATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAATTTGCTGCCAATCCCTCCACCCGTGTGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATaa                                                                                              <  2:1023498/139‑1 (MQ=255)
                     cTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                         <  2:957332/139‑1 (MQ=255)
                     cTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                         <  1:1009047/139‑1 (MQ=255)
                       ggAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgcg                                                                                       <  1:200576/139‑1 (MQ=255)
                           ggCGATGCGATTCAGGACTTTGCCCAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGc                                                                                   <  2:930775/139‑1 (MQ=255)
                                  cgaaTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                            <  1:1360787/139‑1 (MQ=255)
                                          gACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTaa                                                                    >  2:228109/1‑139 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgat                                                                                                            <  1:1420791/94‑1 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgat                                                                                                            >  2:1420791/1‑94 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCtt                                                                                                      >  2:934341/1‑100 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCtt                                                                                                      <  1:934341/100‑1 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgc                                                                                        >  2:1038968/1‑114 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgc                                                                                        <  1:1038968/114‑1 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                               >  2:1605186/1‑139 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                               >  2:216562/1‑139 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                               >  2:1239136/1‑139 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                               >  1:1418639/1‑139 (MQ=255)
                                               tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                               >  2:703050/1‑139 (MQ=255)
                                                       gTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGaa                                                         >  2:279867/1‑137 (MQ=255)
                                                       gTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGaa                                                         <  1:279867/137‑1 (MQ=255)
                                                             tgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCg                                                 <  2:608715/139‑1 (MQ=255)
                                                              ggAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                <  1:703050/139‑1 (MQ=255)
                                                                    tCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAAtt                                          <  2:1418639/139‑1 (MQ=255)
                                                                      cGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                           <  1:553451/72‑1 (MQ=255)
                                                                      cGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                           >  2:553451/1‑72 (MQ=255)
                                                                        gTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTgtgt                                      <  2:1239495/139‑1 (MQ=255)
                                                                         tGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGttt                                    >  2:863399/1‑139 (MQ=255)
                                                                              gTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGc                                <  1:1646971/139‑1 (MQ=255)
                                                                                   tGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                            <  1:785859/90‑1 (MQ=255)
                                                                                   tGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                            >  2:785859/1‑90 (MQ=255)
                                                                                     ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCtt                                                                                                      <  1:2320/62‑1 (MQ=255)
                                                                                     ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCtt                                                                                                      >  2:2320/1‑62 (MQ=255)
                                                                                     ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcg                         >  2:1381515/1‑139 (MQ=255)
                                                                                       aaTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtg                       <  1:1605186/139‑1 (MQ=255)
                                                                                        aTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtga                      <  1:216562/139‑1 (MQ=255)
                                                                                        aTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtga                      <  1:1239136/139‑1 (MQ=255)
                                                                                         tCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtgag                     >  2:1182415/1‑139 (MQ=255)
                                                                                            cTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTaaaa                                                                  >  1:1296306/1‑91 (MQ=255)
                                                                                            cTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTaaaa                                                                  <  2:1296306/91‑1 (MQ=255)
                                                                                             tCCACCCTTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                        <  1:1182798/52‑1 (MQ=39)
                                                                                             tCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                        >  2:1182798/1‑52 (MQ=255)
                                                                                               cACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACcgc                                                                                                                  <  2:487356/40‑1 (MQ=255)
                                                                                               cACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACcgc                                                                                                                  >  1:487356/1‑40 (MQ=255)
                                                                                               cACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                           <  2:898815/47‑1 (MQ=255)
                                                                                               cACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                           >  1:648254/1‑47 (MQ=255)
                                                                                               cACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                           <  2:648254/47‑1 (MQ=255)
                                                                                               cACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                           >  1:898815/1‑47 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCt                                               >  2:838349/1‑102 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCt                                               <  1:838349/102‑1 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTtgt                                  >  2:560884/1‑115 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTtgt                                  <  1:560884/115‑1 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTtgt                                  >  2:1640509/1‑115 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTtgt                                  <  1:1640509/115‑1 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtgagc                    >  1:188204/1‑129 (MQ=255)
                                                                                                    gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtgagc                    <  2:188204/129‑1 (MQ=255)
                                                                                                     tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCAt                                                                           >  1:333684/1‑73 (MQ=255)
                                                                                                     tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCAt                                                                           <  2:333684/73‑1 (MQ=255)
                                                                                                     tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAg         >  2:417120/1‑139 (MQ=255)
                                                                                                       aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGc                                                                                   >  2:12298/1‑63 (MQ=255)
                                                                                                       aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGc                                                                                   <  1:12298/63‑1 (MQ=255)
                                                                                                       aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCt       >  1:589079/1‑139 (MQ=255)
                                                                                                         cgcAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTgg     >  2:910036/1‑139 (MQ=255)
                                                                                                            aaaaaaaTTCATCAGCGGTATTACCGCGCCGATGATCTTGAg                                                                                                   >  2:765430/1‑42 (MQ=255)
                                                                                                            aaaaaaaTTCATCAGCGGTATTACCGCGCCGATGATCTTGAg                                                                                                   <  1:765430/42‑1 (MQ=255)
                                                                                                            aaaaaaaTTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCaa  >  1:963503/1‑139 (MQ=255)
                                                                                                                           |                                                                                                                           
ACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAA  >  NC_000913/2713445‑2713691

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