Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,718,096 A→T 11.9% 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,0960AT11.9% 156.5 / 10.6 59M166K (ATG→AAG) ‡paaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base T (5/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAAT  >  NC_002947/3717950‑3718241
                                                                                                                                                  |                                                                                                                                                   
gcTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATg                                                                                                                                                   >  5:56527/1‑149 (MQ=255)
     aTCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAg                                                                                                                                              <  4:99470/149‑1 (MQ=255)
         gtgtCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTACTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGcc                                                                                                                                          <  3:329195/149‑1 (MQ=255)
          tgtCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCt                                                                                                                                         <  6:56527/149‑1 (MQ=255)
                         atcttgggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGc                                                                                                                          <  3:96505/145‑1 (MQ=255)
                            ctgggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGa                                                                                                                       <  1:157065/149‑1 (MQ=255)
                            ctgggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGa                                                                                                                       <  2:31017/149‑1 (MQ=255)
                              gggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGAc                                                                                                                      >  8:201185/1‑148 (MQ=255)
                              gggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGAGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCCTGc                                                                                                                          >  4:96505/1‑144 (MQ=255)
                                             gggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTTGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCagcagc                                                                                                      >  8:155172/1‑149 (MQ=255)
                                              ggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGcc                                                                                                     <  6:308083/149‑1 (MQ=255)
                                              ggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGcc                                                                                                     <  2:225312/149‑1 (MQ=255)
                                                        aGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGcc                                                                                                     <  3:265416/139‑1 (MQ=255)
                                                        aGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGcc                                                                                                     >  4:265416/1‑139 (MQ=255)
                                                        aGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgcc                                                                                         >  1:282538/1‑149 (MQ=255)
                                                            cTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAg                                                                                     >  1:3791/1‑149 (MQ=255)
                                                             ttCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGt                                                                                    >  8:232843/1‑149 (MQ=255)
                                                             ttCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGt                                                                                    >  8:237334/1‑149 (MQ=255)
                                                               catcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTcc                                                                                  >  6:154222/1‑149 (MQ=255)
                                                                atcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCg                                                                                 >  5:153685/1‑149 (MQ=255)
                                                                atcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCg                                                                                 <  7:201185/149‑1 (MQ=255)
                                                                    cGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCAc                                                                             >  2:212172/1‑149 (MQ=255)
                                                                    cGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCAc                                                                             >  7:262334/1‑149 (MQ=255)
                                                                         cgcTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGc                                                                        <  2:282538/149‑1 (MQ=255)
                                                                            tGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGGCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCgt                                                                     <  5:353221/149‑2 (MQ=255)
                                                                              catcatCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGTGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATc                                                                   >  6:217002/1‑149 (MQ=255)
                                                                              catcatCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATc                                                                   <  6:129562/149‑1 (MQ=255)
                                                                                tcatcCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGTGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCtt                                                                 <  5:217002/149‑1 (MQ=255)
                                                                                  atcCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGc                                                               <  1:212172/149‑1 (MQ=255)
                                                                                          cgcgGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGCAGGcc                                                       >  5:19875/1‑149 (MQ=255)
                                                                                             ggTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTgg                                                    <  7:232843/149‑1 (MQ=255)
                                                                                             ggTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTgg                                                    <  7:237334/149‑1 (MQ=255)
                                                                                             ggTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTg                                                     >  5:169118/1‑148 (MQ=255)
                                                                                                   tCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGaaa                                              >  6:246634/1‑149 (MQ=255)
                                                                                                        ccccATTGCTCGGCCTGTTTGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGc                                         <  7:155172/149‑1 (MQ=255)
                                                                                                        ccccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGc                                         <  5:246634/149‑1 (MQ=255)
                                                                                                         cccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCtt                                                                 >  6:107782/1‑124 (MQ=255)
                                                                                                         cccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCtt                                                                 <  5:107782/124‑1 (MQ=255)
                                                                                                         cccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGcc                                        >  2:70323/1‑149 (MQ=255)
                                                                                                         cccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGcc                                        <  6:153685/149‑1 (MQ=255)
                                                                                                         cccATTGCTCGGCCAGATCGGCCCCAAGGCGATCGCCCAGCATGGAAAACGCCTTGGCCCGGGCCATGCCGACCAACCGCGGCAGCAGACAAGTACCGCCGGAGTCCGGCACAAGGACGATATTGCAGAAGGCCAGGATGAAACTGGcc                                        >  8:333596/1‑149 (MQ=255)
                                                                                                             ttGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGGGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgc                                    <  2:353137/149‑1 (MQ=255)
                                                                                                             ttGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgc                                    <  8:262334/149‑1 (MQ=255)
                                                                                                                     ccTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggc                             >  2:62439/1‑148 (MQ=255)
                                                                                                                      atggcccgcgcCAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGCCCAACCGCGGCAGCAGCCAAGTACCGCCGCAGTCCGGCACCACGCCGAAATTGCAGAAGGCCAGGATGAAACTGGCCGAGCGGGCGGCCa                           >  6:109488/8‑149 (MQ=25)
                                                                                                                                 cAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCGCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCg                >  1:328968/1‑149 (MQ=38)
                                                                                                                                 cAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGc                                                               <  4:117634/102‑1 (MQ=37)
                                                                                                                                 cAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGc                                                               >  3:117634/1‑102 (MQ=37)
                                                                                                                                   aGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCAGACCAAGCGCGGCAGCAGCCAAGTAGGAGCACGGAGTCCCGCACCAAGCCGATATTGCAGAAGCCCTGGATGAAACTGGCCGCGCGGGCGGCCAGCACCAAGTc                 >  1:290286/1‑148 (MQ=255)
                                                                                                                                    ggCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCt                                                                  <  5:69791/96‑1 (MQ=37)
                                                                                                                                    ggCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCt                                                                  >  6:69791/1‑96 (MQ=37)
                                                                                                                                     gcgcTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAgg            >  5:247421/1‑149 (MQ=255)
                                                                                                                                          cGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAg       <  7:333596/149‑1 (MQ=255)
                                                                                                                                            cccAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAggg     <  5:154222/149‑1 (MQ=255)
                                                                                                                                               aGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCtt                                                                 >  7:183544/1‑86 (MQ=255)
                                                                                                                                               aGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAg       <  6:275730/144‑1 (MQ=255)
                                                                                                                                               aGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAg       >  5:275730/1‑144 (MQ=255)
                                                                                                                                               aGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAAt  >  7:120670/1‑149 (MQ=255)
                                                                                                                                               aGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGCCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCtt                                                                 <  8:183544/86‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGGAAT  >  NC_002947/3717950‑3718241

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