Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,759,765 G→T 7.0% P393Q (CCG→CAG)  PP_2413 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,759,7650GT 7.0% 182.5 / 3.9 57P393Q (CCG→CAG) PP_2413GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base T (2/2);  total (30/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 = 3.11e-01

GTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCG  >  NC_002947/2759623‑2759899
                                                                                                                                                |                                                                                                                                       
gTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCgg                                                                                                                                   >  4:174871/1‑149 (MQ=255)
 tCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGt                                                                                                                                  <  1:247655/149‑1 (MQ=255)
     ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              <  6:100615/149‑1 (MQ=255)
     ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              >  2:207172/1‑149 (MQ=255)
     ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              <  2:134840/149‑1 (MQ=255)
     ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              <  1:207172/149‑1 (MQ=255)
     ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGACCGATCCGCTGACCGGCCTgg                                                                                                                              <  7:342949/149‑1 (MQ=255)
      cACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTggg                                                                                                                             >  2:329303/1‑149 (MQ=255)
              gtCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGa                                                                                                                     >  6:136231/1‑149 (MQ=255)
               tCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAt                                                                                                                    >  6:14338/1‑149 (MQ=255)
                  aTGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGa                                                                                                                 <  7:118703/149‑1 (MQ=255)
                      cTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccg                                                                                                             >  4:338985/1‑149 (MQ=255)
                        aTCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggc                                                                                                           >  2:286090/1‑149 (MQ=255)
                         tCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                          <  8:336312/149‑1 (MQ=255)
                                    cTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                               >  4:334421/1‑149 (MQ=255)
                                          gTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTgcagca                                                                                         <  3:338985/149‑1 (MQ=255)
                                             gTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAg                                                                                       <  5:124505/148‑1 (MQ=255)
                                               cGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGtt                                                                                     >  3:391063/1‑148 (MQ=255)
                                                     ccAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggccg                                                                                                         >  4:276248/1‑122 (MQ=255)
                                                     ccAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggccg                                                                                                         <  3:276248/122‑1 (MQ=255)
                                                      cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAtt                                                                                                                   >  5:275212/1‑111 (MQ=255)
                                                      cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAtt                                                                                                                   <  6:275212/111‑1 (MQ=255)
                                                      cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                          <  3:348731/120‑1 (MQ=255)
                                                      cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                          >  4:348731/1‑120 (MQ=255)
                                                      cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTg                                                                             >  5:15726/1‑148 (MQ=255)
                                                            aaTTGCTGCAAACTCCCGCTTTTCACCCTCCAGTAACGACATTACCATATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAg                                                                        <  1:136555/149‑1 (MQ=255)
                                                            aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                               >  5:192481/1‑125 (MQ=255)
                                                            aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                               <  6:192481/125‑1 (MQ=255)
                                                            aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGc                                                                       >  3:382116/1‑148 (MQ=255)
                                                                  tGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGc                                                                       <  8:352858/142‑1 (MQ=255)
                                                                  tGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAAcgcgc                                                                >  5:18978/1‑149 (MQ=255)
                                                                  tGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAAcgcgc                                                                >  6:245763/1‑149 (MQ=255)
                                                                  tGCAAACTCCCGCTATTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGc                                                                       >  7:352858/1‑142 (MQ=255)
                                                                         tCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAg                                                         >  1:271252/1‑149 (MQ=255)
                                                                         tCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAg                                                         >  2:372269/1‑149 (MQ=255)
                                                                          cccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCtt                                                            <  6:40915/145‑1 (MQ=255)
                                                                          cccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCtt                                                            >  5:40915/1‑145 (MQ=255)
                                                                          cccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGc                                                        >  6:384422/1‑149 (MQ=255)
                                                                           ccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCt                                                       >  6:280431/1‑149 (MQ=255)
                                                                           aggattttCATCCTCAAGTTACGAGA‑‑GGCTTATTTCCAGGCTTCGGCGTTTGCCCAGGCCTGTCAGCGGATCGTTCTGGGCGTCCAGATTGAGGCGACCGATACGCTGTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTTAg                                                         <  8:130876/143‑1 (MQ=255)
                                                                            cGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTTCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAAAAGCAACGCGCGCTTGAGCTc                                                      >  5:387912/1‑149 (MQ=255)
                                                                            cGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTc                                                      >  4:282246/1‑149 (MQ=255)
                                                                             gCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCg                                                     <  5:14338/149‑1 (MQ=255)
                                                                                 ttCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGgc                                                   <  3:174871/147‑1 (MQ=255)
                                                                                     gCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTc                                             <  6:18978/149‑1 (MQ=255)
                                                                                          cAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGaa                                          <  5:280431/147‑1 (MQ=255)
                                                                                          cAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCGCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAg                                         <  2:69013/148‑1 (MQ=255)
                                                                                             tAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGCGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTAc                                      <  1:286090/149‑1 (MQ=255)
                                                                                                 gACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGctc                                                                                                 >  4:64793/1‑86 (MQ=255)
                                                                                                 gACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGctc                                                                                                 <  3:64793/86‑1 (MQ=255)
                                                                                                 gACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGc                                 >  1:335701/1‑149 (MQ=255)
                                                                                                                  aGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCaa                  <  6:15726/149‑1 (MQ=255)
                                                                                                                           cggTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTcggc         <  2:335701/149‑1 (MQ=255)
                                                                                                                                gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAg                                   >  7:325621/1‑118 (MQ=38)
                                                                                                                                gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAg                                   <  8:325621/118‑1 (MQ=38)
                                                                                                                                gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTAc    >  4:283715/1‑149 (MQ=255)
                                                                                                                                  ccAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCg  >  1:123060/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                       
GTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCG  >  NC_002947/2759623‑2759899

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: