Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,422,392 A→T 8.4% Y110N (TAC→AAC)  tesB ← acyl‑CoA thioesterase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,422,3920AT 8.4% 274.5 / 9.8 94Y110N (TAC→AAC) tesBacyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand):  ref base A (46/40);  new base T (3/5);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGGCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAGTAGCCG  >  NC_002947/5422246‑5422537
                                                                                                                                                  |                                                                                                                                                 
gATCGGCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACt                                                                                                                                                 >  7:46530/1‑149 (MQ=255)
 aTCGGCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTg                                                                                                                                                >  2:92348/1‑149 (MQ=255)
    ggCTTGGCGCACAGCACCTTGTCGCGGATGTGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGaa                                                                                                                                             <  7:68602/149‑1 (MQ=255)
    ggCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGaa                                                                                                                                             <  6:115753/149‑1 (MQ=255)
    ggCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGaa                                                                                                                                             <  1:40291/149‑1 (MQ=255)
     gCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTgg                                                                                                                                               <  3:334807/146‑1 (MQ=255)
     gCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTgg                                                                                                                                               >  4:334807/1‑146 (MQ=255)
          gcgcACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTTCCAGTacgacga                                                                                                                                       <  5:438817/149‑2 (MQ=255)
          gcgcACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTTCCAGTacgacga                                                                                                                                       <  7:150493/149‑2 (MQ=255)
          gcgcACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTTCCAGTacgacga                                                                                                                                       <  6:329403/149‑2 (MQ=255)
           cgcACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTTCCAGTacgacgag                                                                                                                                      <  4:274458/149‑3 (MQ=255)
               caGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCTAAGCCTTCCTCGTCGTACTGGAACGACGCGCTg                                                                                                                                   >  5:352084/1‑148 (MQ=255)
                aGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTg                                                                                                                                   >  6:317972/1‑147 (MQ=255)
                 gCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAg                                                                                                                                <  1:183245/149‑1 (MQ=255)
                 gCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAg                                                                                                                                <  1:277533/149‑1 (MQ=255)
                  cACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAgg                                                                                                                               >  5:433390/1‑149 (MQ=255)
                     cTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAacgac                                                                                                                                         >  3:335819/1‑136 (MQ=255)
                     cTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGa                                                                                                                            <  6:433390/149‑1 (MQ=255)
                     cTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGa                                                                                                                            <  3:78613/149‑1 (MQ=255)
                     cTTGTCGAGGATGGGCTCTGGCTACTGGTCGTCCATGGCATGGGCCAGTTCGACTTCGGTAGGCTGGTTTTCCGGGGCGACCACGTCTGGCATCTGCGCGTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAacgac                                                                                                                                         <  4:335819/136‑1 (MQ=255)
                         tCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGa                                                                                                                         <  6:352084/148‑1 (MQ=255)
                         tCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCCTCCTCGTCGTAc                                                                                                                                                  >  8:162988/1‑123 (MQ=255)
                           gcgGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATgg                                                                                                                      >  8:280981/1‑149 (MQ=255)
                             ggATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTc                                                                                                                    <  4:164048/149‑1 (MQ=255)
                              gATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAg                                                                                                                                <  3:286238/136‑1 (MQ=255)
                              gATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAg                                                                                                                                >  4:286238/1‑136 (MQ=255)
                              gATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCt                                                                                                                   <  6:232916/149‑1 (MQ=255)
                              gATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCt                                                                                                                   <  4:160244/149‑1 (MQ=255)
                                  cgcgcgGGCAACTGGGCGGCCATGGCATGGGCCAGTTCGACTTCGGTACGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGc                                                                                                                <  6:135971/148‑1 (MQ=255)
                                  cgcTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGc                                                                                                                <  1:153943/148‑1 (MQ=255)
                                   gcTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGccc                                                                                                              >  4:147451/1‑149 (MQ=255)
                                   gcTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGccc                                                                                                              >  7:40831/1‑149 (MQ=255)
                                       gggCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCt                                                                                                          <  5:317972/149‑1 (MQ=255)
                                       gggCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCt                                                                                                          <  2:418049/149‑1 (MQ=255)
                                       gggCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCt                                                                                                          <  1:92348/149‑1 (MQ=255)
                                           aaCTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTggct                                                                                                      >  6:233750/1‑147 (MQ=255)
                                           aaCTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGAt                                                                                                      >  5:284080/1‑149 (MQ=255)
                                           aaCTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGACATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGAt                                                                                                      >  7:329562/1‑149 (MQ=255)
                                             cTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCg                                                                                                    >  5:307904/1‑149 (MQ=255)
                                                      catggcatggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTTCCAGTacgacgag                                                                                                                                      <  3:75450/106‑3 (MQ=255)
                                                      catggcatggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTTCCAGTacgacgag                                                                                                                                      >  4:75450/1‑104 (MQ=255)
                                                      catggcatggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCa                                                                                              >  1:351260/1‑146 (MQ=255)
                                                      catggcatggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCa                                                                                              <  2:351260/146‑1 (MQ=255)
                                                      catggcatggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCAccc                                                                                           <  3:147451/149‑1 (MQ=255)
                                                         ggcatggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACccgcc                                                                                        >  2:110717/1‑149 (MQ=255)
                                                          gcatggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACccgccg                                                                                       >  3:255242/1‑149 (MQ=255)
                                                           catggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACccgccgg                                                                                      <  6:288617/149‑2 (MQ=255)
                                                           catggGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCgt                                                                                      <  1:32417/149‑1 (MQ=255)
                                                                    gTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTTCAGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAg                                                                             >  7:278039/1‑149 (MQ=255)
                                                                       cGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTg                                                                          >  6:346025/1‑149 (MQ=255)
                                                                           ttCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTgccgc                                                                      >  2:278109/1‑149 (MQ=255)
                                                                             cGGTCGGCAGGGTTTCCGGGCCGACCAGGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTgccgc                                                                      <  8:278039/147‑1 (MQ=255)
                                                                              ggTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGGCGCGCTGCAGGTGAGGCGGGTCTGGCCCGTCGGGCTCGCCGTCACCCGCCGTGGAGtgccgctgccgcgct                                                                   >  4:318819/1‑146 (MQ=255)
                                                                                   gCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATcgcg                                                               >  5:49788/1‑148 (MQ=255)
                                                                                   gCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCgcgc                                                              >  8:409401/1‑149 (MQ=255)
                                                                                    cAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCgcgc                                                              <  3:192285/148‑1 (MQ=255)
                                                                                       gTTTTCCGGGCCGACCACGTCAGGCTTCTGCGCGTGGTGTTCGAAGCCTTCCTCGTCGTTCTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTAGTGATGCTGCCGCCATCGCGCAcgg                                                          >  6:347302/1‑148 (MQ=255)
                                                                                       gTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCAcgc                                                          >  3:363896/1‑149 (MQ=255)
                                                                                        ttttCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCAcgcg                                                         >  3:37759/1‑149 (MQ=255)
                                                                                            ccGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTc                                                      >  3:239152/1‑148 (MQ=255)
                                                                                             cGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTa                                                    >  4:419397/1‑149 (MQ=255)
                                                                                               ggCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTAcc                                                  <  1:81782/149‑1 (MQ=255)
                                                                                               ggCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTAcc                                                  <  4:363896/149‑1 (MQ=255)
                                                                                                  cGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAg                                               <  4:239152/149‑1 (MQ=255)
                                                                                                  cGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAg                                               <  8:40831/149‑1 (MQ=255)
                                                                                                     ccACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTac                                            <  7:409401/149‑1 (MQ=255)
                                                                                                      cACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTaca                                           >  3:383354/1‑149 (MQ=255)
                                                                                                            aGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTg                                     >  5:407135/1‑149 (MQ=255)
                                                                                                               cATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGaa                                  <  5:346025/149‑1 (MQ=255)
                                                                                                                   tGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCt                              >  2:294263/1‑149 (MQ=255)
                                                                                                                   tGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCt                              >  4:17270/1‑149 (MQ=255)
                                                                                                                     cgcCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCt                              <  6:307904/147‑1 (MQ=255)
                                                                                                                     cgcCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTg                             <  8:46530/148‑1 (MQ=255)
                                                                                                                      gcCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTgg                            >  5:135384/1‑148 (MQ=255)
                                                                                                                        cTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATGGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCAt                         >  5:375445/1‑149 (MQ=255)
                                                                                                                         tGGTGCTCGAAGCCTTCCTCGTCGTACTGTAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATc                        >  1:254502/1‑149 (MQ=255)
                                                                                                                             gCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGccc                    >  8:2210/1‑149 (MQ=255)
                                                                                                                             gCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGccc                    >  6:333150/1‑149 (MQ=255)
                                                                                                                                     ccTTCCTCGTCGTACTGGAACGAGGAGCTGGAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCGGCCGTGTACTGTAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAg            <  8:196129/149‑1 (MQ=255)
                                                                                                                                       ttCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCAGTGGAAGGCTGGCATCGCCCGGGCGCAgg           >  5:386119/1‑148 (MQ=255)
                                                                                                                                       ttCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCt                              <  6:169785/129‑1 (MQ=255)
                                                                                                                                       ttCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCt                              >  5:169785/1‑129 (MQ=255)
                                                                                                                                       ttCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAgg           >  7:273936/1‑148 (MQ=255)
                                                                                                                                         ccTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAg        >  6:224131/1‑149 (MQ=255)
                                                                                                                                          cTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAg        <  1:110717/148‑1 (MQ=255)
                                                                                                                                           tcgtcgTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAGTa      >  8:268721/1‑149 (MQ=255)
                                                                                                                                             gtcgtACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAGTAGc    >  2:161362/1‑149 (MQ=255)
                                                                                                                                             gtcgtACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTACCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAg            >  8:183612/1‑141 (MQ=255)
                                                                                                                                             gtcgtACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTACCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAg            <  7:183612/141‑1 (MQ=255)
                                                                                                                                              tcgtACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAGTAGcc   <  7:280981/149‑1 (MQ=255)
                                                                                                                                              tcgtACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAGTAGcc   >  2:31818/1‑149 (MQ=255)
                                                                                                                                                gtACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAGTAGCCg  <  6:284080/148‑1 (MQ=255)
                                                                                                                                                 tACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGaa         <  3:147055/140‑1 (MQ=255)
                                                                                                                                                 tACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGaa         >  4:147055/1‑140 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GATCGGCTTGGCGCACAGCACCTTGTCGCGGATGCGCTCGGGCAACTGGTCGGCCATGGCATGGGCCAGTTCGACTTCGGTAGGCAGGTTTTCCGGGCCGACCACGTCAGGCATCTGCGCCTGGTGCTCGAAGCCTTCCTCGTCGTACTGGAACGACGCGCTGCAGGTGAAGATGGTCTGGCCCTTCTGGATCGCCGTCACCCGCCGTGTACTGAAGCTGCCGCCATCGCGCACGCGGTCTACCGAGTACACCACTGGAAGGCTGGCATCGCCCGGGCGCAGGAAGTAGCCG  >  NC_002947/5422246‑5422537

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: