Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,717,845 T→G 16.3% N250H (AAC→CAC)  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,717,8450TG16.3% 145.4 / 21.8 61N250H (AAC→CAC) paaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base G (5/5);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAG  >  NC_002947/3717705‑3717989
                                                                                                                                            |                                                                                                                                                
ggTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCtt                                                                                                                                          >  5:86120/1‑149 (MQ=255)
    gccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgcg                                                                                                                                      >  6:112259/1‑149 (MQ=255)
     ccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcatga                                                                                                                                     <  6:115329/149‑5 (MQ=255)
      cggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                                                                                                                     >  6:229864/1‑148 (MQ=255)
      cggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGt                                                                                                                                    <  6:86120/149‑1 (MQ=255)
        gccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCGCGAAGGCGTGAACGCCTTcatgaac                                                                                                                                   <  1:238285/148‑7 (MQ=255)
        gccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCGCGAAGGCGTGAACGCCTTcatgaac                                                                                                                                   >  2:238285/1‑142 (MQ=255)
                 gcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgcg                                                                                                                                      >  1:237072/1‑136 (MQ=255)
                 gcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgcg                                                                                                                                      <  2:237072/136‑1 (MQ=255)
                   gACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTgcg                                                                                                                       >  3:46208/1‑149 (MQ=255)
                            gATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcatgaac                                                                                                                                   >  1:186634/1‑122 (MQ=255)
                            gATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcatgaac                                                                                                                                   <  2:186634/128‑7 (MQ=255)
                             aTACCGGCGCCCATGACGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAgg                                                                                                             <  5:156354/149‑1 (MQ=255)
                                 cGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAgg                                                                                                             <  8:167153/145‑1 (MQ=255)
                                 cGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAgg                                                                                                             >  7:167153/1‑145 (MQ=255)
                                   gcgcCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGc                                                                                                       <  5:112259/149‑1 (MQ=255)
                                       ccATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGt                                                                                                   >  8:287399/1‑149 (MQ=255)
                                       ccATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGt                                                                                                   >  6:309334/1‑149 (MQ=255)
                                         atcgcgcCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTgcgc                                                                                                                      <  1:63976/125‑1 (MQ=255)
                                         aTGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTAGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTgcgc                                                                                                                      >  2:63976/1‑128 (MQ=255)
                                                   cACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCt                                                                                                                              <  3:139468/110‑1 (MQ=255)
                                                   cACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCt                                                                                                                              >  4:139468/1‑110 (MQ=255)
                                                   cACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagct                                                                                       <  3:186810/149‑1 (MQ=255)
                                                    aCCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagctc                                                                                      <  1:316906/149‑1 (MQ=255)
                                                    aCCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagctc                                                                                      >  6:16889/1‑149 (MQ=255)
                                                       gATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGctccagctccag                                                                                   >  7:22220/1‑149 (MQ=255)
                                                             cGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTgcgc                                                                             >  7:221481/1‑149 (MQ=255)
                                                                cACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTc                                                                          <  4:46208/149‑1 (MQ=255)
                                                                         aTTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCAtt                                                                 <  5:16889/149‑1 (MQ=255)
                                                                            gcttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGt                                                               <  7:287399/148‑1 (MQ=255)
                                                                             cttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccag                                                                                         >  1:138615/1‑121 (MQ=255)
                                                                             cttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccag                                                                                         <  2:138615/121‑1 (MQ=255)
                                                                              ttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCAt                                                                  >  7:101567/1‑143 (MQ=255)
                                                                              ttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCAt                                                                  <  8:101567/143‑1 (MQ=255)
                                                                                   aGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTg                                                       <  8:221481/149‑1 (MQ=255)
                                                                                        ggCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                  >  5:256769/1‑149 (MQ=255)
                                                                                         gCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                  >  6:72416/1‑148 (MQ=255)
                                                                                          ccATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                  >  2:205791/1‑147 (MQ=255)
                                                                                          ccATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                  >  2:210318/1‑147 (MQ=255)
                                                                                                     gcCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCtt                                      >  4:63580/1‑148 (MQ=255)
                                                                                                        ccTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGa                                    <  1:295730/147‑1 (MQ=255)
                                                                                                        ccTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGa                                    >  2:295730/1‑147 (MQ=255)
                                                                                                         cTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGAt                                   >  5:26609/1‑147 (MQ=255)
                                                                                                         cTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCa                                 <  1:205791/149‑1 (MQ=255)
                                                                                                         cTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCa                                 <  1:210318/149‑1 (MQ=255)
                                                                                                           tGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAgt                               >  3:285420/1‑149 (MQ=255)
                                                                                                                    gcgtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCg                      >  8:122004/1‑149 (MQ=255)
                                                                                                                      gtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGa                                                >  4:70412/1‑121 (MQ=255)
                                                                                                                      gtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGa                                                <  7:105573/121‑1 (MQ=255)
                                                                                                                      gtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGa                                                <  3:70412/121‑1 (MQ=255)
                                                                                                                         cgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGt                 <  8:22220/149‑1 (MQ=255)
                                                                                                                                ccgcgAAGGCGTGAACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAgt                               <  3:171099/124‑1 (MQ=255)
                                                                                                                                ccgcgAAGGCGTGAACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAgt                               >  4:171099/5‑128 (MQ=255)
                                                                                                                                ccgcgAAGGCGTGAACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctg           <  4:40359/144‑1 (MQ=255)
                                                                                                                                ccgcgAAGGCGTGAACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctg           >  3:40359/5‑148 (MQ=255)
                                                                                                                                  gcgAAGGCGTGAACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctgggc        >  5:153383/3‑149 (MQ=255)
                                                                                                                                   tGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcg       <  5:229864/149‑1 (MQ=255)
                                                                                                                                    gAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcgg      >  6:229709/1‑149 (MQ=255)
                                                                                                                                       ggCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcggc     >  7:39498/1‑147 (MQ=255)
                                                                                                                                       ggCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcggc     >  6:186135/1‑147 (MQ=255)
                                                                                                                                        gCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAg  <  6:26609/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAG  >  NC_002947/3717705‑3717989

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: ‑