Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,699,749 A→G 8.2% D76G (GAC→GGC) ‡ phaA → poly(3‑hydroxyalkanoate) polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,699,7490AG 8.2% 165.4 / 5.1 61D76G (GAC→GGC) ‡phaApoly(3‑hydroxyalkanoate) polymerase
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (3/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

GCCGCAAGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCGGCCAGTTCG  >  NC_002947/5699602‑5699892
                                                                                                                                                   |                                                                                                                                               
gccgcAAGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCgacgac                                                                                                                                                >  5:138471/1‑149 (MQ=255)
gccgcAAGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCgacgac                                                                                                                                                >  6:34124/1‑149 (MQ=255)
gccgcAAGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCgacgac                                                                                                                                                >  3:227262/1‑149 (MQ=255)
   gcAAGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATACAGCCTGGCCCCGGACAGCGACGACCGt                                                                                                                                             >  5:115650/1‑149 (MQ=255)
     aaGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTc                                                                                                                                            <  4:233293/148‑1 (MQ=255)
      aGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGc                                                                                                                                          <  4:354944/149‑1 (MQ=255)
            tGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAAt                                                                                                                                    >  5:208017/1‑149 (MQ=255)
             gTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGAGAGCGACGACCGTCGCTTCAAt                                                                                                                                    >  6:383044/1‑148 (MQ=255)
              ttGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGGTCGTCGCtgtccgg                                                                                                                                   <  7:160683/148‑7 (MQ=255)
              ttGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGGTCGTCGCtgtccg                                                                                                                                    <  8:260801/147‑6 (MQ=255)
                gAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAAt                                                                                                                                    >  5:379964/1‑145 (MQ=255)
                gAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAAt                                                                                                                                    <  6:379964/145‑1 (MQ=255)
                    tCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGc                                                                                                                            >  1:110311/1‑149 (MQ=255)
                      ggCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGc                                                                                                                            <  8:30191/147‑1 (MQ=255)
                         acgcacCTGGCTGCGCCAGGCCGTGCGCCAAGCGCTGCACAGCGCCAAGCTTGTGGCTCACTATGGCCTGGAGCTGAAGTACGTGTTGCTGGGCAAATCCAGGCTGGCCCCGGACAGCGACGACCCTCGCTTCAATGACCCGGCCTGGa                                                                                                                       <  8:349622/149‑1 (MQ=255)
                            cacCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGca                                                                                                                    <  5:174492/149‑1 (MQ=255)
                             acCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGcaa                                                                                                                   >  3:181171/1‑149 (MQ=255)
                              ccGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGcaac                                                                                                                  <  2:353256/149‑1 (MQ=255)
                               cGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGcaaca                                                                                                                 <  6:208017/149‑1 (MQ=255)
                                  gctgcGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAAcc                                                                                                              >  3:21233/1‑149 (MQ=255)
                                    tgcGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGCCCCGGCCTGGAGCTACACCCCg                                                                                                            >  5:362932/1‑149 (MQ=255)
                                      cgcCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCt                                                                                                          >  1:351608/1‑149 (MQ=255)
                                        ccAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTg                                                                                                         >  8:58419/1‑148 (MQ=255)
                                                     cAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCt                                                                                           >  4:257555/1‑149 (MQ=255)
                                                      aaCCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTg                                                                                          <  5:34124/149‑1 (MQ=255)
                                                      aaCCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTg                                                                                          >  7:62890/1‑149 (MQ=255)
                                                          gCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCaaa                                                                                      >  1:83205/1‑149 (MQ=255)
                                                           cTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAac                                                                                     <  6:138471/149‑1 (MQ=255)
                                                               acaGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCTGCTACCTGCAAacctac                                                                                 >  3:85047/1‑149 (MQ=255)
                                                                      cAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTAcctggcct                                                                          >  5:162124/1‑149 (MQ=255)
                                                                       aaGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTAcctggcctg                                                                         <  8:300870/149‑1 (MQ=255)
                                                                        aGCATTTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCATCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTAcctggcctgg                                                                        <  7:56287/149‑1 (MQ=255)
                                                                             gtgGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGgcgc                                                                     >  4:28278/1‑147 (MQ=255)
                                                                              tgGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAg                                                                  <  4:231863/149‑1 (MQ=255)
                                                                                          gCCTGGAGCTGAAGAACGTGTTGCTGTGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGAc                                                      <  2:351608/149‑1 (MQ=255)
                                                                                             tGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTgg                                                   <  4:21233/149‑1 (MQ=255)
                                                                                               gAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTgg                                                   <  1:90719/147‑1 (MQ=255)
                                                                                                 gCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGa                                               <  5:354304/149‑1 (MQ=255)
                                                                                                   tGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGagc                                             <  7:58419/149‑1 (MQ=255)
                                                                                                    gaagaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGag                                              >  8:106468/1‑147 (MQ=255)
                                                                                                    gaagaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAAAGACCCGGCCTGGAGCAACAACCCGCTGTACCACAGCTACCTGCAAACCTACCTGGCCTGCCGCAAGGAGCTGCAGGACTCGGTGag                                              >  1:144970/1‑147 (MQ=255)
                                                                                                     aagaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGAGAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGagc                                             <  5:383044/147‑1 (MQ=255)
                                                                                                     aagaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGAGGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGag                                              <  6:115650/146‑1 (MQ=255)
                                                                                                     aagaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGagcag                                           >  4:94506/1‑149 (MQ=255)
                                                                                                     aagaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGCCCTGGCGCAACGAGCTGCAGGACTGGGTGagcag                                           >  8:212976/1‑149 (MQ=255)
                                                                                                       gaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTAcctggcctgg                                                                        <  3:335752/118‑1 (MQ=255)
                                                                                                       gaaCGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTAcctggcctgg                                                                        >  4:335752/1‑118 (MQ=255)
                                                                                                             gtTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTggg                                                  <  5:58803/134‑1 (MQ=255)
                                                                                                             gtTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTggg                                                  >  6:58803/1‑134 (MQ=255)
                                                                                                                    ggCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGTGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACGTGGCCTGGCTCAAGGATCTTCAGGACTGTGTGAGCATCAGCGACCTGTcc                              <  6:25208/147‑1 (MQ=255)
                                                                                                                    ggCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTcccc                            >  2:197721/1‑149 (MQ=255)
                                                                                                                        aaTCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAg                        <  4:227262/149‑1 (MQ=255)
                                                                                                                          tCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGt                                >  4:366049/1‑139 (MQ=255)
                                                                                                                          tCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGt                                <  3:366049/139‑1 (MQ=255)
                                                                                                                              gCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCATCGACCTGTCCCCACAGGACAt                   >  3:291755/1‑148 (MQ=255)
                                                                                                                                  ggCCCCGGACAGCGACGACCGACGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTTTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGcc              >  8:121308/1‑149 (MQ=255)
                                                                                                                                    ccccGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCcgc            >  7:271739/1‑149 (MQ=255)
                                                                                                                                    ccccGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCcgc            >  4:368952/1‑149 (MQ=255)
                                                                                                                                     cccGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCcgcg           >  2:232623/1‑149 (MQ=255)
                                                                                                                                      ccGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGTAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCgg          <  2:110311/149‑1 (MQ=255)
                                                                                                                                          gaaGCGACGGTCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCGGCCAg      >  8:240364/3‑149 (MQ=255)
                                                                                                                                           aaGCGACGGTCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGAGCAAGGATCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCACCACAGGACATCAGCCGCGCCCAGt     >  2:213818/2‑149 (MQ=255)
                                                                                                                                           aaGCGACGGTCCTCGCTTCAATGACCCGGCCTGCAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTCGCGCAAGGAGCGTCAGGACTCGCTGAGCAGAAGCGACCTGTACCCCCAGGCAATCAGCCGCGGCCAg      >  7:279391/2‑148 (MQ=255)
                                                                                                                                            aGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCGGCCAGtt    <  3:94506/149‑1 (MQ=255)
                                                                                                                                              cgacggTCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCGGCCAGTTCg  <  7:240364/149‑1 (MQ=255)
                                                                                                                                              cgacgaCCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCGGCCAGTTCg  >  7:125307/1‑149 (MQ=255)
                                                                                                                                              cgacgaCCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTCCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCGGCCAGTTCg  <  3:257555/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCCGCAAGGACCTGTTGAGCTCGGCACGCACCGTGCTGCGCCAGGCCGTGCGCCAACCGCTGCACAGCGCCAAGCATGTGGCTCACTTTGGCCTGGAGCTGAAGAACGTGTTGCTGGGCAAATCCAGCCTGGCCCCGGACAGCGACGACCGTCGCTTCAATGACCCGGCCTGGAGCAACAACCCGCTGTACCGCCGCTACCTGCAAACCTACCTGGCCTGGCGCAAGGAGCTGCAGGACTGGGTGAGCAGCAGCGACCTGTCCCCCCAGGACATCAGCCGCGGCCAGTTCG  >  NC_002947/5699602‑5699892

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: