Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,718,096 A→T 12.4% M166K (ATG→AAG) ‡ paaG ← 2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,718,0960AT12.4% 194.4 / 14.4 72M166K (ATG→AAG) ‡paaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (6/3);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTT  >  NC_002947/3717948‑3718244
                                                                                                                                                    |                                                                                                                                                    
gcgcTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCa                                                                                                                                                      >  5:174841/1‑149 (MQ=255)
      gATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCTCAGCATGGc                                                                                                                                                  >  1:283026/1‑147 (MQ=255)
          aGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCATGGCCCGGGCCAAGGCGCTGGCCATGCTGGGCGAgcgc                                                                                                                                            >  8:23374/1‑149 (MQ=37)
              tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGTTCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCa                                                                                                                                                >  7:118153/1‑141 (MQ=255)
              tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGTTCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCa                                                                                                                                                <  8:118153/141‑1 (MQ=255)
              tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTGGCCATGCTGGGCGAGCGCCtt                                                                                                                                         >  2:346375/1‑148 (MQ=255)
              tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTGGCCATGCTGGGCGAGCGCCtt                                                                                                                                         <  1:346375/148‑1 (MQ=255)
              tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGATGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTg                                                                                                                                        <  8:196081/149‑1 (MQ=255)
                  gCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGccc                                                                                                                                    >  8:357707/1‑149 (MQ=255)
                     gTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCggg                                                                                                                                 <  6:201469/149‑1 (MQ=255)
                     gTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCggg                                                                                                                                 >  2:156938/1‑149 (MQ=255)
                        ggTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGcc                                                                                                                               <  4:162448/148‑1 (MQ=255)
                         gTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGc                                                                                                                                >  6:205454/1‑146 (MQ=255)
                         gTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGc                                                                                                                                <  5:205454/146‑1 (MQ=255)
                           gggctgggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGTCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCATGGCCCGGGCCATGc                                                                                                                           <  5:55695/149‑1 (MQ=255)
                              ctgggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGa                                                                                                                        <  7:357707/149‑1 (MQ=255)
                               tgggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGAc                                                                                                                       >  1:80918/1‑149 (MQ=255)
                                  gcggcgAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCaa                                                                                                                    <  1:156938/149‑1 (MQ=255)
                                     gcgAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCATGGCCCGGGCCAAGGCGCTGGCCATGCTGGGCGAGCTCCTTGGCGCCGAACAGGCCGAGCaaag                                                                                                                  >  7:132645/1‑146 (MQ=17)
                                         ggTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCgg                                                                                                              >  2:65131/1‑148 (MQ=255)
                                               gggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTg                                                                                                                                        >  8:216674/1‑116 (MQ=255)
                                               gggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTg                                                                                                                                        <  7:216674/116‑1 (MQ=255)
                                                  cgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCaa                                                                                                    <  6:174841/149‑1 (MQ=255)
                                                         gAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgc                                                                                             <  1:65131/149‑1 (MQ=255)
                                                            ggCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgccg                                                                                           >  3:163676/1‑148 (MQ=255)
                                                              cTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAg                                                                                        >  1:166181/1‑149 (MQ=255)
                                                               ttCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGt                                                                                       >  2:273537/1‑149 (MQ=255)
                                                               ttCATCACGAAGCGCTGCATCATCCACCACGCTGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCTCTCGCCCATCATGGCCAGCGCCTTGGcc                                                                                                                                     <  8:186835/103‑1 (MQ=255)
                                                               ttCATCACGAAGCGCTGCATCATCAACCACGCGGTAGATAACCCACCATTGCTCGGCCTGTTCGGCGCCAAGGAGCTCGCCCAGCATGGCCAGCGCCTTGGcc                                                                                                                                     >  7:186835/1‑103 (MQ=255)
                                                                tcatcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCAt                                                                                                                             <  2:364212/110‑1 (MQ=255)
                                                                tcatcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCAt                                                                                                                             >  1:364212/1‑110 (MQ=255)
                                                                 catcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGt                                                                                       <  1:214284/147‑1 (MQ=255)
                                                                 catcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTcc                                                                                     >  3:191371/1‑149 (MQ=255)
                                                                      cGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCAc                                                                                <  4:338594/149‑1 (MQ=255)
                                                                       gCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCAcc                                                                               >  8:4276/1‑149 (MQ=255)
                                                                           cgcTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCTCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGc                                                                           <  2:283026/149‑1 (MQ=255)
                                                                           cgcTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGc                                                                           <  7:258764/149‑1 (MQ=255)
                                                                                catcatCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATc                                                                      <  2:80918/149‑1 (MQ=255)
                                                                                    atcCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                                   <  4:191371/148‑1 (MQ=255)
                                                                                    atcCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                                   <  2:166181/148‑1 (MQ=255)
                                                                                    atcCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                                   <  6:254211/148‑1 (MQ=255)
                                                                                         ccacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAg                                                             <  1:273537/149‑1 (MQ=255)
                                                                                           acGCGGTAGATCAGCCCCCATTGCTCCGCCTGTTCGGCGCCTAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGc                                                           <  5:325872/149‑1 (MQ=255)
                                                                                                gTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGaa                                                              <  2:16376/141‑1 (MQ=255)
                                                                                                gTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGaa                                                              >  1:16376/1‑141 (MQ=255)
                                                                                                         cccccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTgg                                             <  7:230105/149‑1 (MQ=255)
                                                                                                         cccccATTGCTCGGCCTGTTCCGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTgg                                             >  2:107365/1‑149 (MQ=255)
                                                                                                          ccccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCGTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGc                                            <  6:115481/149‑1 (MQ=255)
                                                                                                                     ggCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcgggcgg                                 >  8:120641/1‑149 (MQ=255)
                                                                                                                      gCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGTCCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggc                                <  7:263465/149‑1 (MQ=255)
                                                                                                                      gCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggc                                <  3:186479/149‑1 (MQ=255)
                                                                                                                      gCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggc                                >  8:263465/1‑149 (MQ=255)
                                                                                                                         tGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcgggcg                                  <  5:307647/144‑1 (MQ=255)
                                                                                                                         tGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcgggcg                                  >  6:307647/1‑144 (MQ=255)
                                                                                                                          gTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcgggcgg                                 <  8:313403/144‑1 (MQ=255)
                                                                                                                          gTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggcagat                            >  7:313403/1‑145 (MQ=255)
                                                                                                                           ttCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGc                            <  7:4276/148‑1 (MQ=255)
                                                                                                                              ggCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCAc                          >  1:155095/1‑147 (MQ=255)
                                                                                                                                 gcCAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGt                     >  6:170109/1‑149 (MQ=39)
                                                                                                                                 gcCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggc                                <  7:344537/138‑1 (MQ=255)
                                                                                                                                 gcCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggc                                >  8:344537/1‑138 (MQ=255)
                                                                                                                                   cAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCGCCGAACAGGCCGAGCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAg                       <  7:326480/145‑1 (MQ=25)
                                                                                                                                   cAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCGCCGAACAGGCCGAGCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAg                       >  8:326480/1‑145 (MQ=25)
                                                                                                                                   cAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTgg                                                       <  1:163486/113‑1 (MQ=38)
                                                                                                                                   cAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTgg                                                       >  2:163486/1‑113 (MQ=38)
                                                                                                                                    aaGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCg                   >  8:80167/1‑148 (MQ=255)
                                                                                                                                      ggCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAg                >  1:359975/1‑149 (MQ=255)
                                                                                                                                      ggCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAg                >  2:174061/1‑149 (MQ=255)
                                                                                                                                       gcgcTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAgg               >  8:38983/1‑149 (MQ=255)
                                                                                                                                         gcTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGcc             <  2:359975/149‑1 (MQ=255)
                                                                                                                                         gcTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGc              >  2:331135/1‑148 (MQ=255)
                                                                                                                                             gCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAgg         >  4:117647/1‑149 (MQ=255)
                                                                                                                                                   cATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGt   <  1:174061/149‑1 (MQ=255)
                                                                                                                                                    aTGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGtt  <  7:80167/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAATGTT  >  NC_002947/3717948‑3718244

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: