Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,271,146 G→C 9.9% R313R (CGC→CGG) ‡ PP_2869 ← FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,271,1460GC 9.9% 210.7 / 8.3 70R313R (CGC→CGG) ‡PP_2869FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base C (3/4);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCAAGCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGATATGGACCGCGCCA  >  NC_002947/3270999‑3271293
                                                                                                                                                   |                                                                                                                                                   
gggCAAGCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCcgc                                                                                                                                                    >  6:267792/1‑149 (MQ=255)
 ggCAAGCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCcgcg                                                                                                                                                   >  5:110451/1‑149 (MQ=255)
  gCAAGCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCgg                                                                                                                                                  >  7:44250/1‑149 (MQ=255)
      gCCACGATGGAAATGGCCGATTCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGGCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCa                                                                                                                                              <  5:222686/149‑1 (MQ=255)
      gCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCCGCCTGGCCGCGGGc                                                                                                                                                >  6:82032/1‑147 (MQ=255)
       ccACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCCGCGGCCAg                                                                                                                                             >  1:345/1‑149 (MQ=255)
         aCGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAgg                                                                                                                                            <  5:371812/148‑1 (MQ=255)
         aCGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAgg                                                                                                                                            >  6:371812/1‑148 (MQ=255)
         aCGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCCGCGGCCAggc                                                                                                                                           <  2:14118/149‑2 (MQ=255)
         aCGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCCGCGGCCAggc                                                                                                                                           <  2:345/149‑2 (MQ=255)
           gATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGAt                                                                                                                                          >  6:21041/1‑148 (MQ=255)
                 aaTGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCg                                                                                                                                   <  8:44250/149‑1 (MQ=255)
                         atgcatgcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTg                                                                                                                           <  5:21041/149‑1 (MQ=255)
                            catgcCTGGTTGCAGGCGATGCAGGCGCGAACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATAGACGGATCGCAGATCAGCGCCCGGGTCATTCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCGCGGCTTCCTggg                                                                                                                         <  1:283656/148‑1 (MQ=255)
                             atgcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGc                                                                                                                       >  4:379576/1‑149 (MQ=255)
                             atgcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGc                                                                                                                       >  8:218092/1‑149 (MQ=255)
                               gcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTgg                                                                                                                          >  2:52665/1‑144 (MQ=255)
                               gcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTgg                                                                                                                          <  1:52665/144‑1 (MQ=255)
                               gcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTg                                                                                                                     >  7:41080/1‑149 (MQ=255)
                               gcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCCGCGGCCAgg                                                                                                                                            <  1:173374/126‑1 (MQ=255)
                               gcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCCGCGGCCAgg                                                                                                                                            >  2:173374/1‑126 (MQ=255)
                                ccTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTg                                                                                                                     <  7:171139/148‑1 (MQ=255)
                                 cTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACTTCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGt                                                                                                                   <  3:333390/149‑1 (MQ=255)
                                  tGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGtt                                                                                                                  >  6:345863/1‑149 (MQ=255)
                                     ttGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGAt                                                                                                               <  5:345863/149‑1 (MQ=255)
                                     ttGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGAt                                                                                                               <  3:217320/149‑1 (MQ=255)
                                      tGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATg                                                                                                              <  3:314518/149‑1 (MQ=255)
                                          ggCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCTGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGcc                                                                                                          >  4:351189/1‑149 (MQ=255)
                                          ggCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCgg                                                                                                                                  >  7:122222/1‑125 (MQ=255)
                                          ggCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCgg                                                                                                                                  <  8:122222/125‑1 (MQ=255)
                                           gCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGccc                                                                                                         >  5:45604/1‑149 (MQ=255)
                                           gCGATGCAGGCGCACACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGccc                                                                                                         >  8:50285/1‑149 (MQ=255)
                                             gATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCgg                                                                                                       >  2:46709/1‑149 (MQ=255)
                                               tGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTg                                                                                                     <  5:267792/149‑1 (MQ=255)
                                                  aGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGAc                                                                                                   >  6:188928/1‑148 (MQ=255)
                                                            gTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGa                                                                                        >  4:102793/1‑149 (MQ=255)
                                                                  gACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGc                                                                                    >  7:17522/1‑147 (MQ=255)
                                                                  gACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCa                                                                                  <  7:50285/149‑1 (MQ=255)
                                                                   aCCTGGCCGCGTTCGGTCTTGGCTGGGATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGGTGATGCGCCCGGGGACGCACAGGGGGATAGCCAg                                                                                 >  8:346951/1‑149 (MQ=255)
                                                                   aCCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATACGCCCGGTGACGAACAGGGGGATAGCCAg                                                                                 >  7:69948/1‑149 (MQ=255)
                                                                     cTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCt                                                                               >  7:348010/1‑149 (MQ=255)
                                                                     cTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCt                                                                               <  3:186850/149‑1 (MQ=255)
                                                                        gCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGtt                                                                            >  5:370771/1‑149 (MQ=255)
                                                                        gCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGtt                                                                            >  4:324597/1‑149 (MQ=255)
                                                                         ccGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTg                                                                           >  4:265255/1‑149 (MQ=255)
                                                                          cgcgTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGACCAGGATCAGCTCGGCTTCCTGGGGCTGGTTTATGTGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGc                                                                          <  2:79144/149‑1 (MQ=255)
                                                                            cgTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCtt                                                                        <  8:41080/149‑1 (MQ=255)
                                                                              ttCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTg                                                                       >  3:106591/1‑148 (MQ=255)
                                                                                cGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGaa                                                                     <  3:324597/148‑1 (MQ=255)
                                                                                cGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGaa                                                                     <  6:45604/148‑1 (MQ=255)
                                                                                   tCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTg                                                                 >  5:199453/1‑149 (MQ=255)
                                                                                     ttGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGc                                                                >  8:294954/1‑148 (MQ=255)
                                                                                             cATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTcgc                                                       <  4:37742/149‑1 (MQ=255)
                                                                                              aTCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTcgcg                                                      >  1:175839/1‑149 (MQ=255)
                                                                                              aTCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGAAGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGctcgac                                                       >  4:195975/1‑143 (MQ=255)
                                                                                                 gACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGc                                                   <  2:175839/149‑1 (MQ=255)
                                                                                                  aCGGATCGCAGATCAGCGCCCGGGTGATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGcc                                                  <  6:110451/149‑1 (MQ=255)
                                                                                                         gCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACg                                           >  1:167204/1‑149 (MQ=255)
                                                                                                          cAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGc                                          >  4:336373/1‑149 (MQ=255)
                                                                                                             aTCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCgg                                       >  6:263619/1‑149 (MQ=255)
                                                                                                                      cGGGTCATGCCGCACACATCAGCCTGGCCGCTGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGc                              <  3:102793/149‑1 (MQ=255)
                                                                                                                            aTGCCGCACACATGAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCgg                        <  7:348039/149‑1 (MQ=255)
                                                                                                                             tGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGtt                                     >  6:88926/1‑135 (MQ=255)
                                                                                                                             tGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGtt                                     <  5:88926/135‑1 (MQ=255)
                                                                                                                              gCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTg                      >  6:359962/1‑149 (MQ=255)
                                                                                                                                  cacacaTCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCAc                  <  1:46709/149‑1 (MQ=255)
                                                                                                                                     acaTCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGCTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGat               <  5:263619/149‑1 (MQ=255)
                                                                                                                                        tgatccTGGCCCGTGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCATGTACGCCGGTTCGATGGCCATCGGTGGCACGATATg            >  6:202267/5‑149 (MQ=255)
                                                                                                                                          atccTGGCCCGCGGACAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTTCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGata              <  7:121647/143‑1 (MQ=255)
                                                                                                                                          aGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGATATGGa          <  7:15576/149‑1 (MQ=255)
                                                                                                                                            ccTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGATATGGAcc        >  1:303717/1‑149 (MQ=255)
                                                                                                                                                gCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGATATGGACcgcg     <  8:17522/148‑1 (MQ=255)
                                                                                                                                                  cgcgGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATACGCCCGGTGACGAACAGGGGGATAGCCAGCGGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTGCGATGGCCATCGGGGGCACGATATGCACCGCGCCa  <  8:69948/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGCAAGCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGATATGGACCGCGCCA  >  NC_002947/3270999‑3271293

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: