Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,820,598 G→C 13.0% intergenic (+399/+113) PP_2473 → / ← PP_2474 hypothetical protein/glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,820,5980GC13.0% 160.0 / 11.4 58intergenic (+399/+113)PP_2473/PP_2474hypothetical protein/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base C (4/3);  total (31/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 = 9.34e-01

AGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTTGT  >  NC_002947/2820457‑2820742
                                                                                                                                             |                                                                                                                                                 
aGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCgg                                                                                                                                            <  8:120769/149‑1 (MQ=255)
aGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCg                                                                                                                                             <  4:160276/148‑1 (MQ=255)
  ggATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGt                                                                                                                                          <  8:192767/149‑1 (MQ=255)
  ggATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGt                                                                                                                                          >  4:286545/1‑149 (MQ=255)
      gCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCTGGTCGAACCCGGGTTTTc                                                                                                                                      >  2:334094/1‑149 (MQ=255)
       cccAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTTCGACCCGGGttttt                                                                                                                                      <  8:146564/148‑2 (MQ=255)
       cccAGGCGTTAGTTCGAAAGTGTCTCCTCTTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAATCAAAAACCCGGGTCGAACCCGGGTTTTCg                                                                                                                                     <  1:292306/149‑1 (MQ=255)
        ccAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGtttt                                                                                                                                       <  2:216275/146‑1 (MQ=255)
        ccAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGtttt                                                                                                                                       >  1:216275/1‑146 (MQ=255)
        ccAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCgg                                                                                                                                    >  1:192702/1‑149 (MQ=255)
            gggTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGa                                                                                                                                <  3:286545/149‑1 (MQ=255)
                   ttCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTc                                                                                                                         >  3:176442/1‑149 (MQ=255)
                     cGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGTGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCg                                                                                                                        <  4:176442/148‑1 (MQ=255)
                         aGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGa                                                                                                                   <  7:142415/149‑1 (MQ=255)
                                 ctcGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCtttt                                                                                                           >  1:272580/1‑149 (MQ=255)
                                    gTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTggg                                                                                                        <  5:223451/149‑1 (MQ=255)
                                                tatCCCATGATCGGTGGTCCGTGAGGTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTTATAAAGCATAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGa                                                                                             <  8:58602/148‑1 (MQ=255)
                                                tatCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGa                                                                                             <  2:272580/148‑1 (MQ=255)
                                                  tCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGc                                                                                           >  1:264361/1‑148 (MQ=255)
                                                     cATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCg                                                                                                                        >  7:17900/1‑116 (MQ=255)
                                                     cATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCg                                                                                                                        <  8:17900/116‑1 (MQ=255)
                                                     cATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGCTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAAtgc                                                                                      >  1:113606/1‑149 (MQ=255)
                                                      aTGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAAtgc                                                                                      >  2:30854/1‑149 (MQ=255)
                                                            ggtggtCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGctct                                                                                <  2:50971/149‑1 (MQ=255)
                                                             gtggtCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTg                                                                                                                                 <  4:2495/99‑1 (MQ=255)
                                                             gtggtCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTg                                                                                                                                 >  3:2495/1‑99 (MQ=255)
                                                             gtggtCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCtctc                                                                               <  2:192702/149‑1 (MQ=255)
                                                                          tAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTgg                                                                 >  7:99455/1‑149 (MQ=255)
                                                                              aGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTc                                                             >  2:35964/1‑149 (MQ=255)
                                                                              aGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTc                                                             >  1:177931/1‑149 (MQ=255)
                                                                                      tATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTCGTTTCTTGGAGGc                                                     >  2:187941/1‑149 (MQ=255)
                                                                                               gctgATATAGAGAAAGTGGGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCCGGCGTTGGTTTCTTGGAGGCttttttt                                              <  1:187941/148‑1 (MQ=255)
                                                                                               gctgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGc                                            <  1:334094/149‑1 (MQ=255)
                                                                                                     ataGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGc                                       >  1:312054/1‑148 (MQ=255)
                                                                                                        gagaAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGCACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAgg                                   >  3:260216/1‑149 (MQ=255)
                                                                                                        gagaAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAgg                                   >  3:264750/1‑149 (MQ=255)
                                                                                                           aaaGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa                                                        <  4:166587/125‑1 (MQ=255)
                                                                                                           aaaGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa                                                        >  3:166587/1‑125 (MQ=255)
                                                                                                               tGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAAtt                            <  1:35964/149‑1 (MQ=255)
                                                                                                               tGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAAtt                            <  2:264361/149‑1 (MQ=255)
                                                                                                                   tCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCa                        <  2:177931/149‑1 (MQ=255)
                                                                                                                   tCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCa                        >  2:144631/1‑149 (MQ=255)
                                                                                                                         aaaaaGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCtata                  >  8:234025/1‑149 (MQ=255)
                                                                                                                          tcaaccgaaaaCCCGGGTTCGACCCGGGTTTTTGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTatat                 >  3:34408/8‑149 (MQ=255)
                                                                                                                            aaccgaaaaCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATc                >  1:112047/6‑148 (MQ=255)
                                                                                                                             aGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATc                >  2:2679/1‑147 (MQ=255)
                                                                                                                             aGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCa              <  1:2679/149‑1 (MQ=255)
                                                                                                                              ccgaaaaCCCGGGTTAGACCCGGGTTTTAGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTAGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGGTTAATTGGCAGCTATATACAg             >  8:237425/4‑149 (MQ=255)
                                                                                                                                gaaaaCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGtt                                                                                                                     <  5:225216/43‑1 (MQ=17)
                                                                                                                                gaaaaCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGtt                                                                                                                     >  6:225216/2‑44 (MQ=17)
                                                                                                                                aaaaaCCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTg           >  6:15986/1‑149 (MQ=255)
                                                                                                                                 aaaaCCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTATAATTGGCATCTATATCCAGTGc          >  2:303856/1‑149 (MQ=255)
                                                                                                                                    aCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGtt                                                               <  2:13396/93‑1 (MQ=255)
                                                                                                                                    aCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGtt                                                               >  1:13396/1‑93 (MQ=255)
                                                                                                                                      ccGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTg           <  4:130874/143‑1 (MQ=255)
                                                                                                                                      ccGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTg           >  3:130874/1‑143 (MQ=255)
                                                                                                                                        gggTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTTg   <  1:144631/149‑1 (MQ=255)
                                                                                                                                          gTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTATTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTTGt  >  5:178412/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
AGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTTGT  >  NC_002947/2820457‑2820742

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