Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,082:1 +G 9.5% intergenic (+36/‑277) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,0821.G 9.5% 224.9 / 20.5 63intergenic (+36/‑277)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base . (26/31);  new base G (3/3);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

GGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCG  >  NC_002947/954937‑955177
                                                                                                                                                     |                                                                                               
gggCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGcc                                                                                              >  7:60934/1‑149 (MQ=255)
              tttggCTTCTCCATCCTCACCGTCACTGCCGTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGttttttttt                                                                                 <  6:201801/145‑1 (MQ=255)
              ttCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGttttttttg                                                                                 <  2:246275/148‑2 (MQ=255)
               tCGGCCTCTCCATCCTCCCCGGCACTGCCCTGAGCGTCATGATGGGTGCGGACAGCCAGACCGCTCCGGATGCCCTGATGACCTCCGC‑‑ACGCTGAAACTGCTGCCTGCTTCACGAAAACCCCCGCCACTGTGCCGGGGTTTTTTTTTgg                                                                               <  6:60761/149‑1 (MQ=255)
                 ggCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTgggg                                                                             <  7:169900/149‑1 (MQ=255)
                 ggCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTggg                                                                              <  5:35331/148‑1 (MQ=255)
                    ttCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTa                                                                          >  7:252514/1‑149 (MQ=255)
                         cATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGttttc                                                                                     >  5:187486/1‑132 (MQ=255)
                         cATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGttttc                                                                                     <  6:187486/133‑2 (MQ=255)
                         cATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAttttt                                                                     >  5:51769/1‑149 (MQ=255)
                           tCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTcc                                                                   <  6:51769/149‑1 (MQ=255)
                             atCACCTTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCt                                                                  <  4:266312/147‑1 (MQ=255)
                              tCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGGttttcgt                                                                                   >  5:340828/1‑127 (MQ=255)
                              tCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGGttttcgt                                                                                   <  6:340828/130‑4 (MQ=255)
                               cACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTgg                                                                               >  3:137344/1‑133 (MQ=255)
                               cACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTgg                                                                               <  4:137344/133‑1 (MQ=255)
                                  cGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCTGAGCGCTGAAACTGGTGCCTGTTTCACAAAAACCCCGGCCAATGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCtttg                                                               >  1:108804/1‑149 (MQ=255)
                                  cGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCg                                                            >  2:134356/1‑149 (MQ=255)
                                          cccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTc                                                    >  8:306888/1‑149 (MQ=255)
                                            cTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTggg                                                                              >  3:198054/1‑121 (MQ=255)
                                            cTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTggg                                                                              <  4:198054/121‑1 (MQ=255)
                                                gCGTGAGGACGGGTGTGGACAGCCAGTCCGCGGCGGAGGTGCGGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGt                                              <  8:141852/149‑1 (MQ=255)
                                                     atgatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc                                                  >  4:241500/1‑140 (MQ=255)
                                                     atgatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc                                                  <  3:241500/140‑1 (MQ=255)
                                                         tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGc                                                              >  2:329364/1‑124 (MQ=255)
                                                         tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGc                                                              <  1:329364/124‑1 (MQ=255)
                                                         tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCt                                     >  7:42101/1‑149 (MQ=255)
                                                             tgtgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                 <  7:306888/149‑1 (MQ=255)
                                                                ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAtttt                                                                      <  6:312630/109‑1 (MQ=255)
                                                                ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAtttt                                                                      >  5:312630/1‑109 (MQ=255)
                                                                 gACAGCCAGACCGCTGCGGATGTGCTTATGACCTACGC‑‑ACGCTGAAACTTCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAg                             <  7:193720/149‑1 (MQ=255)
                                                                 gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                              <  8:42101/148‑1 (MQ=255)
                                                                  aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACagg                             >  5:279411/1‑146 (MQ=255)
                                                                    aGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTgg                                                                               >  6:353538/1‑96 (MQ=255)
                                                                     gCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaggaacc                         <  3:270581/149‑7 (MQ=255)
                                                                      ccAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTATTAGCGGGTAAAACCCCTCCTACAATGttac                        >  3:103900/1‑147 (MQ=255)
                                                                        aGACCGCTGCGGATGTGCTGTTGACCTACGC‑‑CCGCTGAAACTGCTGCCTGTTTCACGAAAACCGCGGCCACTGTGCCGGGGTTTTTTTTTtgggg                                                                            <  6:16025/95‑1 (MQ=255)
                                                                        aGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                      <  2:108804/149‑1 (MQ=255)
                                                                          aCCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGGt                                                                                         <  5:209257/80‑1 (MQ=38)
                                                                          aCCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGGt                                                                                         <  5:204272/80‑1 (MQ=38)
                                                                          aCCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGGt                                                                                         >  6:204272/1‑80 (MQ=38)
                                                                          aCCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGGt                                                                                         >  6:209257/1‑80 (MQ=38)
                                                                           ccGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                   <  1:134356/149‑1 (MQ=255)
                                                                             gCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                 <  8:252514/149‑1 (MQ=255)
                                                                                  ggATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAAGGTTCCTGTAGGAGCggg            >  1:235775/1‑149 (MQ=255)
                                                                                  gcctgtgCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTATCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg            <  8:60934/146‑1 (MQ=255)
                                                                                          tgatgaCCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCGGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGtt                                             >  4:44676/1‑107 (MQ=255)
                                                                                          tgatgaCCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACa                               <  6:251088/122‑1 (MQ=255)
                                                                                          tgatgaCCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGt                           >  7:167882/1‑126 (MQ=255)
                                                                                          tgatgaCCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGt                           <  8:167882/126‑1 (MQ=255)
                                                                                               aCCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg            <  2:45276/136‑1 (MQ=37)
                                                                                               aCCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg            >  1:45276/1‑136 (MQ=37)
                                                                                                 cTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg            <  7:335394/134‑1 (MQ=37)
                                                                                                 cTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg            >  8:335394/1‑134 (MQ=37)
                                                                                                 cTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg             <  4:127436/133‑1 (MQ=37)
                                                                                                 cTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg             >  3:127436/1‑133 (MQ=37)
                                                                                                                         ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTaaa                                           >  6:265304/1‑82 (MQ=255)
                                                                                                                                 aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg            >  3:195223/1‑105 (MQ=21)
                                                                                                                                 aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg            <  4:195223/105‑1 (MQ=21)
                                                                                                                                          ccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt           <  4:127123/97‑1 (MQ=255)
                                                                                                                                          ccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt           >  3:127123/1‑97 (MQ=255)
                                                                                                                                          ccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgcg  <  3:237901/106‑1 (MQ=255)
                                                                                                                                          ccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgcg  >  4:237901/1‑106 (MQ=255)
                                                                                                                                                     |                                                                                               
GGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGC‑‑ACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCG  >  NC_002947/954937‑955177

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