Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,199,332 C→T S181S (AGC→AGT nfi → endonuclease V, deoxyinosine 3' endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,199,3320CT100.0% 213.3 / NA 63S181S (AGC→AGTnfiendonuclease V, deoxyinosine 3' endonuclease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/30);  total (33/30)

AAATTCGAACCGCTCTCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGCGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCAGCCCTAATTC  >  NC_000913/4199201‑4199464
                                                                                                                                   |                                                                                                                                    
aaaTTCGAACCGCTCTCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACa                                                                                                                               <  2:1120856/139‑1 (MQ=255)
 aaTTCGAACCGCTCTCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAg                                                                                                                              >  1:683067/1‑139 (MQ=255)
      gAACCGCTCTCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAgcgcgc                                                                                                                         >  1:1319672/1‑139 (MQ=255)
          cGCTCTCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGc                                                                                                                     <  2:1440737/139‑1 (MQ=255)
           gCTCTCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCg                                                                                                                    >  2:856749/1‑139 (MQ=255)
            ctctCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGt                                                                                                                   <  1:1415701/139‑1 (MQ=255)
                ccAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGt                                                                                                               <  1:856749/139‑1 (MQ=255)
                      aaCCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACa                                                                                                                               >  1:348999/1‑117 (MQ=255)
                      aaCCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACa                                                                                                                               <  2:348999/117‑1 (MQ=255)
                      aaCCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACg                                                                                                         >  1:164342/1‑139 (MQ=255)
                      aaCCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACg                                                                                                         >  2:784699/1‑139 (MQ=255)
                       aCCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGc                                                                                                        >  1:401671/1‑139 (MQ=255)
                        ccGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCt                                                                                                       <  1:784699/139‑1 (MQ=255)
                        ccGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCt                                                                                                       <  1:687226/139‑1 (MQ=255)
                         cGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTg                                                                                                      >  1:336276/1‑139 (MQ=255)
                             cgcgcTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATg                                                                                                  <  1:1122448/139‑1 (MQ=255)
                               cgcTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGaa                                                                                                <  2:1128706/139‑1 (MQ=255)
                                             aTGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGc                                                                                  <  2:164342/139‑1 (MQ=255)
                                                     aGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCa                                                                          >  1:583173/1‑139 (MQ=255)
                                                       gCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGt                                                                                      >  2:567265/1‑125 (MQ=255)
                                                       gCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGt                                                                                      <  1:567265/125‑1 (MQ=255)
                                                       gCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAAc                                                                        >  1:1176197/1‑139 (MQ=255)
                                                          agcagcTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACaa                                                                                                           <  1:967909/101‑1 (MQ=255)
                                                          agcagcTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACaa                                                                                                           >  2:967909/1‑101 (MQ=255)
                                                           gcagcTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACgcgc                                                                    >  2:686258/1‑139 (MQ=255)
                                                           gcagcTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACgcgc                                                                    >  1:1056177/1‑139 (MQ=255)
                                                               cTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTggg                                                                >  1:1393002/1‑139 (MQ=255)
                                                                   ccTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGa                                                            <  2:401671/139‑1 (MQ=255)
                                                                       ggTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGAcgcg                                                        >  2:396393/1‑139 (MQ=255)
                                                                            ggCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGc                                                   >  1:300626/1‑139 (MQ=255)
                                                                            ggCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGc                                                   >  2:238521/1‑139 (MQ=255)
                                                                            ggCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGc                                                   >  1:1242476/1‑139 (MQ=255)
                                                                                  gcaAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGt                                                                                                                   <  2:1060293/69‑1 (MQ=255)
                                                                                  gcaAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGt                                                                                                                   >  1:1060293/1‑69 (MQ=255)
                                                                                   caAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGaa                                            >  2:682263/1‑139 (MQ=255)
                                                                                        cgcgcTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGAGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTCTCGTCTGCCGGAGCCAACGCGATGGGCAGACGCGGCGGCCTCGGAACGtac                                       >  2:1116677/1‑137 (MQ=255)
                                                                                         gcgcTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCg                                      <  2:1176197/139‑1 (MQ=255)
                                                                                             tGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAAcgcg                                                                     <  1:303634/104‑1 (MQ=255)
                                                                                             tGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAAcgcg                                                                     >  2:303634/1‑104 (MQ=255)
                                                                                             tGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGAcgcg                                                        >  2:1158691/1‑117 (MQ=255)
                                                                                             tGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGAcgcg                                                        <  1:1158691/117‑1 (MQ=255)
                                                                                             tGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGt                                  <  2:683067/139‑1 (MQ=255)
                                                                                              gTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGtt                                 >  1:579439/1‑139 (MQ=255)
                                                                                               tAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTc                                >  1:625859/1‑139 (MQ=255)
                                                                                               tAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTc                                <  1:953552/139‑1 (MQ=255)
                                                                                                    cGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTgcg                           <  2:1319672/139‑1 (MQ=255)
                                                                                                    cGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTgcg                           <  2:712749/139‑1 (MQ=255)
                                                                                                    cGTTGTTTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTgcg                           <  2:583173/139‑1 (MQ=255)
                                                                                                     gttgttTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTAc                                                                                                             >  2:1200243/1‑56 (MQ=255)
                                                                                                     gttgttTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTAc                                                                                                             <  1:1200243/56‑1 (MQ=255)
                                                                                                       tgttTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCta                        <  1:101089/139‑1 (MQ=255)
                                                                                                       tgttTATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCta                        >  1:363203/1‑139 (MQ=255)
                                                                                                         tttATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCtata                      >  1:1278815/1‑139 (MQ=255)
                                                                                                           tATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCgg                                     >  1:1399102/1‑122 (MQ=255)
                                                                                                           tATCGCTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCgg                                     <  2:1399102/122‑1 (MQ=255)
                                                                                                                cTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTc                                <  1:630375/122‑1 (MQ=255)
                                                                                                                cTACCGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTc                                >  2:630375/1‑122 (MQ=255)
                                                                                                                   ccGGCCATCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCa            >  1:447137/1‑139 (MQ=255)
                                                                                                                          tCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCAGCCCTaa     <  2:300626/139‑1 (MQ=255)
                                                                                                                          tCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCAGCCCTaa     >  2:452009/1‑139 (MQ=255)
                                                                                                                          tCGGGTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCAGCCCTaa     <  2:314265/139‑1 (MQ=255)
                                                                                                                            gggTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCAGCCCTAAtt   <  1:238521/139‑1 (MQ=255)
                                                                                                                             ggTCAGTGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCAGCCCTAATTc  <  2:447137/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
AAATTCGAACCGCTCTCCAGCGAACCGGGCGCGCTGGCCCCACTGATGGATAAAGGCGAGCAGCTGGCCTGGGTCTGGCGCAGCAAAGCGCGCTGTAACCCGTTGTTTATCGCTACCGGCCATCGGGTCAGCGTGGACAGCGCGCTGGCGTGGGTACAACGCTGCATGAAAGGCTATCGTCTGCCGGAGCCAACGCGCTGGGCGGACGCGGTGGCCTCGGAACGTCCGGCGTTCGTGCGCTATACAGCAAATCAGCCCTAATTC  >  NC_000913/4199201‑4199464

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: