Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,271,148 G→C 9.1% R313G (CGC→GGC) ‡ PP_2869 ← FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,271,1480GC 9.1% 196.4 / 8.7 69R313G (CGC→GGC) ‡PP_2869FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base C (2/4);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GGCAAGCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGATATGGACCGCGCCA  >  NC_002947/3271000‑3271293
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                       >  8:218092/1‑149 (MQ=255)
                            atgcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGc                                                                                                                       >  4:379576/1‑149 (MQ=255)
                              gcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTgg                                                                                                                          <  1:52665/144‑1 (MQ=255)
                              gcCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTgg                                                                                                                          >  2:52665/1‑144 (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                                                                                                                                  <  8:122222/125‑1 (MQ=255)
                                         ggCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCgg                                                                                                                                  >  7:122222/1‑125 (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                                                                            >  4:324597/1‑149 (MQ=255)
                                                                       gCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGtt                                                                            >  5:370771/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)
                                                                                                                                                    |                                                                                                                                                 
GGCAAGCCACGATGGAAATGGCCGATGCATGCCTGGTTGCAGGCGATGCAGGCGCGCACGTCTTCGACCTGGCCGCGTTCGGTCTTGGCTGGCATCGACGGATCGCAGATCAGCGCCCGGGTCATGCCGCACACATCAGCCTGGCCGCGGGCCAGGATCAGCTCGGCTTCCTGGGGCTGGTTGATGCGCCCGGTGACGAACAGGGGGATAGCCAGCTGTTGCTTGAAGGTGCCAGCCTCGCGGGCCAGGTACGCCGGTTCGATGGCCATCGGTGGCACGATATGGACCGCGCCA  >  NC_002947/3271000‑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: