Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,718,093 A→T 12.5% L167Q (CTG→CAG) ‡ 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,0930AT12.5% 187.3 / 14.9 72L167Q (CTG→CAG) ‡paaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (6/3);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGG  >  NC_002947/3717948‑3718237
                                                                                                                                                 |                                                                                                                                                
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                                                                                                                                  <  1:346375/148‑1 (MQ=255)
              tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTGGCCATGCTGGGCGAGCGCCtt                                                                                                                                  >  2:346375/1‑148 (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                                                                                                                         <  5:205454/146‑1 (MQ=255)
                         gTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGc                                                                                                                         >  6:205454/1‑146 (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                                                                                                                      >  1:364212/1‑110 (MQ=255)
                                                                tcatcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCAt                                                                                                                      <  2:364212/110‑1 (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                                                            <  6:254211/148‑1 (MQ=255)
                                                                                    atcCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                            <  2:166181/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         >  2:174061/1‑149 (MQ=255)
                                                                                                                                      ggCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAg         >  1:359975/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)
                                                                                                                                                 |                                                                                                                                                
GCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGG  >  NC_002947/3717948‑3718237

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: