Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,883 T→C 5.8% intergenic (‑82/+62) PP_4634 ← / ← PP_4635 alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,257,8830TC 5.8% 264.0 / 4.5 87intergenic (‑82/+62)PP_4634/PP_4635alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (34/47);  new base C (3/2);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAG  >  NC_002947/5257736‑5258031
                                                                                                                                                   |                                                                                                                                                     
gccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTa                                                                                                                                                      >  4:53465/1‑149 (MQ=255)
      gCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAggg                                                                                                                                                 <  1:113004/148‑1 (MQ=255)
        tGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAgggg                                                                                                                                                >  6:65860/1‑147 (MQ=255)
          tATCCGTCGGCCGCGGAGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAggggg                                                                                                                                               >  2:5646/1‑146 (MQ=255)
          tATCCGTCGGCCGCGGAGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAggggg                                                                                                                                               <  1:5646/146‑1 (MQ=255)
            tCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAAc                                                                                                                                          <  3:53465/149‑1 (MQ=255)
               gTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCtg                                                                                                                                       <  1:248041/149‑1 (MQ=255)
                  ggCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGGGAACCtgtg                                                                                                                                     <  1:258970/148‑1 (MQ=255)
                  ggCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCtgtg                                                                                                                                     <  5:159537/148‑1 (MQ=255)
                  ggCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTgg                                                                                                                                    <  7:57018/149‑1 (MQ=255)
                    ccGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGa                                                                                                                                  <  1:267976/149‑1 (MQ=255)
                    ccGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCGGTgg                                                                                                                                    >  8:248588/1‑147 (MQ=255)
                    ccGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCACTCCTACAGGGGGAGAACCTGTgg                                                                                                                                    >  5:215703/1‑147 (MQ=255)
                     agcgGGGTTGCGAAATGGGGTAGCTGTTTTGATGATGCTACTAATGGCTCATGTCATGTTGTTGTGCGCAAAAACTAGAAAAGATTTTGTAGCCTTTACTGGCCTCTTCGCGGGTGAACCCTCTCCTACAGGGGGAGAACCTGTGGGa                                                                                                                                  <  8:267853/147‑1 (MQ=37)
                        ggggggGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCgg                                                                                                                              <  5:50050/149‑1 (MQ=255)
                        ggggTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCgg                                                                                                                              <  6:48974/149‑1 (MQ=255)
                        ggggTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCgg                                                                                                                              <  2:319478/149‑1 (MQ=255)
                              gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCgggttcac                                                                                                                        <  6:319524/149‑6 (MQ=255)
                              gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCATGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGc                                                                                                                        >  1:232616/1‑149 (MQ=255)
                                gAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGccc                                                                                                                      >  8:179829/1‑149 (MQ=255)
                                      gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAag                                                                                                                <  8:179770/149‑1 (MQ=255)
                                      gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCACTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAag                                                                                                                <  6:215703/149‑1 (MQ=255)
                                       ggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                               <  4:97095/149‑1 (MQ=255)
                                        gTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                              <  5:314512/149‑1 (MQ=255)
                                        gTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAagag                                                                                                              <  6:241977/149‑1 (MQ=17)
                                           gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGcc                                                                                                           >  8:291227/1‑149 (MQ=255)
                                           gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                            <  2:146218/148‑1 (MQ=255)
                                           gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                            <  2:141304/148‑1 (MQ=255)
                                            cGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                          <  6:289362/149‑1 (MQ=255)
                                            cGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGGTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGGGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                          <  5:327433/149‑1 (MQ=255)
                                               gTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTa                                                                                                       <  7:179829/149‑1 (MQ=255)
                                                   gATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAgg                                                                                                   <  8:87578/149‑1 (MQ=255)
                                                   gATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                    <  5:179408/148‑1 (MQ=255)
                                                   gATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                    <  5:65860/148‑1 (MQ=255)
                                                         gattcTCATGCCGCATGTCCGGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTCCTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTAcc                                                                                             <  8:220023/145‑1 (MQ=255)
                                                              tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGAt                                                                                         >  3:231497/1‑148 (MQ=255)
                                                                aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                       <  6:290106/148‑1 (MQ=255)
                                                                aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                       >  3:240645/1‑148 (MQ=255)
                                                                  ggCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAg                                                                                    >  2:97051/1‑149 (MQ=255)
                                                                   gcgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAg                                                                                    <  7:291227/148‑1 (MQ=255)
                                                                       aTGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                    >  2:175668/1‑128 (MQ=255)
                                                                       aTGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                    <  1:175668/128‑1 (MQ=255)
                                                                            cTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                    >  6:18847/1‑123 (MQ=255)
                                                                            cTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                    <  5:18847/123‑1 (MQ=255)
                                                                            cTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCAGATTACAGCTCGAAGCAg                                                                          >  3:200791/1‑149 (MQ=255)
                                                                             tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGGGAACCTGTGGGAGCGGGCTTGCCCGCCAAGAGGCCGGTACCGGCTTACCCGATTACAGCAAGCCGCAGt                                                                         >  6:216463/1‑149 (MQ=38)
                                                                                  gtgtgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGTTGAACCCGCTCCTACTGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAAc                                                                    <  3:274627/148‑1 (MQ=255)
                                                                                     tgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACtt                                                                  >  7:233708/1‑148 (MQ=255)
                                                                                     tgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTc                                                                 <  7:145593/149‑1 (MQ=255)
                                                                                         ccAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                       >  2:182873/1‑123 (MQ=255)
                                                                                         ccAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                       <  1:182873/123‑1 (MQ=255)
                                                                                          cAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATc                                                            >  7:267364/1‑149 (MQ=255)
                                                                                          cAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACACGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATc                                                            >  5:256910/1‑149 (MQ=255)
                                                                                             aaCTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGt                                                         <  1:74148/149‑1 (MQ=255)
                                                                                               cTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCAACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCCCGAAGAGGCCGGTACAGGCATACCCGATTACAGCACGACGCAGTCGAACGTCACATCCGTgg                                                      >  2:317940/1‑149 (MQ=38)
                                                                                               cTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCCCTGCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCATAACCCAATACAGCACGACGCAGTCGAACTTCACATCGGTgg                                                      >  2:313119/1‑149 (MQ=38)
                                                                                                    aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                                 >  5:282481/1‑149 (MQ=255)
                                                                                                    aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                                 >  1:152427/1‑149 (MQ=255)
                                                                                                    aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTACTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCGGCCAAGAGGCCGGTACAGGCTTACCCGATTTCATCTCGACGCCGTCGAACTTCACATCCGAGGCCACg                                                 >  6:188177/1‑149 (MQ=38)
                                                                                                         ttttGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCttt                                            <  4:240645/149‑1 (MQ=255)
                                                                                                            tGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCGTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGt                                          <  3:132254/148‑1 (MQ=255)
                                                                                                                 ccTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGt                                    <  2:232616/149‑1 (MQ=255)
                                                                                                                       cTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTAcc                                                                                             <  3:195616/87‑1 (MQ=255)
                                                                                                                       cTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTAcc                                                                                             >  4:195616/1‑87 (MQ=255)
                                                                                                                        tGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGt                              >  5:120608/1‑148 (MQ=255)
                                                                                                                          gCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAg                           >  2:44613/1‑149 (MQ=255)
                                                                                                                          gCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAg                           >  1:302821/1‑149 (MQ=255)
                                                                                                                            ctctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcg                         >  2:211654/1‑149 (MQ=17)
                                                                                                                             tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgc                        >  2:208189/1‑149 (MQ=17)
                                                                                                                             tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGCACTTCACATCGGTGGCCACGTCTTTGTCg                                        >  5:299639/1‑133 (MQ=17)
                                                                                                                             tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCg                                        <  6:299639/133‑1 (MQ=18)
                                                                                                                                  gcgGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGGCGTAGTCAACGTCAGCGCGCTc                    >  7:159197/1‑148 (MQ=255)
                                                                                                                                    gggTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGaa                 <  8:233708/149‑1 (MQ=255)
                                                                                                                                     tgTGAACCCGCTCCTACAGGGGGATAACCTGTGGGAGCGGGCTTGCCGGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAg                <  6:256910/148‑1 (MQ=255)
                                                                                                                                     ggTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTc                    >  4:227392/1‑145 (MQ=255)
                                                                                                                                     ggTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTc                    <  3:227392/145‑1 (MQ=255)
                                                                                                                                      gTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTc                                                                 >  7:264320/1‑100 (MQ=255)
                                                                                                                                      gTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTc                                                                 <  8:264320/100‑1 (MQ=255)
                                                                                                                                         aaCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGc               >  7:161327/1‑146 (MQ=255)
                                                                                                                                         aaCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGc               <  8:161327/146‑1 (MQ=255)
                                                                                                                                         aaCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCg             <  4:231497/148‑1 (MQ=255)
                                                                                                                                         aaCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCg             >  2:79879/1‑148 (MQ=255)
                                                                                                                                           cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGAGGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAAc          <  1:97051/149‑1 (MQ=255)
                                                                                                                                           cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc              <  8:159732/145‑1 (MQ=255)
                                                                                                                                           cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc              >  7:159732/1‑145 (MQ=255)
                                                                                                                                                tCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCAGTCAGCCGAACAGGtt     >  1:111613/1‑149 (MQ=255)
                                                                                                                                                  cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGCGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTg    <  1:79879/149‑1 (MQ=255)
                                                                                                                                                   tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg  >  2:152869/1‑149 (MQ=255)
                                                                                                                                                   tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg  >  4:116859/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAG  >  NC_002947/5257736‑5258031

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: