Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,996,011 C→G 14.2% intergenic (‑257/+28) PP_mr26 ← / ← rfbC C4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,996,0110CG14.2% 167.4 / 14.3 62intergenic (‑257/+28)PP_mr26/rfbCC4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base G (5/4);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCAATGCAG  >  NC_002947/1995863‑1996156
                                                                                                                                                    |                                                                                                                                                 
cgcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgcgc                                                                                                                                                   >  6:200493/1‑149 (MQ=255)
     gggTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCg                                                                                                                                              >  2:306378/1‑149 (MQ=255)
              tGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGcctt                                                                                                                                       <  4:281082/147‑3 (MQ=255)
                 gAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTcc                                                                                                                                  >  2:19799/1‑149 (MQ=255)
                   aaTGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGc                                                                                                                                <  1:227017/149‑1 (MQ=255)
                       gCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTaa                                                                                                                             >  6:177132/1‑148 (MQ=255)
                         cTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGccttgtct                                                                                                                                   >  6:293688/1‑134 (MQ=255)
                         cTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGccttgtct                                                                                                                                   <  5:293688/140‑7 (MQ=255)
                         cTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGccttgt                                                                                                                                     >  1:207816/1‑134 (MQ=255)
                         cTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGccttgt                                                                                                                                     <  2:207816/138‑5 (MQ=255)
                           aCAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGccc                                                                                                                                        >  6:329959/1‑133 (MQ=255)
                           aCAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGccc                                                                                                                                        <  5:329959/133‑1 (MQ=255)
                            cAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGt                                                                                                                       <  2:196572/149‑1 (MQ=255)
                                            tCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTcc                                                                                                                                  <  1:120018/122‑1 (MQ=255)
                                            tCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTcc                                                                                                                                  >  2:120018/1‑122 (MQ=255)
                                            tCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCgg                                                                                                       >  1:253921/1‑149 (MQ=255)
                                              tGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCg                                                                                                     >  2:187605/1‑149 (MQ=255)
                                              tGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCg                                                                                                     >  6:53898/1‑149 (MQ=255)
                                                   attTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGc                                                                                                >  8:147816/1‑149 (MQ=255)
                                                   attTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGc                                                                                                >  4:145159/1‑149 (MQ=255)
                                                     ttGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAg                                                                                                                                           <  6:332190/104‑1 (MQ=255)
                                                     ttGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAg                                                                                                                                           >  5:332190/1‑104 (MQ=255)
                                                        ccAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATg                                                                                           <  5:333311/149‑1 (MQ=255)
                                                         cAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGc                                                                                          <  8:126182/149‑1 (MQ=255)
                                                           ggCCATGCTGGGGTTGCTTCTTAGTGGCGTGTGAAAAACGGCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCg                                                                                                        <  6:229136/133‑1 (MQ=255)
                                                           ggCCATGCAGCGGTTGCTTCTTAGTCGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTAAGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCACCGCTAAAACCGTAAGTCGAACAGTTCg                                                                                                        >  5:229136/1‑133 (MQ=255)
                                                            gCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGa                                                                                        <  3:26215/148‑1 (MQ=255)
                                                            gCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGa                                                                                        >  4:26215/1‑148 (MQ=255)
                                                                 gCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCaa                                                                                  >  2:127047/1‑149 (MQ=255)
                                                                   aGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAg                                                                                 <  7:321477/148‑1 (MQ=255)
                                                                   aGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAg                                                                                 >  8:321477/1‑148 (MQ=255)
                                                                   aGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGc                                                                                <  2:253921/149‑1 (MQ=255)
                                                                     cGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTg                                                                              <  5:237299/149‑1 (MQ=255)
                                                                        ttGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTc                                                                           <  5:53898/149‑1 (MQ=255)
                                                                               ttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                                                                      >  5:154155/1‑147 (MQ=255)
                                                                                 aGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCa                                                                    <  1:187605/147‑1 (MQ=255)
                                                                                     gCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGAc                                                                 >  3:241536/1‑146 (MQ=255)
                                                                                        tgtgAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCtt                                                           <  4:178304/149‑1 (MQ=255)
                                                                                          tgAAAAACGCCACTATCGTCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTg                                                          <  5:177132/148‑1 (MQ=255)
                                                                                          tgAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTgg                                                         <  6:324080/149‑1 (MQ=255)
                                                                                          tgAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTg                                                          <  4:241536/148‑1 (MQ=255)
                                                                                            aaaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTa                                                       >  4:169799/1‑149 (MQ=255)
                                                                                             aaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAc                                                      <  1:19799/149‑1 (MQ=255)
                                                                                                aCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGa                                                   <  5:200493/149‑1 (MQ=255)
                                                                                                 cGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAg                                                  <  3:145159/149‑1 (MQ=255)
                                                                                                    cACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTc                                                >  8:306291/1‑148 (MQ=255)
                                                                                                         tCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATg                                            <  3:169799/147‑1 (MQ=255)
                                                                                                           ggCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTAACGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGc                                         >  8:212342/1‑148 (MQ=255)
                                                                                                             cAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCAt                                      >  7:175248/1‑149 (MQ=255)
                                                                                                                gACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTc                                    >  4:139642/1‑148 (MQ=255)
                                                                                                                 aCCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGa                                                   >  7:290547/1‑132 (MQ=255)
                                                                                                                 aCCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGa                                                   <  8:290547/132‑1 (MQ=255)
                                                                                                                             aTGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCg                      >  5:258235/1‑149 (MQ=255)
                                                                                                                              tgtgTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCgg                     <  8:175248/149‑1 (MQ=255)
                                                                                                                                    cAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCAt                <  1:127047/148‑1 (MQ=255)
                                                                                                                                           aaGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTg                                                          <  6:6319/99‑1 (MQ=255)
                                                                                                                                           aaGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTg                                                          >  5:6319/1‑99 (MQ=255)
                                                                                                                                            gggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCtt                                                                        <  2:173109/83‑1 (MQ=255)
                                                                                                                                            gggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCtt                                                                        >  1:173109/2‑84 (MQ=255)
                                                                                                                                            gggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCaa       >  2:211270/2‑149 (MQ=255)
                                                                                                                                            gggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCaa       >  2:87550/2‑149 (MQ=255)
                                                                                                                                             ggCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATgg                                           <  3:191297/112‑1 (MQ=255)
                                                                                                                                             ggCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATgg                                           >  4:191297/1‑112 (MQ=255)
                                                                                                                                                 gcgcTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCAATGCAg  <  7:306291/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCAATGCAG  >  NC_002947/1995863‑1996156

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: ‑