Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,872 G→A 6.1% intergenic (‑71/+73) 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,8720GA 6.1% 243.8 / 4.4 82intergenic (‑71/+73)PP_4634/PP_4635alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base G (34/43);  new base A (3/2);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

AGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAG  >  NC_002947/5257726‑5258017
                                                                                                                                                  |                                                                                                                                                 
agcagAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGaa                                                                                                                                                 <  6:75150/149‑1 (MQ=255)
 gcagAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAAc                                                                                                                                                >  5:319524/1‑149 (MQ=255)
   agAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAAccc                                                                                                                                              >  5:289362/1‑149 (MQ=255)
     aaGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGc                                                                                                                                            >  2:267976/1‑149 (MQ=255)
          gccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTa                                                                                                                                       >  4:53465/1‑149 (MQ=255)
                gCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAggg                                                                                                                                  <  1:113004/148‑1 (MQ=255)
                  tGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAgggg                                                                                                                                 >  6:65860/1‑147 (MQ=255)
                    tATCCGTCGGCCGCGGAGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAggggg                                                                                                                                <  1:5646/146‑1 (MQ=255)
                    tATCCGTCGGCCGCGGAGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAggggg                                                                                                                                >  2:5646/1‑146 (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:141304/148‑1 (MQ=255)
                                                     gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                             <  2:146218/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:65860/148‑1 (MQ=255)
                                                             gATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                     <  5:179408/148‑1 (MQ=255)
                                                                   gattcTCATGCCGCATGTCCGGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTCCTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTAcc                                                                              <  8:220023/145‑1 (MQ=255)
                                                                        tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGAt                                                                          >  3:231497/1‑148 (MQ=255)
                                                                          aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                        >  3:240645/1‑148 (MQ=255)
                                                                          aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                        <  6:290106/148‑1 (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                                                                        <  1:182873/123‑1 (MQ=255)
                                                                                                   ccAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                        >  2:182873/1‑123 (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                                   >  1:152427/1‑149 (MQ=255)
                                                                                                              aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                   >  5:282481/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)
                                                                                                                                                  |                                                                                                                                                 
AGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAG  >  NC_002947/5257726‑5258017

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: