Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,259,880 C→T 6.8% intergenic (‑109/+35) 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,259,8800CT 6.8% 312.7 / 8.5 102intergenic (‑109/+35)PP_4634/PP_4635alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base C (56/39);  new base T (4/3);  total (60/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 102 total aligned reads displayed.

GGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGT  >  NC_002947/5259735‑5260027
                                                                                                                                                 |                                                                                                                                                   
ggCGACATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGc                                                                                                                                                 <  1:3782698/150‑1 (MQ=255)
    acaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCGTGCCcgc                                                                                                                                             <  2:3641828/148‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                          <  2:3202624/150‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                          <  2:13926/150‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                          <  1:2489848/150‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                          <  1:4083285/150‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                          <  1:167834/150‑1 (MQ=255)
        ggggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAag                                                                                                                                         >  1:5023024/1‑150 (MQ=255)
        ggggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAag                                                                                                                                         >  1:5023068/1‑150 (MQ=255)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                        >  1:1784277/1‑150 (MQ=255)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                        <  1:3409969/150‑1 (MQ=255)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                        <  2:1086435/150‑1 (MQ=17)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                        <  2:1348832/150‑1 (MQ=17)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                        <  2:4382773/150‑1 (MQ=17)
          ggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                                                       <  1:747479/150‑1 (MQ=255)
          ggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                                                       <  1:3319393/150‑1 (MQ=255)
          ggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                                                       <  1:3879889/150‑1 (MQ=255)
            tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                     >  2:2336135/1‑150 (MQ=255)
            tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                     <  1:333952/150‑1 (MQ=255)
                      gATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGc                                                                                                                           >  1:3641593/1‑150 (MQ=255)
                         gAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTa                                                                                                                        <  2:3658032/150‑1 (MQ=255)
                              tACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGa                                                                                                                   <  2:1461603/150‑1 (MQ=255)
                                 tCATGGCGCCTGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATGTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCGGGGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                >  2:2422195/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                >  1:3621559/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                >  1:5034913/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                >  2:3224517/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                <  1:3824067/150‑1 (MQ=255)
                                   aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACa                                                                                                              >  1:2043127/1‑150 (MQ=255)
                                       cgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCGGTCCGGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTc                                                                                                          <  1:3349877/150‑1 (MQ=255)
                                       cgcATGTCCGTTTGTTGTGCGCAAAAACTACAAAAGATTTTGTAGCCTGTACTGGCCTCTTCGGGGGTGAACCCGCTCCTACCGGGGTAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTc                                                                                                          <  1:4344492/150‑1 (MQ=38)
                                              ccTGTTTTTGTGCGCAAAAACTAGAAAAGTTTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTGGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAg                                                                                                   <  1:4806616/150‑1 (MQ=38)
                                              ccTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAg                                                                                                   <  2:4336192/150‑1 (MQ=255)
                                                tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTc                                                                                                 >  2:2713691/1‑150 (MQ=255)
                                                tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                           <  2:5119364/108‑3 (MQ=255)
                                                tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                           >  1:5119364/1‑106 (MQ=255)
                                                 gttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCg                                                                                                >  1:1189256/1‑150 (MQ=255)
                                                 gttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGGGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCg                                                                                                >  1:1189261/1‑150 (MQ=255)
                                                    gttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAAc                                                                                             <  2:1566889/150‑1 (MQ=255)
                                                             cAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCg                                                                                    >  1:4978378/1‑150 (MQ=255)
                                                                       aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGt                                                                          >  1:2805371/1‑150 (MQ=255)
                                                                       aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGt                                                                          >  2:3970816/1‑150 (MQ=255)
                                                                       aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGt                                                                          >  2:4198754/1‑150 (MQ=255)
                                                                        aaGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGAACCCGCACCTACAGGGGGAGAACCTGCGGGAGCGGGCTTGCCCGCGCAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTc                                                                         >  1:866859/1‑150 (MQ=255)
                                                                         aGATTTGGTAGCCTGTAGTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGGGGCGACCCTGTGGGCGCGGGCTTGCCCGCGACGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCt                                                                        <  1:654641/150‑1 (MQ=39)
                                                                           aTTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCttt                                                                      <  2:1001024/150‑1 (MQ=255)
                                                                              ttGTAGCCTGTACTGGCCTCTTCGCGGGTGACCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTc                                                                   <  2:4831023/150‑1 (MQ=255)
                                                                                           tGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCa                                                      >  2:371975/1‑150 (MQ=255)
                                                                                             gCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGGGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgc                                                    >  2:4390519/1‑150 (MQ=255)
                                                                                              ccTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcg                                                   >  1:3077469/1‑150 (MQ=255)
                                                                                              ccTCTTCGCGGGTGAACCCGCTCCTACAGCGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcg                                                   >  1:3734721/1‑150 (MQ=255)
                                                                                                tctTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                 >  2:1941199/1‑150 (MQ=255)
                                                                                                tctTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTATGTCGTAGTCAACGTCAgcgcg                                                 >  1:3346310/1‑150 (MQ=255)
                                                                                                tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                 >  2:366698/1‑150 (MQ=17)
                                                                                                tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                 >  2:2924204/1‑150 (MQ=17)
                                                                                                tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                 >  1:108873/1‑150 (MQ=17)
                                                                                                  ttCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCt                                               >  2:121699/1‑150 (MQ=255)
                                                                                                     gcgGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGa                                            >  1:656526/1‑150 (MQ=255)
                                                                                                      cgGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGaa                                           <  2:3638605/150‑1 (MQ=255)
                                                                                                        ggTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGc                                         >  1:1050213/1‑150 (MQ=255)
                                                                                                         gTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc                                        >  2:4019164/1‑150 (MQ=255)
                                                                                                         gTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc                                        >  2:4019482/1‑150 (MQ=255)
                                                                                                             aCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAAc                                    <  2:5023024/150‑1 (MQ=255)
                                                                                                             aCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAAc                                    <  2:5023068/150‑1 (MQ=255)
                                                                                                              cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACa                                   >  1:1313923/1‑150 (MQ=255)
                                                                                                              cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACa                                   >  1:3893932/1‑150 (MQ=255)
                                                                                                               ccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAg                                          <  2:1008500/142‑1 (MQ=255)
                                                                                                               ccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAg                                          >  1:1008500/1‑142 (MQ=255)
                                                                                                                 gCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGt                                >  2:913028/1‑150 (MQ=255)
                                                                                                                   tCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTg                              >  1:717336/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                            >  2:897331/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                            >  1:856864/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                            >  1:856723/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                            >  1:4938550/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                            >  1:4938462/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGTGCGCTCGAAGCCGAACAGGTTGAgg                           >  2:2002544/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                           >  1:1312958/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                           >  1:1313223/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                           >  2:4646865/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                           >  2:4646812/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                           >  2:4145282/1‑150 (MQ=255)
                                                                                                                        cAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGaa                         >  1:1052957/1‑150 (MQ=255)
                                                                                                                          gggggAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACt                       <  2:3445535/150‑1 (MQ=255)
                                                                                                                           ggggAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTg                      >  1:3754823/1‑150 (MQ=255)
                                                                                                                             ggAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGct                    >  1:466542/1‑150 (MQ=255)
                                                                                                                               tcccccTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttc                  >  1:5120712/5‑150 (MQ=38)
                                                                                                                               agaACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGCTGCCAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttc                  <  2:3864998/150‑1 (MQ=255)
                                                                                                                               agaACCTGGGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGGAGTCAACGGCTGCGCGCTCGAAGCCGCCCAGGTTGAGGAACTGcttc                  >  1:325485/1‑150 (MQ=39)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGTAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                <  2:2714981/150‑1 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                >  1:2447783/1‑150 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                <  2:5034913/150‑1 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                <  1:2046787/150‑1 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttcta                <  1:1671638/150‑2 (MQ=255)
                                                                                                                                  aCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGCACTGCTTCTTg               <  2:1207333/150‑1 (MQ=255)
                                                                                                                                     tgtgGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTaa            >  2:2521406/1‑150 (MQ=255)
                                                                                                                                         ggAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCgg        >  1:3310947/1‑150 (MQ=255)
                                                                                                                                             cGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTa    <  2:3621559/150‑1 (MQ=255)
                                                                                                                                             cGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTa    >  2:1923020/1‑150 (MQ=255)
                                                                                                                                               ggCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTagt  <  2:755287/150‑1 (MQ=255)
                                                                                                                                                gCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                                                          >  2:3398174/1‑93 (MQ=255)
                                                                                                                                                gCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                                                          <  1:3398174/93‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGT  >  NC_002947/5259735‑5260027

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