Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,703,671 C→T 6.2% S502S (AGC→AGT phaC‑II → poly(3‑hydroxyalkanoate) polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,703,6710CT 6.2% 197.2 / 2.3 64S502S (AGC→AGTphaC‑IIpoly(3‑hydroxyalkanoate) polymerase
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base T (2/2);  total (34/30)
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.67e-01

CATCACACCATGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCGCTAGGCCCCG  >  NC_002947/5703524‑5703819
                                                                                                                                                   |                                                                                                                                                     
caTCACACCATGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGagca                                                                                                                                                      <  2:172169/149‑1 (MQ=255)
  tCACACCATGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGagcagc                                                                                                                                                    <  2:245100/149‑1 (MQ=255)
      acCATGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCACCCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGAcc                                                                                                                                                >  1:144015/1‑149 (MQ=255)
       ccATGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGAccc                                                                                                                                               <  3:264107/149‑1 (MQ=255)
        cATGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGTAGCGACCCa                                                                                                                                              >  6:51564/1‑149 (MQ=255)
          tGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCAcg                                                                                                                                            >  8:97881/1‑149 (MQ=255)
          tGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCAcg                                                                                                                                            >  8:300289/1‑149 (MQ=255)
            ggATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCAcgc                                                                                                                                           <  5:60850/148‑1 (MQ=255)
                 cGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGTAGCGACCCACGCGCCTgg                                                                                                                                     <  5:51564/149‑1 (MQ=255)
                           cTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGa                                                                                                                                                  >  4:309382/1‑126 (MQ=255)
                           cTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGa                                                                                                                                                  <  3:309382/126‑1 (MQ=255)
                           cTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCacgacg                                                                                                                          <  7:213839/149‑1 (MQ=255)
                             cGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGc                                                                                                                         >  2:101640/1‑148 (MQ=255)
                              ggCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAAAAGCGGGCATATCCAGAGCATCATCAACCAGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGcc                                                                                                                        >  5:260700/1‑148 (MQ=255)
                                   tgctgcTGGATGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGtt                                                                                                                                   >  4:256370/1‑133 (MQ=255)
                                   tgctgcTGGATGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGtt                                                                                                                                   <  3:256370/133‑1 (MQ=255)
                                           ggTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGaa                                                                                                            <  3:343047/147‑1 (MQ=255)
                                                cGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCt                                                                                                     <  6:79357/149‑1 (MQ=255)
                                                         cTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGt                                                                                            >  7:87212/1‑149 (MQ=255)
                                                           tCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCtggtg                                                                                                 <  4:277700/142‑1 (MQ=255)
                                                           tCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCtggtg                                                                                                 >  3:277700/1‑142 (MQ=255)
                                                            cGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGTCAAGCGCAGTGAAGGCAGCTGGTGGCCGtt                                                                                           <  2:144015/147‑1 (MQ=255)
                                                             gTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGcc                                                                                                                                        <  2:219435/102‑1 (MQ=255)
                                                             gTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGcc                                                                                                                                        >  1:219435/1‑102 (MQ=255)
                                                                cTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCt                                                                                     >  7:5918/1‑149 (MQ=255)
                                                                 tGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTg                                                                                    >  8:322013/1‑149 (MQ=255)
                                                                     cAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCACCCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCa                                                                                                                 >  7:275963/1‑116 (MQ=255)
                                                                        cAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTgg                                                                                   >  2:180609/1‑143 (MQ=255)
                                                                        cAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTgg                                                                                   <  1:180609/143‑1 (MQ=255)
                                                                           cGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCATACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCAAGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCa                                                                          >  2:311207/1‑149 (MQ=255)
                                                                              gCATAACCAGACCATCCTCAACCCGCCACGCAACCCCAACGCCTACTACCTTGCCAACCCCAAGCTGACCAGCGACCCACGCGCCTGGTTCCACCGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCAcc                                                                        <  1:101640/149‑1 (MQ=255)
                                                                                atatCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGaa                                                                                                            >  3:78057/1‑110 (MQ=255)
                                                                                atatCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGaa                                                                                                            <  4:78057/110‑1 (MQ=255)
                                                                                atatCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGa                                                                             >  2:177039/1‑141 (MQ=255)
                                                                                atatCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGa                                                                             <  1:177039/141‑1 (MQ=255)
                                                                                  atCCATAGCATCATAAAACTTCAAGGCAACCACAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCAc                                                                    <  6:260700/148‑1 (MQ=255)
                                                                                   tCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGc                                                                  >  1:188200/1‑149 (MQ=255)
                                                                                     cAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTc                                                                >  3:291017/1‑149 (MQ=255)
                                                                                      agagCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTcc                                                               >  6:11625/1‑149 (MQ=255)
                                                                                       gagCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCg                                                              >  8:69131/1‑149 (MQ=255)
                                                                                       gagCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCg                                                              <  2:312717/149‑1 (MQ=255)
                                                                                        agCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCgg                                                             >  3:107213/1‑149 (MQ=255)
                                                                                         gCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGc                                                            <  8:5918/149‑1 (MQ=255)
                                                                                            tcatcaACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGTAGCGACCCACGCGCCTGGTTCCACGACGCCaa                                                                                                                      >  6:348205/1‑88 (MQ=255)
                                                                                            tcatcaACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGTAGCGACCCACGCGCCTGGTTCCACGACGCCaa                                                                                                                      <  5:348205/88‑1 (MQ=255)
                                                                                            tcatcaACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTg                                                         <  2:188200/149‑1 (MQ=255)
                                                                                                 aaCCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTc                                                      >  5:311435/1‑147 (MQ=255)
                                                                                                             aaCCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTAc                                        <  7:97881/149‑1 (MQ=255)
                                                                                                                cccAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTg                                      >  1:203902/1‑148 (MQ=255)
                                                                                                                  cAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGa                                                                             >  4:241265/1‑107 (MQ=255)
                                                                                                                  cAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGa                                                                             <  3:241265/107‑1 (MQ=255)
                                                                                                                      gCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTggg                               <  4:107213/149‑1 (MQ=255)
                                                                                                                              ccTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTg                                 <  3:159278/139‑1 (MQ=255)
                                                                                                                              ccTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTg                                 >  4:159278/1‑139 (MQ=255)
                                                                                                                              ccTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCa                       <  5:11625/149‑1 (MQ=255)
                                                                                                                              ccTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCa                       >  6:76589/1‑149 (MQ=255)
                                                                                                                                tGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACt                     >  4:349384/1‑149 (MQ=255)
                                                                                                                                     aaCCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTg                                                         <  4:223555/108‑1 (MQ=255)
                                                                                                                                        cccAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCAcc             <  6:311435/149‑1 (MQ=255)
                                                                                                                                         ccAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCg            >  2:122592/1‑149 (MQ=255)
                                                                                                                                         ccAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCg            >  1:189923/1‑149 (MQ=255)
                                                                                                                                         ccAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCg            >  8:85716/1‑149 (MQ=255)
                                                                                                                                          cAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCGc           <  7:300289/149‑1 (MQ=255)
                                                                                                                                            aGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCGCTa         <  7:69131/149‑1 (MQ=255)
                                                                                                                                            aGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCGCTa         <  1:122592/149‑1 (MQ=255)
                                                                                                                                             gCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCGCTAg        >  5:324675/1‑149 (MQ=255)
                                                                                                                                                   cagcGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCGCTAGGCCCcg  >  1:18813/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATCACACCATGGGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATCAACCCGCCAGGCAACCCCAAGGCCTACTACCTGGCCAACCCCAAGCTGAGCAGCGACCCACGCGCCTGGTTCCACGACGCCAAGCGCAGTGAAGGCAGCTGGTGGCCGTTGTGGCTGGAGTGGATCACCGCACGCTCCGGCCTGCTCAAGGCACCGCGTACTGAACTGGGCAACGCCACTTACCCACCGCTAGGCCCCG  >  NC_002947/5703524‑5703819

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: