Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 497,139:1 (T)5→6 100% intergenic (+145/‑36) htpG → / → adk protein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase/adenylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913497,1341.T92.6% 207.9 / 11.0 68intergenic (+140/‑41)htpG/adkprotein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base T (35/28);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGCTCA  >  NC_000913/497035‑497227
                                                                                                    |                                                                                              
tttcATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGc                                                                                                 <  1:1948541/100‑1 (MQ=255)
  tcatcaTCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtt                                                                                              <  1:562208/100‑1 (MQ=255)
  tcatcaTCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtt                                                                                              <  1:1669569/100‑1 (MQ=255)
   catcatCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttt                                                                                             >  2:394717/1‑100 (MQ=255)
   catcatCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttt                                                                                             >  2:1473199/1‑100 (MQ=255)
     tcatcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttt                                                                                           >  1:334196/1‑100 (MQ=255)
      catcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtttttt                                                                                          <  1:325656/100‑2 (MQ=255)
        tcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttca                                                                                        >  2:136715/1‑97 (MQ=255)
        tatGCACTTTCCGGAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttca                                                                                        >  1:1727286/3‑97 (MQ=255)
         cTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttcaa                                                                                       <  2:1532542/100‑5 (MQ=255)
          tGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttcaaa                                                                                      <  1:876248/100‑6 (MQ=255)
          tGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttcaaa                                                                                      >  1:924493/1‑95 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTcaaaa                                                                                    <  1:2179471/100‑1 (MQ=255)
            tccgatctgCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCaaaaaa                                                                                   <  1:2539046/92‑1 (MQ=255)
              cTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAAt                                                                                  >  2:506075/1‑100 (MQ=255)
                ttCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTc                                                                                <  1:1467846/100‑1 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTTAAAAAATTCGac                                                                             >  2:467601/1‑100 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAAttt                                                                                >  2:1735560/1‑96 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGac                                                                             <  2:1880178/100‑1 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCACAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGac                                                                             >  1:2425901/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGcca                                                                            >  1:1065313/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCaaaaaaagatcgga                                                                           >  2:2539031/1‑92 (MQ=255)
                     cAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGacac                                                                           >  1:295743/1‑100 (MQ=255)
                      aaaTTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGAcaca                                                                          >  2:811899/1‑100 (MQ=255)
                       aaTTATGTCGCCATTAACCGTTTCCGCCCCAGGTGCCTTTCTTTAGTCAATCGCCTGTTTGTGGTATCTTTTATCGCTTTTTTCAAAAAATTCGacac                                                                           <  2:1976131/98‑1 (MQ=255)
                                gCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAgggg                                                                >  1:1093921/1‑100 (MQ=255)
                                  cATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAggggct                                                              >  2:1549023/1‑98 (MQ=255)
                                   aTTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAtt                                                             >  1:755162/1‑100 (MQ=255)
                                     tAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCATTAAATTTGAAACAGATTACGGGGAtttt                                                           >  2:1752955/1‑100 (MQ=255)
                                     tAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAtttt                                                           <  1:1623955/100‑1 (MQ=255)
                                     tAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAtttt                                                           >  2:2439299/1‑100 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTc                                                          >  2:2424896/1‑99 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                         >  1:2244778/1‑100 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                         <  1:1983214/100‑1 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                         <  1:1452315/100‑1 (MQ=255)
                                        ccGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                        <  1:2177105/100‑1 (MQ=255)
                                         cgtgttaGCCTCAGGTGTCTTTTTTGAGGCAATCGCATGTTGCTGGGATCGTGTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                        <  2:2150622/93‑1 (MQ=255)
                                         cGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCa                                                       >  1:680841/1‑100 (MQ=255)
                                         cGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCa                                                       <  1:1190891/100‑1 (MQ=255)
                                            ttCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATg                                                    <  1:1283209/100‑1 (MQ=255)
                                             tCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGc                                                   >  1:2328736/1‑100 (MQ=255)
                                              cAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                                  <  1:1715083/100‑1 (MQ=255)
                                                gCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTa                                                >  1:82618/1‑100 (MQ=255)
                                                gCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTa                                                >  1:1777984/1‑100 (MQ=255)
                                                   cctggTGCGTTTCTTTAGGCAATCTCCTGTTGGTGGTATCTTTTATCGCTTTTTGCAAAAAATTCGACACATTTTATTGGGATTTTCGCAATGCGTAt                                               <  2:1447048/95‑1 (MQ=255)
                                                          ccTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCt                                      >  2:700708/1‑100 (MQ=255)
                                                            tttCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTg                                    <  1:229655/100‑1 (MQ=255)
                                                            tttCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTg                                    <  2:583821/100‑1 (MQ=255)
                                                              tCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgc                                  <  1:2378642/100‑1 (MQ=255)
                                                               cTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgcg                                 >  2:2463685/1‑100 (MQ=255)
                                                               cTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgcg                                 <  1:639483/100‑1 (MQ=255)
                                                                  gAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTTTTTTTTTTTTACACATTTTAAGGGGATTTTTGCATTGCGTTTTATTCTGCTTGGCGCt                               >  2:1755601/1‑99 (MQ=255)
                                                                   aGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTcc                             <  2:316086/100‑1 (MQ=255)
                                                                       aaTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgc                         >  1:2245949/1‑100 (MQ=255)
                                                                       aaTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgc                         <  1:394721/100‑1 (MQ=255)
                                                                       aaTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgc                         >  2:2425284/1‑100 (MQ=255)
                                                                       aaTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATGTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgc                         >  2:444326/1‑100 (MQ=255)
                                                                        aTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAATGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgcg                        <  2:2064849/100‑1 (MQ=255)
                                                                        aTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgcg                        <  1:1830010/100‑1 (MQ=255)
                                                                          cGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGcgcg                      <  1:1303601/100‑1 (MQ=255)
                                                                           gCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgg                     >  2:1468393/1‑100 (MQ=255)
                                                                              tGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgggg                   >  1:189111/1‑99 (MQ=255)
                                                                              tGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGa                  <  1:757545/100‑1 (MQ=255)
                                                                               gTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGaa                 >  1:1838223/1‑100 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg              <  2:2403244/100‑1 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg              <  2:629487/100‑1 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg              <  1:1159347/100‑1 (MQ=255)
                                                                                   gtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAggg             <  1:1941480/100‑1 (MQ=255)
                                                                                   gtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAggg             >  1:155203/1‑100 (MQ=255)
                                                                                   gtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAggg             >  1:1453228/1‑100 (MQ=255)
                                                                                   gtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATGTTAAGGGGATTTTCGCAATGCGTATCATTCTGGTTGGCGCTCCGGGCGCGGGGAAAggg             <  1:1179743/100‑1 (MQ=255)
                                                                                   gtggtATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGa            >  2:2127968/1‑100 (MQ=255)
                                                                                     ggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGAc           >  2:2473382/1‑100 (MQ=255)
                                                                                        aTCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCa        <  2:299953/100‑1 (MQ=255)
                                                                                          cGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAgg      >  2:2077375/1‑100 (MQ=255)
                                                                                           gTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGTGCTCCGGGCGCGGGGAAAGGGACTCAGGc     >  2:1666674/1‑100 (MQ=255)
                                                                                           gTTTATCGCTTTTTTCAAAAAATTCGACACATATTAAGGGGATTTTTGCAATGCCTATCATTCTGCTCGGCGCTCCAGACGCGCGGGATGGGACTCa        >  2:2065904/1‑97 (MQ=255)
                                                                                              tATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGCTCa  >  1:428195/1‑100 (MQ=255)
                                                                                              tATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGCTCa  <  1:1216504/100‑1 (MQ=255)
                                                                                               aTCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGAGAGGGGAAAGGGACTCAGGCTCa  >  2:2385947/1‑99 (MQ=255)
                                                                                                    |                                                                                              
TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGCTCA  >  NC_000913/497035‑497227

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