Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 508,220 T→A 5.5% N264I (AAC→ATC)  estP ← esterase EstP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947508,2200TA 5.5% 218.4 / 2.2 73N264I (AAC→ATC) estPesterase EstP
Reads supporting (aligned to +/- strand):  ref base T (36/32);  new base A (2/2);  total (38/35)
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.59e-01

GTAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGT  >  NC_002947/508075‑508366
                                                                                                                                                 |                                                                                                                                                  
gTAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGCCCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCt                                                                                                                                                 >  5:162787/1‑149 (MQ=255)
 tAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTg                                                                                                                                                <  5:257901/149‑1 (MQ=255)
 tAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGGCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTg                                                                                                                                                <  5:318031/149‑1 (MQ=255)
    tCGGCAATCAGCTGCTGACCCGCGATGGGCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTc                                                                                                                                             <  7:112420/149‑1 (MQ=255)
      ggCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGc                                                                                                                                           >  8:231520/1‑149 (MQ=255)
      ggCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGc                                                                                                                                           >  5:151497/1‑149 (MQ=255)
        cAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTg                                                                                                                                         <  7:231520/149‑1 (MQ=255)
           tCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTc                                                                                                                                             <  8:241787/142‑1 (MQ=255)
           tCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTc                                                                                                                                             >  7:241787/1‑142 (MQ=255)
           tCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGc                                                                                                                                       >  5:336698/1‑148 (MQ=255)
             aGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGtt                                                                                                                                                   <  1:254509/134‑1 (MQ=255)
             aGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGtt                                                                                                                                                   >  2:254509/1‑134 (MQ=255)
             aGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAg                                                                                                                                    <  5:89500/149‑1 (MQ=255)
                 gctgACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGcc                                                                                                                                 >  5:332691/1‑148 (MQ=255)
                     aCCCGCGATGGTCGGGTGGACCGAGTTGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAact                                                                                                                            <  7:292139/149‑1 (MQ=255)
                     aCCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAac                                                                                                                             <  5:107171/148‑1 (MQ=255)
                     aCCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAac                                                                                                                             <  6:151497/148‑1 (MQ=255)
                       ccGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGccccagtcag                                                                                                                          >  1:62929/1‑142 (MQ=255)
                        cgcgATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGAACCTGGTCggc                                                                                                                                          <  5:64250/132‑3 (MQ=255)
                        cgcgATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGAACCTGGTCggc                                                                                                                                          >  6:64250/1‑130 (MQ=255)
                           gATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTgg                                                                                                                      >  2:200140/1‑149 (MQ=255)
                           gATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTgg                                                                                                                      <  6:332691/149‑1 (MQ=255)
                                            aGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAac                                                                                                                             <  1:88794/125‑1 (MQ=255)
                                            aGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAac                                                                                                                             >  2:88794/1‑125 (MQ=255)
                                            aGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTgg                                                                                                                      >  1:128425/1‑132 (MQ=255)
                                            aGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTgg                                                                                                                      <  2:128425/132‑1 (MQ=255)
                                             gTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGTTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGc                                                                                                    <  4:24451/149‑1 (MQ=255)
                                                 ttGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTgg                                                                                                                                               >  1:39198/1‑102 (MQ=255)
                                                 ttGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTgg                                                                                                                                               <  2:39198/102‑1 (MQ=255)
                                                   gAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAAc                                                                                               <  6:336698/148‑1 (MQ=255)
                                                    aaCAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAg                                                                                             >  3:326517/1‑149 (MQ=255)
                                                      cagcagTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGc                                                                                            >  3:143413/1‑148 (MQ=255)
                                                        gcagTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCa                                                                                           >  2:203425/1‑147 (MQ=255)
                                                              tGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGAt                                                                                    <  1:200140/148‑1 (MQ=255)
                                                                gTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTc                                                                                                  >  7:18998/1‑132 (MQ=255)
                                                                gTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTc                                                                                                  <  8:18998/132‑1 (MQ=255)
                                                                 tCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGtt                                                                                 >  6:92865/1‑148 (MQ=255)
                                                                   gggTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGTTCAACAGCAAAGGGATGTTCAg                                                                              <  6:277888/149‑1 (MQ=255)
                                                                   gggTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAg                                                                              >  6:234515/1‑149 (MQ=255)
                                                                     gTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTg                                                                            >  7:316084/1‑149 (MQ=255)
                                                                        tGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGaa                                                                         <  6:162787/149‑1 (MQ=255)
                                                                          gCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGaa                                                                                                                              <  1:337247/94‑1 (MQ=255)
                                                                          gCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGaa                                                                                                                              >  2:337247/1‑94 (MQ=255)
                                                                          gCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGTACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGaa                                                                         <  5:92865/147‑1 (MQ=255)
                                                                          gCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATg                                                                       <  7:114427/149‑1 (MQ=255)
                                                                                tgCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTgg                                                                 <  4:143413/149‑1 (MQ=255)
                                                                                  ccGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTggg                                                                >  5:234896/1‑148 (MQ=255)
                                                                                       gATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGa                                                          >  2:176658/1‑149 (MQ=255)
                                                                                       gATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGa                                                          >  2:31260/1‑149 (MQ=255)
                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATc                                                        >  4:15483/1‑149 (MQ=255)
                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATc                                                        >  6:174911/1‑149 (MQ=255)
                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATc                                                        >  4:10776/1‑149 (MQ=255)
                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATc                                                        >  8:101012/1‑149 (MQ=255)
                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGAAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATc                                                        >  7:263812/1‑149 (MQ=255)
                                                                                             gTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGcc                                                    >  6:343656/1‑149 (MQ=255)
                                                                                                 acCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCAGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGCCCAGCGCATCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAg                                                >  8:180166/1‑149 (MQ=255)
                                                                                                   cGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAgct                                              >  4:308442/1‑149 (MQ=255)
                                                                                                     ggTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAgctgg                                            <  4:326517/149‑1 (MQ=255)
                                                                                                      gTTTTCCACGCAGCTATCGACGCTATAGCACGTGCCGACCAGGTTCTGGTCGCAGGCCAGGCCGAACTGACTGGGGATGGACAGCGCCTCGCTCAAAAGCAAACGGATGATCAGTGGAATGATCTGGGCATCGATCTGCCCTAgctgga                                           >  4:5008/1‑148 (MQ=255)
                                                                                                           ccACGCAGCTATCGCCGCTAAAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAAAAGCAAAGGTATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGAt                                       <  8:125298/148‑1 (MQ=255)
                                                                                                            cACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATc                                      <  5:234515/148‑1 (MQ=255)
                                                                                                                 aGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGAc                                >  8:24256/1‑149 (MQ=255)
                                                                                                                   cTATCGCCGCTTTAGCAGGTGCCGTCCAGGGTCTTGTCGCTGGCCAGGCCTAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACATCAATGGGATGTTCAGTGGAATGATCTGTGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTg                              <  7:180166/149‑1 (MQ=255)
                                                                                                                    tATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTg                                                                  >  4:89381/1‑112 (MQ=255)
                                                                                                                    tATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTg                                                                  <  3:89381/112‑1 (MQ=255)
                                                                                                                     aTCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCg                                                                                                     >  7:203705/1‑76 (MQ=255)
                                                                                                                     aTCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCg                                                                                                     <  8:203705/76‑1 (MQ=255)
                                                                                                                         ccgcTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGaa                        <  5:174911/149‑1 (MQ=255)
                                                                                                                           gcTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACg                      <  6:234896/149‑1 (MQ=255)
                                                                                                                           gcTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACg                      >  1:213027/1‑149 (MQ=255)
                                                                                                                                   aGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCg                 >  8:191814/1‑146 (MQ=255)
                                                                                                                                      cagcGACCAGAACCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGaa                                                                         >  2:150381/4‑87 (MQ=255)
                                                                                                                                      cagcGACCAGAACCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGaa                                                                         <  1:150381/84‑1 (MQ=255)
                                                                                                                                               ggTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGt  <  1:203425/149‑1 (MQ=255)
                                                                                                                                               ggTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGt  <  1:176658/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTATAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAATGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGT  >  NC_002947/508075‑508366

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: