Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,593,940 G→C 7.6% R311R (CGG→CGC PP_3171 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,593,9400GC 7.6% 242.9 / 8.5 79R311R (CGG→CGCPP_3171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base C (3/3);  total (36/43)
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.49e-01

CACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCGGT  >  NC_002947/3593794‑3594086
                                                                                                                                                  |                                                                                                                                                  
cACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAg                                                                                                                                                 >  2:523139/1‑150 (MQ=255)
         gcCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGCCCCGGGGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                        <  1:87742/150‑1 (MQ=255)
         gcCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                        <  2:870827/150‑1 (MQ=255)
         gcCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                        >  2:924222/1‑150 (MQ=255)
                 cGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGc                                                                                                                                <  2:1178679/150‑1 (MQ=255)
                  gATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGcc                                                                                                                               <  1:577834/150‑1 (MQ=255)
                    tCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCcgc                                                                                                                             >  2:298275/1‑150 (MQ=255)
                        gTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAg                                                                                                                         >  2:89658/1‑150 (MQ=255)
                            tGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGAAGCGGTCACAGGCCTATGCCGCTCAGgcgc                                                                                                                     >  2:1017033/1‑150 (MQ=255)
                             ggCCTGCTCCGTAGCCGGAGGGCCAGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgct                                                                                                                    <  2:233512/150‑1 (MQ=255)
                              gCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGAACGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctg                                                                                                                   <  1:990237/150‑1 (MQ=255)
                               ccTGCTCCGTAGCCGGGGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctgg                                                                                                                  >  1:409363/1‑150 (MQ=255)
                                 tGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctggct                                                                                                                >  2:924620/1‑150 (MQ=255)
                                    tCCGTAGCCGGAGTGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGaa                                                                                                             <  1:291460/150‑1 (MQ=255)
                                    tCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGaa                                                                                                             >  2:1421809/1‑150 (MQ=255)
                                       gTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTg                                                                                                          <  2:50056/150‑1 (MQ=255)
                                        tAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACTCGCGGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTgg                                                                                                         <  2:1058675/150‑1 (MQ=255)
                                          gCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGc                                                                                                       <  1:1024508/150‑1 (MQ=255)
                                          gCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCGGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGc                                                                                                       <  1:612190/150‑1 (MQ=255)
                                           ccGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCt                                                                                                      <  1:131923/150‑1 (MQ=255)
                                           ccGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                              <  1:1042558/110‑1 (MQ=255)
                                           ccGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                              >  2:1042558/1‑110 (MQ=255)
                                            cGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTg                                                                                                     >  1:1420394/1‑150 (MQ=255)
                                              gAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGac                                                                                                   <  2:105482/150‑1 (MQ=255)
                                               aGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGaca                                                                                                  <  2:1063658/150‑1 (MQ=255)
                                                    ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                             <  1:1146295/150‑1 (MQ=255)
                                                    ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                             >  2:693399/1‑150 (MQ=255)
                                                    ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                             >  1:445449/1‑150 (MQ=255)
                                                    ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                              >  1:151136/1‑101 (MQ=255)
                                                    ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                              <  2:151136/101‑1 (MQ=255)
                                                      gAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATGGa                                                                                           <  1:1421809/150‑1 (MQ=255)
                                                      gAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                           <  2:1005773/150‑1 (MQ=255)
                                                            gTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGAGACGGTCGAGTaggg                                                                                     >  2:865138/1‑149 (MQ=255)
                                                            gTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTagga                                                                                     >  2:1390754/1‑150 (MQ=255)
                                                              gAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTaggagg                                                                                   <  1:391891/150‑1 (MQ=255)
                                                               aGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGc                                                                                  <  2:302157/150‑1 (MQ=255)
                                                                 ggAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGa                                                                                               <  1:947262/135‑1 (MQ=255)
                                                                 ggAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGa                                                                                               >  2:947262/1‑135 (MQ=255)
                                                                 ggAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCgg                                                                                >  1:575846/1‑150 (MQ=255)
                                                                 ggAAGGGCTGCGCTCGCTGATCGAGCAAACCGGGGCGGAACAGCGGCTGATAGAGTGCCGCATCCATGCCCCGCTGGCACGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCgg                                                                                >  1:574852/1‑150 (MQ=255)
                                                                  gAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGa                                                                               <  1:1279342/150‑1 (MQ=255)
                                                                       gCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTcac                                                                          <  1:924620/150‑1 (MQ=255)
                                                                         tGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                           <  2:1136853/131‑1 (MQ=255)
                                                                         tGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                           >  1:1136853/1‑131 (MQ=255)
                                                                                gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctggct                                                                                                                >  2:983563/1‑103 (MQ=255)
                                                                                gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctggct                                                                                                                <  1:983563/103‑1 (MQ=255)
                                                                                gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                            >  1:483543/1‑139 (MQ=255)
                                                                                gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                            <  2:483543/139‑1 (MQ=255)
                                                                                                  ggCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGc                                               >  2:563157/1‑150 (MQ=255)
                                                                                                         gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAg                                                                                          >  1:1215124/1‑100 (MQ=255)
                                                                                                         gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAg                                                                                          <  2:1215124/100‑1 (MQ=255)
                                                                                                         gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTcaca                                                                         >  1:1364187/1‑117 (MQ=255)
                                                                                                         gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTcaca                                                                         <  2:1364187/117‑1 (MQ=255)
                                                                                                          aGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGCCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGtt                                       <  1:89658/150‑1 (MQ=255)
                                                                                                             tgatgaTCGACTGCCACATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATGCGCACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCg                                    >  1:727688/1‑150 (MQ=255)
                                                                                                                tgatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTac                                 >  2:597701/1‑150 (MQ=255)
                                                                                                                 gatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGc                                                                                  <  1:957228/100‑1 (MQ=255)
                                                                                                                 gatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGc                                                                                  >  2:957228/1‑100 (MQ=255)
                                                                                                                  atCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                        >  2:489928/1‑45 (MQ=255)
                                                                                                                  atCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                        <  1:489928/45‑1 (MQ=255)
                                                                                                                       cTGCCGCCTCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                             >  2:420499/1‑83 (MQ=255)
                                                                                                                       cTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                             <  1:420499/83‑1 (MQ=255)
                                                                                                                                ccATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTa                 >  2:1343889/1‑150 (MQ=255)
                                                                                                                                 cATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTAc                >  1:1248930/1‑150 (MQ=255)
                                                                                                                                  aTGCCCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGt                                                                                         <  1:1387339/76‑1 (MQ=255)
                                                                                                                                  aTGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGt                                                                                         >  2:1387339/1‑76 (MQ=255)
                                                                                                                                     gaccGCTGCCGCGCCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                            >  1:946187/3‑86 (MQ=255)
                                                                                                                                     gaccGCTGCCGCGCCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                            <  2:946187/84‑1 (MQ=255)
                                                                                                                                     aCCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCgg            <  1:523139/150‑1 (MQ=255)
                                                                                                                                       ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                           >  1:323313/1‑69 (MQ=255)
                                                                                                                                       ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                           <  2:323313/69‑1 (MQ=255)
                                                                                                                                       ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTcac                                                         <  2:792021/103‑1 (MQ=255)
                                                                                                                                       ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTcac                                                         >  1:792021/1‑103 (MQ=255)
                                                                                                                                       ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTa                                  <  2:1099762/126‑1 (MQ=255)
                                                                                                                                       ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTa                                  >  1:1099762/1‑126 (MQ=255)
                                                                                                                                            ggCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGcc     <  1:597701/150‑1 (MQ=255)
                                                                                                                                             gcgcGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCg    <  1:924222/150‑1 (MQ=255)
                                                                                                                                             gcgcGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCg    <  1:1343889/150‑1 (MQ=255)
                                                                                                                                               gcgGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCGGt  <  1:693399/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCGGT  >  NC_002947/3593794‑3594086

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