Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,717,846 C→A 8.7% V249V (GTG→GTT 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,8460CA 8.7% 228.3 / 8.9 77V249V (GTG→GTTpaaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base A (3/4);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

GGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGC  >  NC_002947/3717705‑3717994
                                                                                                                                             |                                                                                                                                                      
ggTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCtt                                                                                                                                                 <  4:344895/149‑1 (MQ=255)
   ggccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcat                                                                                                                                              <  1:294330/149‑3 (MQ=255)
    gccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgcg                                                                                                                                             <  5:166232/149‑1 (MQ=255)
    gccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgcg                                                                                                                                             <  8:302827/149‑1 (MQ=255)
    gccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgcg                                                                                                                                             >  7:198843/1‑149 (MQ=255)
    gccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcccg                                                                                                                                             >  8:384640/1‑149 (MQ=255)
      cggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                                                                                                                            >  3:4349/1‑148 (MQ=255)
       ggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTa                                                                                                                                          <  2:181638/149‑1 (MQ=255)
        gccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGc                                                                                                                                                    >  1:419921/1‑138 (MQ=255)
        gccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGc                                                                                                                                                    <  2:419921/138‑1 (MQ=255)
            gggCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCt                                                                                                                                     >  4:334990/1‑149 (MQ=255)
                cgcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCCGCGCCCCGTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCt                                                                                                                               <  5:117499/149‑1 (MQ=255)
                cgcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCt                                                                                                                               >  4:90784/1‑149 (MQ=255)
                    aCCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTgcgc                                                                                                                           <  8:112184/149‑1 (MQ=255)
                        gggCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTg                                                                                                                       >  4:191277/1‑149 (MQ=255)
                         ggCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTg                                                                                                                       <  4:4349/148‑1 (MQ=255)
                           cGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGc                                                                                                                      >  7:403439/1‑147 (MQ=255)
                             gtgccGGCGCCATGGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAgg                                                                                                                  >  4:415955/4‑149 (MQ=255)
                                ccGGCGCCATGGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGgcgc                                                                                                               >  4:44029/1‑149 (MQ=255)
                                 cGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTa                                                                                                                                          >  7:298879/1‑123 (MQ=255)
                                 cGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTa                                                                                                                                          <  8:298879/123‑1 (MQ=255)
                                 cGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCt                                                                                                              >  4:9951/1‑149 (MQ=255)
                                         aTGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGGGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTCCCAGGCTCTGCAGGAcg                                                                                                      >  8:360830/1‑149 (MQ=255)
                                           ggCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTcgcg                                                                                                    >  3:264203/1‑149 (MQ=255)
                                                  gCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagc                                                                                             >  7:46649/1‑149 (MQ=255)
                                                     ccGATCCCCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagctc                                                                                           >  2:281062/1‑148 (MQ=255)
                                                        aTCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCtccagctccagc                                                                                       >  4:334944/1‑149 (MQ=255)
                                                        aTCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCtccagctccagc                                                                                       >  3:155705/1‑149 (MQ=255)
                                                           aCCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTgc                                                                                    <  1:83187/149‑1 (MQ=255)
                                                            ccGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTgcg                                                                                   >  6:12731/1‑149 (MQ=255)
                                                             cGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTgcgc                                                                                  <  3:415955/149‑1 (MQ=255)
                                                                      cacaTTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCaaaa                                                                           <  3:334990/147‑1 (MQ=255)
                                                                         aTTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCAtt                                                                      <  8:82066/149‑1 (MQ=255)
                                                                         aTTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCAtt                                                                      <  4:264203/149‑1 (MQ=255)
                                                                          ttgcttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTg                                                                     <  8:198843/149‑1 (MQ=255)
                                                                           ttcttgcGAGTGCGGCCATTTTCAGCGCCCCTTGACCTCTGGCGCGCGTGTGTGCATGTATCCGTGCTCTCCCTTGCGGGTTTCCTCCGGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCAtt                                                                      <  1:343649/147‑1 (MQ=255)
                                                                              ttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTcg                                                                                                      <  4:398028/112‑1 (MQ=255)
                                                                              ttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCg                                                                  >  8:353741/1‑148 (MQ=255)
                                                                              ttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCg                                                                  >  8:203346/1‑148 (MQ=255)
                                                                                    gTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTgg                                                           <  4:412755/149‑1 (MQ=255)
                                                                                        ggCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                       <  3:90784/149‑1 (MQ=255)
                                                                                        ggCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                       >  6:238066/1‑149 (MQ=255)
                                                                                        ggCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                       <  5:189874/149‑1 (MQ=255)
                                                                                           cATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAg                                                    >  1:116625/1‑149 (MQ=255)
                                                                                               ttCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTg                                                >  2:8173/1‑149 (MQ=255)
                                                                                                  aGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTgcgc                                             <  3:9951/149‑1 (MQ=255)
                                                                                                   gcgcCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCt                                            >  8:279932/1‑149 (MQ=255)
                                                                                                       cccTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGAt                                        >  1:244835/1‑149 (MQ=255)
                                                                                                         cTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCa                                      <  1:340472/149‑1 (MQ=255)
                                                                                                            gAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAgt                                    <  7:384640/148‑1 (MQ=255)
                                                                                                             aaGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAgtgt                                  <  5:12731/149‑1 (MQ=255)
                                                                                                                gCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAg                               >  7:339160/1‑149 (MQ=255)
                                                                                                                    gcgtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTgg                                                           >  8:192386/1‑117 (MQ=255)
                                                                                                                    gcgtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTgg                                                           <  7:192386/117‑1 (MQ=255)
                                                                                                                     cgtgcgtTTGTTCATGAAGGCGTTCTCGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGt                          <  7:360830/149‑1 (MQ=255)
                                                                                                                      gtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCg                           >  8:22221/1‑147 (MQ=255)
                                                                                                                       tgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAg                        <  3:44029/149‑1 (MQ=255)
                                                                                                                       tgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAg                        <  2:116625/149‑1 (MQ=255)
                                                                                                                         cgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGGTCAGTGTCAGGTCGTAGGt                      <  8:403439/149‑1 (MQ=255)
                                                                                                                          gtTTGTTCATGATGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGt                      <  3:239975/148‑1 (MQ=255)
                                                                                                                          gtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGt                      >  4:239975/1‑148 (MQ=255)
                                                                                                                            ttGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtggg                   >  4:307020/1‑149 (MQ=255)
                                                                                                                                  gcgAAGGCGTGAACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCa                                                                              >  3:167324/3‑84 (MQ=255)
                                                                                                                                  gcgAAGGCGTGAACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCa                                                                              <  4:167324/82‑1 (MQ=255)
                                                                                                                                  aTGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCACCTCCATTGTCGATAGTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctggg              <  1:281062/148‑1 (MQ=255)
                                                                                                                                   cgAAGGCGTGAACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCgcgc                                                                                                   <  6:418682/61‑1 (MQ=255)
                                                                                                                                   cgAAGGCGTGAACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTg                                          <  5:442075/118‑1 (MQ=255)
                                                                                                                                   cgAAGGCGTGAACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTg                                          >  6:442075/2‑119 (MQ=255)
                                                                                                                                   cgAAGGCGTGAACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCAGCCAGGCGCTGCAGGTCgcgc                                                                                                   >  5:418682/2‑62 (MQ=38)
                                                                                                                                      aGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcggc          >  3:218465/1‑148 (MQ=255)
                                                                                                                                      aGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGgcggcg         >  7:362867/1‑149 (MQ=255)
                                                                                                                                       ggCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGa        <  8:46649/149‑1 (MQ=255)
                                                                                                                                          gTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctgg               >  7:112657/1‑139 (MQ=255)
                                                                                                                                          gTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctgg               <  8:112657/139‑1 (MQ=255)
                                                                                                                                          gTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGt     >  7:43690/1‑149 (MQ=255)
                                                                                                                                          gTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGt     >  1:248528/1‑149 (MQ=255)
                                                                                                                                           ttCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTg    >  4:140494/1‑149 (MQ=255)
                                                                                                                                             cACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTggc  >  3:219761/1‑149 (MQ=255)
                                                                                                                                             cACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTgg   <  2:248528/148‑1 (MQ=255)
                                                                                                                                             cACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTAACAACCATTGTCGAAACTGGCAGTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGCTggc  >  5:52827/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
GGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTC‑‑GCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGC  >  NC_002947/3717705‑3717994

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: