Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,025,587 G→A 7.1% T153T (ACG→ACA PP_5276 → phospholipase D family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,025,5870GA 7.1% 161.8 / 2.5 56T153T (ACG→ACAPP_5276phospholipase D family protein
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base A (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

CTGCCGACCGTGGCGTGCGCGTGCGCATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGCCCAACCTCAACTTCA  >  NC_002947/6025439‑6025732
                                                                                                                                                    |                                                                                                                                                 
cTGCCGACCGTGGCGTGCGCGTGCGCATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTAcg                                                                                                                                                   >  2:216561/1‑149 (MQ=255)
   ccGACCGTGGCGTGCGCGTGCGCATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTAcgc                                                                                                                                                  <  4:297798/147‑1 (MQ=255)
   ccGACCGTGGCGTGCGCGTGCGCATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTAcgc                                                                                                                                                  <  5:106915/147‑1 (MQ=255)
      aCCGTGGCGTGCGCGTGCGCATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGcc                                                                                                                                             <  8:243851/149‑1 (MQ=255)
      aCCGTGGCGTGCGCGTGCGCATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGcc                                                                                                                                             >  1:331020/1‑149 (MQ=255)
               tgcgcgtgcgcATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCgg                                                                                                                                    <  2:331020/149‑1 (MQ=255)
                gcgcgtgcgcATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGc                                                                                                                                   <  7:50539/149‑1 (MQ=255)
                  gcgtgcgcATCCTGCTCGACGACACCACCAGCGACGTCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTg                                                                                                                                 <  4:209402/149‑1 (MQ=255)
                       cgcATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCaa                                                                                                                            <  4:120235/149‑1 (MQ=255)
                             cTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTc                                                                                                                      >  5:210765/1‑149 (MQ=255)
                             cTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTc                                                                                                                      <  1:176341/149‑1 (MQ=255)
                              tGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTcg                                                                                                                     >  1:88852/1‑149 (MQ=255)
                                      gacACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCa                                                                                                             >  8:38172/1‑149 (MQ=255)
                                       acacCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCt                                                                                                                                  >  6:74340/1‑127 (MQ=255)
                                       acacCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCt                                                                                                                                  <  5:74340/127‑1 (MQ=255)
                                         accaccAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCg                                                                                                          >  3:237015/1‑149 (MQ=255)
                                          ccaccaGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCgg                                                                                                         <  1:352276/149‑1 (MQ=255)
                                          ccaccaGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCgg                                                                                                         >  8:149116/1‑149 (MQ=255)
                                           caccaGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGgc                                                                                                        >  2:165428/1‑149 (MQ=255)
                                              caGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCa                                                                                                     >  7:178334/1‑149 (MQ=255)
                                                 cGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGc                                                                                                  <  6:210765/149‑1 (MQ=255)
                                                  gACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGTACCGGCGCATGca                                                                                                 <  5:333103/149‑1 (MQ=255)
                                                           gACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAg                                                                                          <  1:216561/147‑1 (MQ=255)
                                                           gACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGc                                                                                         <  1:165428/148‑1 (MQ=255)
                                                                cGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTcc                                                                                   >  5:28147/1‑149 (MQ=255)
                                                                    aTGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCtggt                                                                               <  7:81008/149‑1 (MQ=255)
                                                                            ccTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGacaa                                                                         >  4:34867/1‑147 (MQ=255)
                                                                            ccTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGacaa                                                                         >  7:91400/1‑147 (MQ=255)
                                                                              tCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAg                                                                      >  4:130642/1‑148 (MQ=255)
                                                                                  tGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATgg                                                                 <  6:28147/149‑1 (MQ=255)
                                                                                      cATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCAt                                                             >  5:244092/1‑149 (MQ=255)
                                                                                      cATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCAt                                                             >  7:32983/1‑149 (MQ=255)
                                                                                          cGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCAtcgtc                                                         >  5:378836/1‑149 (MQ=255)
                                                                                            aaCATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTcgg                                                       >  4:225354/1‑149 (MQ=255)
                                                                                               aTCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTcggcgg                                                    >  7:76509/1‑149 (MQ=255)
                                                                                                   acaTTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCtggtg                                                                              <  5:48674/119‑1 (MQ=255)
                                                                                                   acaTTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCtggtg                                                                              >  6:48674/1‑119 (MQ=255)
                                                                                                   acaTTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCggcgggcg                                                 >  1:344499/1‑148 (MQ=255)
                                                                                                       tCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCGAcc                                            >  8:57985/1‑149 (MQ=255)
                                                                                                         gcgTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACACGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTg                                          >  5:213666/1‑149 (MQ=255)
                                                                                                         gcgTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACACGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTg                                          >  5:208678/1‑149 (MQ=255)
                                                                                                            tgtTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGc                                       <  4:237015/149‑1 (MQ=255)
                                                                                                                cAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGcgac                                    <  3:130642/148‑1 (MQ=255)
                                                                                                                 aaCCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGcgac                                    >  2:115809/1‑147 (MQ=255)
                                                                                                                    ccACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCt                                                                                  >  7:125380/1‑98 (MQ=255)
                                                                                                                    ccACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCt                                                                                  <  8:125380/98‑1 (MQ=255)
                                                                                                                    ccACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTa                               >  6:37061/1‑149 (MQ=255)
                                                                                                                       cTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTAttt                            <  2:344499/149‑1 (MQ=255)
                                                                                                                        tGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTc                           >  3:33928/1‑149 (MQ=255)
                                                                                                                           aTCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGAt                        <  8:32983/149‑1 (MQ=255)
                                                                                                                           aTCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGAt                        <  6:244092/149‑1 (MQ=255)
                                                                                                                                  gcgcAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGc                 <  7:149116/149‑1 (MQ=255)
                                                                                                                                    gcAGTACCGGCGTTACACGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGANccc               <  6:213666/149‑1 (MQ=255)
                                                                                                                                    gcAGTACCGGCGTTACACGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGccc               <  6:208678/149‑1 (MQ=255)
                                                                                                                                          ccGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATTGCCATCGTCTGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGCCCAACCTc         <  6:378836/149‑1 (MQ=255)
                                                                                                                                          ccGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGCCCAACCt          >  1:244262/1‑148 (MQ=255)
                                                                                                                                          ccGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGCCCAACCt          >  8:89914/1‑148 (MQ=255)
                                                                                                                                            ggCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGCCCAACCTCaa       <  7:57985/149‑1 (MQ=255)
                                                                                                                                                  aCGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGCCCAACCTCAACTTCa  <  8:91400/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTGCCGACCGTGGCGTGCGCGTGCGCATCCTGCTCGACGACACCACCAGCGACGGCCTGGACACCGTCATGGGCACCCTCGATGCCCATCCGAACATCCACATTCGCGTGTTCAACCCACTGCATCTGGGGCGCAGTACCGGCGTTACGCGCGCCGTGGGCCGGCTGTTCAACTTGTCGCGCCAGCACCGGCGCATGCACAACAAGCTGTTCCTGGTGGACAACAGCATGGCCATCGTCGGCGGGCGCAACCTGGGCGACGAGTATTTCGATGCCGAGCCCAACCTCAACTTCA  >  NC_002947/6025439‑6025732

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