Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,968,712 G→C 6.8% A50A (GCC→GCG gph ← phosphoglycolate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,968,7120GC 6.8% 370.5 / 14.2 118A50A (GCC→GCGgphphosphoglycolate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (57/53);  new base C (4/4);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 122 total aligned reads displayed.

CGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCAT  >  NC_002947/1968564‑1968861
                                                                                                                                                    |                                                                                                                                                     
cGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCggg                                                                                                                                                      <  1:1777819/150‑1 (MQ=255)
  aTGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGgcg                                                                                                                                                    >  2:670839/1‑150 (MQ=255)
      cGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCAc                                                                                                                                                <  1:1580630/150‑1 (MQ=255)
          cagcagCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatggtt                                                                                                                                            <  1:1426336/150‑6 (MQ=255)
            cacctgcGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCAGGGTCCATGGCGAAAGAAGTGGAGACCATGGGGCGGGCGCCACCCgaga                                                                                                                                          >  1:1289797/7‑150 (MQ=255)
             cagcTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGAt                                                                                                                                         <  1:1157505/150‑1 (MQ=255)
               gcTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCa                                                                                                                                       >  2:1766111/1‑150 (MQ=255)
                ctctGCCGTGCCCTCGAACAGCGTGCTGGTGACCGCGCAGTCGCGCTGGTAACGCTCGAGGAAGTCCAGGCGCAGTGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACGATGGCGCGGGCGCCACCCGAGATCAc                                                                                                                                      <  2:144768/150‑1 (MQ=255)
                 tCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCGGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACg                                                                                                                                     >  2:877598/1‑150 (MQ=255)
                     ccATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCAc                                                                                                                                 >  1:1572320/1‑150 (MQ=255)
                     ccATGCCCTCGAACAGCTTGCTGTGCACCGCACAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCAc                                                                                                                                 >  1:1347704/1‑150 (MQ=255)
                       aTGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACgg                                                                                                                               >  2:2326466/1‑150 (MQ=255)
                        tGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGa                                                                                                                              >  2:771239/1‑150 (MQ=255)
                           ccTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTGCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGCCGCCACCCAAGATCACGCGACGGATCa                                                                                                                           >  2:361136/1‑150 (MQ=255)
                           ccTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCggg                                                                                                                                                      >  1:953846/1‑123 (MQ=255)
                           ccTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCggg                                                                                                                                                      <  2:953846/123‑1 (MQ=255)
                           ccTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCa                                                                                                                           >  1:968195/1‑150 (MQ=255)
                            cTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAg                                                                                                                          <  2:2184013/150‑1 (MQ=255)
                               gAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCAAGATCACGCCACGGATCAGGtt                                                                                                                       <  2:2144519/150‑1 (MQ=255)
                                aaCAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACgg                                                                                                                               <  1:1085656/141‑1 (MQ=255)
                                aaCAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACgg                                                                                                                               >  2:1085656/1‑141 (MQ=255)
                                  cAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGCAGCTCCAGGCGCAGGGCTTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTgtc                                                                                                                    <  1:311034/150‑1 (MQ=255)
                                    gCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCTGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTgtcgt                                                                                                                  <  1:521855/150‑1 (MQ=255)
                                    gCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCATGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTgtcgt                                                                                                                  >  1:1377974/1‑150 (MQ=255)
                                      ttGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCa                                                                                                                >  2:1758516/1‑150 (MQ=255)
                                       tGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatg                                                                                                                                               <  1:1756367/118‑3 (MQ=255)
                                       tGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatg                                                                                                                                               >  2:1756367/1‑116 (MQ=255)
                                          tgtgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACTCAAGCTCACGCCACTGGTCATGTTGTCGTCCAcgg                                                                                                            >  2:22182/1‑150 (MQ=255)
                                            tgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggcc                                                                                                          >  1:1935605/1‑150 (MQ=255)
                                             gCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggccg                                                                                                         >  2:1493256/1‑150 (MQ=255)
                                             gCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatggt                                                                                                                                             >  2:1510017/1‑110 (MQ=255)
                                             gCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatggt                                                                                                                                             <  1:1510017/114‑5 (MQ=255)
                                              cACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATc                                                                                                                                        >  1:1826676/1‑118 (MQ=255)
                                               aCCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggccggc                                                                                                       >  1:1427437/1‑150 (MQ=255)
                                                     cAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCcg                                                                                                 >  2:821533/1‑150 (MQ=255)
                                                          gcgcTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCAc                                                                                                                                      >  1:1718816/1‑108 (MQ=255)
                                                             cTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCGCCCGTGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGc                                                                                         <  2:1642906/150‑1 (MQ=255)
                                                               ggTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                                       <  1:2326466/150‑1 (MQ=255)
                                                               ggTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                                       <  1:2093852/150‑1 (MQ=255)
                                                               ggTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                                       <  1:1995437/150‑1 (MQ=255)
                                                                gTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAg                                                                                      >  1:79549/1‑150 (MQ=255)
                                                                 tAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGc                                                                                     <  1:1493256/150‑1 (MQ=255)
                                                                   aCGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                                       >  1:830874/1‑146 (MQ=255)
                                                                   aCGCTAGAGGAACTCCAGTCGCAGGGCTTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGAACCGAGCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCAt                                                                                   <  1:762021/150‑1 (MQ=255)
                                                                      cTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGc                                                                                <  1:282883/150‑1 (MQ=255)
                                                                       tCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCATGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGcc                                                                               >  2:477333/1‑150 (MQ=255)
                                                                       tCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGcc                                                                               >  1:1241120/1‑150 (MQ=255)
                                                                        cGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCt                                                                              <  1:2380819/150‑1 (MQ=255)
                                                                        cGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCt                                                                              <  2:453778/150‑1 (MQ=255)
                                                                         gAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTg                                                                             <  2:1332737/150‑1 (MQ=255)
                                                                                tCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTgtc                                                                                                                    >  2:716933/1‑104 (MQ=255)
                                                                                tCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTgtc                                                                                                                    <  1:716933/104‑1 (MQ=255)
                                                                                tCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTg                                                                             <  1:911982/143‑1 (MQ=255)
                                                                                tCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTg                                                                             >  2:911982/1‑143 (MQ=255)
                                                                                tCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATg                                                                      <  2:2204308/150‑1 (MQ=255)
                                                                                  cAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                                       <  2:352520/131‑1 (MQ=255)
                                                                                  cAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                                       >  1:352520/1‑131 (MQ=255)
                                                                                  cAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGc                                                                    >  1:1512805/1‑150 (MQ=255)
                                                                                   aGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCg                                                                   >  1:380633/1‑150 (MQ=255)
                                                                                    ggCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGa                                                                  <  1:1019932/150‑1 (MQ=255)
                                                                                        cAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGaa                                                              <  2:1313723/150‑1 (MQ=255)
                                                                                             ccTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatgg                                                                                                                                              <  2:1464874/65‑4 (MQ=255)
                                                                                              cTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCATGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCgg                                                        <  1:477333/150‑1 (MQ=255)
                                                                                              cTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCgg                                                        <  1:670839/150‑1 (MQ=255)
                                                                                               tCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGACCGCAGCAACCGGTCTGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGgc                                                       >  1:160779/1‑150 (MQ=255)
                                                                                                 gAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGgcgc                                                     >  1:2355038/1‑150 (MQ=255)
                                                                                                  aaGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCaaa                                                                        >  2:2295981/1‑130 (MQ=255)
                                                                                                       gTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTc                                               >  2:728694/1‑150 (MQ=255)
                                                                                                            ccTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCa                                          <  2:1935605/150‑1 (MQ=255)
                                                                                                             cTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCaa                                         <  1:2262988/150‑1 (MQ=255)
                                                                                                             cTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCaa                                         >  2:2262988/1‑150 (MQ=255)
                                                                                                             cTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGGCCAGCaa                                         <  2:1601862/150‑1 (MQ=255)
                                                                                                                 gggTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTg                                     <  2:753261/150‑1 (MQ=255)
                                                                                                                   gTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATg                                                                      >  1:2025004/1‑115 (MQ=255)
                                                                                                                   gTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATg                                                                      <  2:2025004/115‑1 (MQ=255)
                                                                                                                       aTGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAgg                                                                                                                         >  1:790138/1‑60 (MQ=255)
                                                                                                                       aTGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAgg                                                                                                                         <  2:790138/60‑1 (MQ=255)
                                                                                                                       aTGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGgc                                                       <  1:1339225/126‑1 (MQ=255)
                                                                                                                       aTGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGgc                                                       >  2:1339225/1‑126 (MQ=255)
                                                                                                                        tGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCa                              >  2:1636398/1‑150 (MQ=255)
                                                                                                                         ggCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCAt                             <  2:1347704/150‑1 (MQ=255)
                                                                                                                         ggCGAACGCAGTGGCAACCATCGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCTAGGTGCCGTCCAt                             <  2:1241120/150‑1 (MQ=255)
                                                                                                                          gCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTgtc                                                                                                                    <  2:1967419/62‑1 (MQ=255)
                                                                                                                            gAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTc                          <  1:728694/150‑1 (MQ=255)
                                                                                                                            gAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTc                          <  1:771239/150‑1 (MQ=255)
                                                                                                                             aaCGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCGGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCg                         >  2:1772141/1‑150 (MQ=255)
                                                                                                                              aCGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCa                              <  1:1495066/144‑1 (MQ=255)
                                                                                                                              aCGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCa                              >  2:1495066/1‑144 (MQ=255)
                                                                                                                                gCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTcc                                                          <  2:2387829/114‑1 (MQ=255)
                                                                                                                                gCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTcc                                                          >  1:2387829/1‑114 (MQ=255)
                                                                                                                                 cAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACTGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAaga                     <  2:261089/150‑1 (MQ=255)
                                                                                                                                  aGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAaga                     <  2:1784431/149‑1 (MQ=255)
                                                                                                                                  aGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAaga                     >  1:1784431/1‑149 (MQ=255)
                                                                                                                                  aGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAAGTGCCGTCCATGTCGAagag                    >  2:514852/1‑150 (MQ=255)
                                                                                                                                      gCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGCCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACt                >  2:2100137/1‑150 (MQ=255)
                                                                                                                                        tcgggtGGCGCCCGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGt                                      >  2:957419/6‑126 (MQ=255)
                                                                                                                                            aTGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGGCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATTAAGTCCGGCGCCGTGCCAGGTGAGGTGCCGTCCATGTCGAAGAGTACTGCTCGc          >  2:803693/1‑150 (MQ=255)
                                                                                                                                                 gcgGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCgtg                                                 <  2:205901/106‑1 (MQ=255)
                                                                                                                                                 gcgGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCgtg                                                 >  1:205901/1‑106 (MQ=255)
                                                                                                                                                    ggCGCCACCCGAGATCACGCCACGGATCAGGTTGGCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGTGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCAt  <  2:1670719/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCAT  >  NC_002947/1968564‑1968861

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